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Reflection, Transmission, and Mode Conversion 
in a Corrugated Feed 


By C. DRAGONE 


(Manuscript received December 14, 1976) 


Microwave antennas are often required to carry signals simulta- 
neously over a broad range of frequencies—e.g., the combined TD-2 and 
TH common carrier bands encompass a total frequency ratio of about 
1.8 to 1 as do the combined 18- and 30-GHz bands. To achieve these 
bandwidths, an efficient broadband feed horn is required. The corru- 
gated (hybrid-mode) horn is a leading candidate, but it is not immune 
to some cross-polarization coupling, input reflection, and pattern 
asymmetry. These problems are introduced mainly by two phenomena: 
variation of the dominant mode shape with frequency and mode con- 
version along the horn taper and at waveguide transitions at the horn 
input. Simple formulas for computing the magnitude of these phe- 
nomena and their effects on return loss and radiation patterns are 
given. 


1. INTRODUCTION 


Corrugated feeds (also called hybrid-mode feeds) are widely used in 
reflector-type antennas because of their excellent radiation character- 
istics.!-16 At the frequency wo at which the surface reactance X, of the 
corrugations becomes infinite, the radiation pattern of a properly de- 
signed feed is circularly symmetric, is free of cross-polarized components, 
and has low sidelobes. In principle these properties can be obtained over 
a frequency range of more than an octave. In fact, one can show that the 
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Fig. 1—Corrugated horn of Ref. 17. 


field over an aperture illuminated by the fundamental mode has the 
remarkable property that in the limit, as 


kao, (k= 2zx/d), 


where a is the aperture radius, the field distribution becomes indepen- 
dent of the surface reactance X, (provided X, ~ 0). A corollary of this 
behavior is that a feed of sufficiently large aperture will have the above 
radiation characteristics over a wide range of frequencies provided only 
the fundamental mode is excited in the horn. To verify this behavior, 
an experiment, described in a companion article,!? was made. A very long 
horn (see Fig. 1) was fabricated carefully, using a special fabrication 
technique to minimize geometrical imperfections in the corrugated walls, 
and the radiation characteristics were measured from 17 GHz to 35 GHz. 
From 17 GHz to 29 GHz the far field was found, as expected, to be es- 
sentially polarized in one direction. At frequencies above 29 GHz, 
however, a cross-polarized component was found to be caused by a cer- 
tain undesirable mode, which will be called the HE;,-mode. This mode 
was excited primarily at the input, where the corrugated waveguide was 
connected directly to a smooth waveguide, as shown in Fig. 1. A calcu- 
lation, given in eqs. (102) and (103) of this article, was therefore made 
to determine the total amount of power converted from the TE;;-mode, 
incident at the input, to the HE,,-mode. 

A peculiarity of corrugated feeds is that there is some mode conversion 
even in a conical horn of constant surface reactance X, (unless X, = 0 
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or X, = ~). An evaluation of this effect is given in Section VII. However, 
in the experiment, the taper angle a (Fig. 1) was chosen sufficiently small 
(a ~ 4°) so that this effect was negligible. 

The analysis starts in Section III with a derivation of the asymptotic 
properties for large ka of the modes of a corrugated waveguide. The re- 
sults provide a simple and accurate representation of the modes in a feed 
aperture of more than a few wavelengths in diameter. Then, in Sections 
V and VI, a first-order derivation of the scattering parameters of a 
junction between two waveguides of slightly different characteristics 
is given. A simple relation [see, eqs. (83), (84), and (115) to (117)] is found 
between the scattering parameters and the coupling coefficients between 
the modes on the two sides of the junction. Each coefficient is given, 
except for a constant, by an integral of the type 


f fi (E, X Ho*) -i, dxdy, 


where S is the junction area, and E, and Hg are the electric and magnetic 
vectors of the two modes, respectively. In Section IV, this surface integral 
is converted to a line integral, thus reducing the calculation of the cou- 
pling coefficients to a straightforward exercise. This result is useful also 
to calculate the far field of an aperture S illuminated by a mode Ej, since 
the far field at a given observation point is, except for a constant, the 
coupling coefficient over S between EK and the field He of a plane wave 
having the direction of the observation point. The far-field calculation 
is thus reduced to a contour integration. 

The calculation of the scattering parameters is carried out in Sections 
IV to VI, using the above contour integral. It is found, for instance, that 
the input reflection of a corrugated feed connected to a smooth wave- 
guide of the same diameter is simply given by the coefficient 
Cie 

By + By’ 

@, and #; being the propagation constants in the two waveguides. An 
identical formula was derived by Brown!® from the principle of conser- 
vation of momentum, but that derivation is not applicable to the present 
problem, which involves hybrid modes. 

Finally, Section VII deals with the problem of spurious mode gener- 
ation in a nonuniform waveguide whose parameters (radius and surface 
reactance) vary along the axis, as in Fig. 1. The differential scattering 
parameters that give, at any point in a nonuniform waveguide, the local 
coupling between the incident mode and the spurious modes are ob- 
tained from the analysis of Sections V and VI. By solving the differential 
equations specified by the above scattering parameters, we can thus 
determine the amplitudes of the spurious modes. An example is provided 


Pi 
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of a first-order calculation of mode conversion in a conical waveguide 
such as the one in Fig. 1 for z > z,. The result, eq. (154), is again quite 
simple. 


ll. PRELIMINARY CONSIDERATIONS 


For a smooth waveguide, there is a simple relation between the 
propagation constant 8 of a mode and the waveguide diameter, but no 
such simple relation exists in the case of a corrugated waveguide [see eq. 
(20)]. For this reason, the properties of the corrugated waveguide modes 
cannot generally be determined as simply as in the case of a smooth 
waveguide. Also, the field configuration of each mode varies with 
waveguide diameter. There is, however, an important exception. When 
the radius a of the corrugated waveguide is sufficiently large, the prop- 
agation constant 8 for some of the modes is simply given by 


8 = V(ka)? = usm » (1) 
where Ugm is the mth zero of the Bessel function Jo of order zero, 
Jo(Uom) = 0. (2) 


For all the other modes except one (for this special mode @ is independent 
of a; see Appendix B) one has 


B = V(ka)?—u3,, (3) 
where Ug is the mth root of the Bessel function of order two, 
Jo(Uam) = 0. (4) 


Equations (1) and (3) are valid provided a > X, a condition which is 
satisfied to a good approximation by most feed apertures. Thus, the 
case 


ka>1 (5) 


is of considerable practical interest. One finds that as ka — o, the 
properties of a mode become independent of the surface reactance X, 
of the corrugated walls, except for the mode of Appendix B. Thus, the 
field distribution over the aperture of a feed illuminated by a single mode 
will be little affected by the surface reactance X, (which varies fairly 
rapidly with frequency) provided ka is sufficiently large. This result, first 
pointed out by Thomas,? is very important for it implies that the aper- 
ture field distribution becomes frequency independent for large ka. The 
main purpose of this section is to determine the asymptotic behavior of 
the hybrid modes for large ka. It is shown that if ka # ~ there is over 
the aperture of a feed a certain undesirable cross-polarized component, 
even if the aperture is illuminated by a single mode, unless of course X; 
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= «, A simple expression for the amplitude of this component is 
given. 


lil. ASYMPTOTIC BEHAVIOR FOR LARGE ka 

Consider a disk-loaded waveguide centered around the z-axis, as in 
Fig. 1, and assume its parameters a, b, and h are independent of z. Let 
r, ¢, 2 be cylindrical coordinates defined by x = r cosd¢ and y = r sing. 
The separation of the disks, which occupy the region a <r < }, is as- 
sumed to be much smaller than a wavelength X, 


kh «1. (6) 


The region between two consecutive disks forms a radial line whose input 
reactance jX atr =a isa function of the radial length / = b — a; for ka 
> 1, one has approximately 


JX = jZ tankl, 


where Zo = Vuo/eg - For a finite number of teeth per wavelength, the 
value of | must be corrected.* Because of condition (6) the effect of the 
disks can be accounted for adequately by introducing an effective surface 
reactance®:!2,19 


ere Ge) (7) 


where t is the thickness of the disks, and by requiring that the field for 
r <a satisfy the boundary conditions 


forr =a, (8) 


where Ey, Hy, E, are the ¢ and z components of the electric and magnetic 
field. 
Let 8 be the propagation constant in the z direction, 


B =k cosh, (9) 


and assume 6; is real, so that 6 < k. The case where 0 is imaginary is 
considered in Appendix B. Assume the ¢ dependence of E, is given by 
cos¢. Then, the field components of a mode that propagates in the z 
direction with propagation constant 8 are eiven by 


E, = Ad,(xr) cos¢ e~JF2, (10) 


* See Ref. 17 for the effect of a finite number of teeth per wavelength, which causes a 
reduction of the effective depth, !. 
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ZoH, = BJ (xr) sing e7/F2, (11) 
J1(xr) 
Kr 





Ey= a 5 | Bel (xr) + A cos6; 
S 


| sing e—J8z, (12) 
ind, 


J (xr) 
Kr 





ZoHs = — ay | Avr) + B cos6; 


cosd e762, (18) 
inf; 


.=— = aes a ser) + A cos#; Titer) cos@ e ~I8z, (14) 
sin@; Kr 
ZoH, = — eee IP lA dite) —— + B cos6; Js(or) sing e~J62, (15) 
sind; Kr 
(r <a), 
where 
K=k sing). 


The boundary conditions (8) give!? 


oth J J(u) 

= 16 
ss ~ cos; Jy(u)’ oo 

cos#; 1 1 J;(u) 
— ro ~~ + Ses 17 
sind; uy  sinO; J,(u) a 

where 
u = ka sin6y, 
Zo 

=—-— 18 
y=-3 (18) 


and y is the ratio between the TM and TE components of the hybrid 


mode, 


A 
Ve (19) 


By eliminating y from eqs. (16) and (17), one obtains the eigenvalue 


equation 
y _ 1 Jy(u) | 2 u? 
ie ee ey Se a 20 
ka u2ud,(u) any, (ka)2 20) 
which is eq. (10) of Ref. 12. 
The solutions of this equation are now studied for large ka. Both u and 
y are assumed to be finite. Then, in the limit as ka — ~, eq. (20) reduces 
to 





Jiu) \2 °41 
Ga") 1=0. (21) 
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We distinguish two cases: 
Ji(u)u _ 


= 22 
J(u) 22 
and 
Ji(u)u 
ee 23 
J\(u) ve) 


According to eq. (16) (with cos#; ~ 1, since 6; > 0 as ka > ~), these two 
cases correspond respectively to 


ee (24) 
and 
y=-l1. (25) 


Using well known recurrence relations between the Bessel functions and 
their derivatives, and using conditions (22) and (23), we find 


Jo(u) = 0 (26) 
and 
J(u) = 0. (27) 


We conclude that for large ka, eq. (20) possesses the two sets of solu- 
tions 


U ~ Uom (m = 1, 2, etc.) (28) 
and 
U ~ Uo (m = 1, 2, etc.), (29) 


Uom and Uy, being respectively the mth root of Jo(u) and Jo(u). Solu- 

tions (28) and (29) are characterized by y ~ 1 and y ~ —1; the corre- 

sponding modes will be designated,* respectively, HE; and HEj,. 
Asymptotic series for u and y in terms of 


1 
ae (30) 
ka 

* This mode classification differs from the one by Clarricoats!? and it was chosen for the 
following reason. Here, and in Ref. 17, we are interested i in horns whose inner radius a varies 
gradually with z, while the wall susceptance y is approximately constant, as in Fig. 1, from 
zo toZ}. Consider therefore a mode propagating in Fig. 1 from z2 towards 21. Clarricoats’ 
classification assigns in some cases a different name to this mode in different regions of 
the horn, even though there will be no discontinuous variation of the mode-field config- 
uration, as it propagates in the horn. On the other hand, our classification based on the 
Bessel function roots Upm and Ugm, assigns a single name everywhere i in the horn. If instead 
the frequency is gradually changed the mode of a waveguide of given dimensions will retain 
the same name with Clarricoats’ classification, whereas this is not always true with our 
classification. To understand better these considerations, see also Ref. 25. 
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are derived in Appendix A under the assumption y # ~. For the HE,» 
modes, characterized by y — 1 as ka — ~, it is found that 


1 1 2 1 \2 
u = tm = om [1-5-5 [1-5 + wba) | (Z-) 


21-2 ena] (2) ] 0 
and 


Yaak. aye 1\2 
=] -: JOP ES (Ae ae Ge 
t = | 8 ( Wom) ka 


2 ka 
Sls biel ae IG] 2 
ae 5 Hom q (Bt + 2) re . (82) 


For the HE,,,-modes, characterized by y ~ —1, u is given by 


; ly 
7 = Uom {1 -=— 
u Uo | cs | (33) 
and 
ly 
=-—-|—-— Bo. —~— eo 8 ; 
. Us Ee | (34) 


The x and y components of the electric field are now derived. First 
consider the HE,,, modes. One finds from eqs. (10) to (15), with cos6; 
= 1 and y given by eq. (82), that for large ka the transverse component 
of E is given by 


B= -j A Jo (cu )is 
u a 
1 


Pe ee ayo ; re | 
4 u ra Jo ( u) (cos2¢ i, + sin2@iy) |, (35) 


omitting the factor e~/6?. Amplitude A is determined by power P carried 
by the mode. From eq. (67) with du/dy given by eq. (92) and nj = A 


sy) Zee 


‘Ale 
[4] a nx aBka J?(u) 





(36) 


if P = \. 
For the HE,,, modes with y ~ —1, on the other hand, 


E; aime E ( u) (cos2¢ i, + sin2¢i,) +--- i (37) 


r 
a 
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where the dots represent terms that vanish as ka — ©. The amplitude 
|A| for P = ¥ is still given by eq. (36).* 

An important property of the field distribution (35) is that Ey — 0 as 
ka — «. Thus, in the limit’as ka — ~, the field becomes polarized in one 
direction, regardless of the value of the surface reactance X, (unless, 
of course X, = 0). From eq. (35), the amplitude of E, is porportional to 
the ratio 


y 

be (38) 
Therefore, in order that E, be negligible over the aperture of a feed, it 
is sufficient that the aperture diameter be large and the thickness ¢ of 
the disks (see Fig. 1) be small compared with their separation.t The far 
field of an aperture illuminated by the fundamental mode, the HE; 
mode given by eq. (35) for u = uo; = 2.4048, is discussed in Ref. 17. From 
a comparison of the radiation patterns of FE, and E,, we find that the 
ratio C2 between the maximum value attained by |E,|? and |E,,|? (which 
occurs on axis) is given by 


oreo (Z)’=—__*4 __, 
ka rw |]2 
| — t/h)ka tan aa 


WoO 


where w, denotes the frequency for which y = 0. One can easily verify 
using this formula that C2 remains less than 0.000316 (—35 dB) over a 
frequency range w; < w < 1.93 w;, provided ka > 10 and t/h < 0.1. 

Thus, good performance over a wide frequency range is possible, 
provided all the power incident at the input of the feed is converted to 
the HE,; mode. If, however, some of the input power is converted into 
some of the HE;,, modes, then, according to eq. (37), the field over the 
feed aperture will contain a cross-polarized component whose amplitude 
is essentially independent of the ratio y/ka. The resulting cross-polarized 
component of the far field is discussed in Ref. 17. If w1 < w < we denotes 
the frequency range over which only the fundamental mode (HE) 
propagates, it is pointed out in Ref. 17 that the largest value that wo/w, 
can assume is 1.6839; this value is attained for b/a = 1.8309. Cutoff fre- 
quency formulas are derived in Appendix D. 

In Appendix B, the properties of a surface-wave mode that can exist 
in a corrugated waveguide, in addition to the modes of eqs. (85) and (37), 
are briefly described. 


* In eqs. (35) and (37), only the leading terms for the symmetrical, asymmetrical, and 
cross-polarized components are retained. 
+ Note that from eas. (7), (8), and (18), y increases with t/h. 
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IV. COUPLING COEFFICIENT BETWEEN TWO MODES 

Suppose the electric field E, at the input of a corrugated waveguide 
is known, and we want to determine the resulting amplitude of one of 
the modes excited in the corrugated waveguide. We have to evaluate a 
surface integral of the form 


f f, (E, X Hi) -i, dxdy, (40) 


where Hg is the magnetic field of the mode whose amplitude is to be 
determined. This integral, identical to that involved in determining the 
far field radiated in a given direction by an aperture containing the field 
Kj, is in general difficult to evaluate. However, in many cases, we can 
assume that 
oe (81K, (41) 
Oz 

where 8; is a constant. This condition is approximately satisfied,* for 
instance, in the case of a feed aperture illuminated by a single mode 
propagating in the z direction with propagation constant 8 ,. We will show 
that the above surface integral can be reduced to a line integral which 
can be evaluated straightforwardly. We use the symbol (E,, He) for the 
integral (40), and call it the scalar product of the two modes Ey and 
Hp. 

If E,, H; and Eg, He are two solutions of Maxwell’s equations, in free 
space,2° 


V-(E, X H; + E; xX Hj) =0 (42) 


in the absence of sources. Now, let the z dependence of the two solutions 
be given by 


eJiz and = e~J622, (43) 
Then, in eq. (42) 

V = Vi — j(Bi — Bediz, (44) 
where V; is the transverse part of V. Therefore eq. (42) gives 


Vi (By X Hg + Ey X Hy) = J (81 — Be)[(Ei X Ha) -i, + (Eg X Hy) -i,]. 
(45) 


Next, consider a new solution E), H{ with propagation constant 6, = 
—8, and with z components given by 


Ey, = —E,z (46) 
* Of course, condition (41) is satisfied exactly by a mode in a cylindrical waveguide. 
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Hi, = +H,. (47) 


Then the x, y components of E}, H, simply coincide?! with the x, y 
components of E,, —Hj, 


E\, a Fy, Ey a Fy, (48) 
Hi. = —-A,,, Ay = sad GE (49) 


Therefore, replacing in eq. (45) 61, E,, H, with —6;, E,, H, and making 
use of eqs. (48) and (49), we obtain 


V,- (E, X Hj + E, X Hj) = —j(61 + Be)[(Ei X Hy) 
— (E3-H))] X iz. (50) 


By adding eq. (45) to eq. (50), we obtain 
(E, X H})-i, = V;- F, (51) 


where 
———\|E, x H#5+E,xH 
a erm 2 2 ul 


crn [(E, X H5) + (Ey X Hy)]. (52) 


We now integrate eq. (51) over a finite area S of the plane z = 0, making 
use of the divergence theorem, 


Sf, Ei XH) is dedy = f° F-nds, (53) 


where C is the contour of S and n is the outward normal. To determine 
F - n, let 7 be a unit vector tangent to C, 


rT=i, Xn. (54) 
Then, if A and B are two arbitrary vectors, 
(A X B) - n=A,B, — A,B,, (55) 
where A,, B, are the components of A and B in the direction of r. 
Therefore, from eq. (52), taking into account eqs. (46) to (49), 
Bo 
Bi - B 





F-n=-j 5 [EL -Ho, + E3,A12] 
2 





+ J 63 [E12H3, + E3,H,,]. (56) 


ae 
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Finally, from eqs. (53) and (56), we obtain the desired result, 


(E,,H2) = f f, (E, X Hi) «i, dxdy 


Bo 
Bi — B5 


| 





nes E,,H’, + E%.Hj,)d 
J Sg. irl 9 orf 1, )ds 


By 
Bt — B5 

Thus, the scalar product (coupling coefficient) of two modes E, and 
E, can be determined straightforwardly from the values of E, and Eo 
on the contour of the aperture S. This result has a number of applica- 
tions. It can be used, as already pointed out, to determine the far field 
radiated by an aperture S with known field distribution E), in which 
case Hp is the magnetic field* of a plane wave with propagation vector 
k and eq. (57) gives, except for a constant independent of k, the field 
component radiated in the direction of k with the polarization of H». In 
this article, we are interested in the special case where S is a circular area 
of radius a, in which case we can replace in eq. (57) 7 with ¢, since 





g. (E,,H>, + E3,Hy,)ds. (57) 


T = ig. 


If E;, H; (i = 1, 2) represents a mode of a corrugated waveguide of radius 
a,so that forr =a 





Eig = 0, ZoHig = —Jyikiz, (58) 
then eq. (57) simplifies to 
1 By * 
E,, He) = - — — E,,E5, ds. 
(Ey, He) ZB — BB? yw) f, 1zlvo, As (59) 
Since the modes are characterized for z = 0 andr = a by 
Eiz = nid 1(u;) cos¢, (60) 
where 7; is the coefficient A of the ith mode, then 
an m 
(Ky, Hy) = — ®t yyy) miniilund(u). 61) 
0 Bi — Bs 
Note that 
a?(BT — 6) = us — uj. (62) 
If we assume 
yo = yi + dy, (63) 
ug =U, + du, (64) 


" * There are two cases (two polarizations) that must be considered, for each value of 
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from eq. (61) we obtain for the power carried by the mode E, 
1 3 
= 9 (Ey, Hy) = — —- > nd f(y). (65) 


The derivative dy/du;, which appears in this expression, is calculated 
in Appendix C. In the following sections, we choose 





(E;, H;) = 1, (66) 
in which case from eq. (65) 
2Z uj du; 1 
| aN re | ae : 67 
In a°x Bia | dy | Ji (uj) oP 





These results are now applied to the problem of a junction between two 
different waveguides. 


V. JUNCTION BETWEEN TWO WAVEGUIDES OF DIFFERENT SURFACE 
REACTANCE 

Let two waveguides of different surface reactance, but the same di- 
ameter, be jointed at z = 0. Assume a single mode incident on the plane 
of the junction from the region z < 0 and let E;, H; denote the transverse 
field components. To determine the amplitudes of the reflected and 
transmitted modes, we expand E; and H; on either side of the junction 
in an infinite series of modes, and then require continuity of E; and H; 
at the junction. A simple solution for the amplitudes of the scattered 
modes is then obtained assuming the difference in surface reactance is 
small. This result will be extended in Section VI to the more general case 
of two waveguides of slightly different diameter. 

Let the transverse fields for z < 0 be represented by a superposition 
of the modes of the waveguide occupying the region z < 0, 


E; = Aje\-#42 + Ree, (2 <0) (68) 
1 


Hy = Arye“ — ¥ Rihyei##, (2 <0), (69) 
where 
Ajeye 812, Ayhye~J812  (Re(B1) > 0) 
are the transverse field components of the incident mode, and 
R;e;e/82, —R; hes | 
are those of the reflected modes. 
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Similarly, for z > 0, 
E, = 5 Tyee, (z > 0) (70) 
1 


H, = > T;hie-i8iz, (2 > 0), (71) 


where TJ; are the amplitudes of the transmitted modes. 
We assume that e;, h; are normalized so that 


(e;, hy) = f f, (ey Xhi) + i, dxdy = din. (72) 


Similarly, 
(ei, h,) = din. (73) 
Since (e;, h;) represents twice the power carried by the th mode e;, this 


power becomes imaginary if the mode is cutoff, in which case eq. (73) for 
t = 1 should be replaced with 


(e;, hj) =). 
However, in this article the calculation of R;, T; is restricted to the modes 


that are not cut off by the two waveguides. 
From eq. with n;, nz given by eq. (67), 





(ei, h,,) = (74) 


a2 ne Poe Voc eo 


and 


7 


Cae 2 (ef (75) 


where (e;, h,,) and (e,, h;) are scalar products defined as in eq. (72) 
and 


= (ka)? — (6;a)?, (76) 


a being the radius of the two waveguides. In eq. (74) 1/y and 1/y’ are the 
normalized surface reactances of the two waveguides. 
Now assume y’ — y is very small and let 


dy =y’—y. (77) 


To determine R;, T;, we require continuity of E, and H; for z = 0, 


Avie Rie, = Ties + Y Te", (78) 
1 
Ah, - ¥ Raby = Th; + Tuhy, (79) 
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Take the scalar product of the first equation with h; and of the second 
with e;. One obtains, taking into account eq. (73) and assuming that the 
mode e; is not cutoff, so that 8; is real, 


Ax(ex, h;) + Ri(e;, hj) + XS Ralen, hj) = Ti, (80) 
n#ly 

Ax(e;, hy) — Ri(e;, hy) -— SX Rnale;, hy) = T;. (81) 
n#l1,i 


Now, assume for the moment that y, y’ ~ ©. Furthermore assume none 
of the modes under consideration is at cutoff. Then 


(en, hi), (e, hn) (i #7) (82) 


are small quantities of the same order of dy. Furthermore, as we show 
below, this is true also for R,,. It follows that the two sums involving R, 
in eqs. (80) and (81) are of order higher than dy. Therefore, subtracting 
these two equations and neglecting terms of order higher than dy, 


= (ex, hy) ~ (ei, hy) 
a OR rr aa =) 


Adding eqs. (80) and (81), and neglecting terms of order higher than 
dy and solving for T;, we obtain 


(e1, h;) + (e; hy) 








T; = Ay 5 (84) 
Using eqs. (74) and (75), we rewrite eqs. (83) and (84) in the form 

jor oF / bi ae dy/du; 

R; = —Ay 85 

By + 3 Vee v [ae dy/du;|’ 7) 
a u,u; 1 

T; = Ay >. .186) 

a? By - =r ee ie | (dy/duy) (dy’/du’)| 


The derivatives dy/du, and dy’/du; are derived in Appendix C. 

It is interesting to note from eq. (85) that the reflection coefficient for 
the mode i = 1 is simply 

pi= pes Bi Bi ; (87) 
Ay By + By 

which coincides with a formula derived by Brown!8 from a principle of 
conservation of momentum. However, that derivation is not applicable 
to the present problem, which involves hybrid modes. Measurements 
of p; described in Ref. 17 show that this formula, although derived as- 
suming y’ ~ y, is quite accurate even for relatively large differences be- 
tween y and y’. 
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It is also interesting to note that the following interpretation can be 
given to eq. (86). If E, for z = 0 were known, we could determine T; 
simply using the formula 


T; = (E;, hj), (z= 0), (88) 


which follows from eq. (70), in view of the orthogonality relations (73). 
Now, if y — y’~ 0, E; does not differ much from Ae; and, therefore, we 
might be tempted to write in eq. (88) E,; ~ A ,e;, in which case we would 
get 


T; ~ Ax(e1, hj). (89) 
Alternatively, since from eq. (71) we also have 
T; = (e;, H;) forz =0, (90) 
we might be tempted to assume H; ~ A,h, for z = 0, in which case 
T; ~ Ax(e;, hy). (91) 


Neither of the two formulas is correct* even if dy ~ 0. However, according 
to eq. (84), a correct expression for small dy is obtained by taking the 
average of the two formulas. We now treat two special cases. 


5.1 Limiting case ka >> 1 | 
Assume that both y and y’ are finite, but the radius a is very large, 


ka > 1, 


a condition which is often satisfied near the aperture of a feed. From eqs. 
(31) and (33) 
du lu 


lim —=--—. 92 
pay 2 ka me) 


Furthermore, for large ka, 


lu? 
~ka--—, 93 
Ba ~ ka Sha (93) 
since (Ga)? = (ka)? — u2. Therefore, 
Lu;?—u?__1u;\2 
/— gB,~ —--—+ + ~-(—) Sy, 94 
ie a lea ae oe ” oe 
since from eq. (92) 
; lu 
‘SZ uj — =— by. 95 
be (95) 


* They are often used, however.?2 
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Using these results, from eqs. (85) and (86) we obtain 


T, = Sy ati Sy) (96) 


uz—u?_ ka 
R= A, Li yo=y 
po AY 5 
4 (ka)? ka 
One can show that these formulas are valid even if dy is not small, pro- 
vided both 


(97) 


are small. 
An application of eq. (96) is considered in Section 7.1. 


5.2 Case 1/y =0 


At the input of the feed of Fig. 1, the corrugated waveguide is con- 
nected to a smooth waveguide (1/y = 0) of the same diameter. We now 
wish to calculate the reflection and transmission coefficients of such a 
junction. Thus, assume y ~ ~ for z < 0. For y ~ ~, there are two types 
of modes: TE modes, in which case y ~ 0, and TM modes, in which case 
+ ~ ~. In the former case, from eq. (179) of Appendix C 


. dy kau = 
lim — = -—y2(u2 — 1) ———__., ~ 0). 98 
nae du yy ) (ka)? — u2 a) (98) 
In the latter case, from eq. (180) 
dy_ su ' ie 
a Bree tere ~ oo), 99 
ae (y ) (99) 


Now let the incident mode be a TE;; mode. We distinguish two cases 
depending on whether the ith mode is a TM mode or a TE mode. In the 
former case, from eqs. (85) and (98) 


B; — B; BiBy 1 
HiRes Ay ay (100) 
yore * Bit B; Vil Dk 
where U is the first root of J;(u;) = 0, 
u; = 1.8411. (101) 


If, on the other hand, the ith mode is also a TE-mode, from eggs. (85) 
and (99), 
sot f6 ibe T 
fake Bie E 
aie By +B; ° Be uy-1 





(102) 
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Fig. 2—Reflection, transmission, and coupling coefficients for input junction of Fig. 1. 


From eq. (86) we obtain, using eq. (98), 
1 : 
lim T; = ee oe ee ee (103) 
yoo aB,—aB; * kaB(uz- 1) \/ dy’ 
du, 





where dy’/du; can be determined using eq. (178), unless y’ > 1, in which 
case we can use eq. (98) or (99) with y,u replaced by y’,u;. 

Equations (100) to (103) have been used to calculate the behavior of 
a junction with b = 1.8309a. Consideration has been restricted to the 
TE; mode and the TM); mode of the smooth waveguide, and the cor- 
responding modes (HE; and HE;,) of the corrugated waveguide. The 
results are shown in Fig. 2, wherei = 2 refers to the TM); mode (or the 
HE}, mode), 








R,|2 
2 1 
=i 104 
ars (104) 
is the input reflection, 
2 — |R2|? (105) 
2 VA 
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gives the power converted into the TM; mode, and 
T 5 
Ay 
gives the power converted into the HE,, mode. In Fig. 2, we is the fre- 


quency at which y’ = ~. The corrugated guide at this frequency behaves 
like a smooth guide and, hence, 


p1 = p2 = t2=0. (107) 


The curves of Fig. 2 are useful in determining the practical bandwidth 
of the junction of Fig. 1. 


9 2 
to= 


(106) 








VI. JUNCTION BETWEEN TWO WAVEGUIDES OF DIFFERENT DIAMETER 


For some applications, to minimize the input reflection of a corrugated 
feed, it may be convenient to choose for the smooth waveguide a diameter 
different from that of the corrugated waveguide. In this section, the 
analysis of Section V is extended to the general case of a junction between 
two corrugated waveguides of different diameter. Let a and a’ be the two 
diameters for z > 0 and z < 0, respectively, and assume again a single 
mode is incident on the junction, from the region z < 0. 

If E, for z = 0 were known, then the transmission coefficients T; which 
appear in eqs. (70) and (71) could be determined at once using the for- 
mula* 


T. = ey (E, X h?)-i, dS, forz =0, (108) 


which follows directly from eq. (70) in view of the orthogonality of the 
modes e;, h; [see eq. (73)]. In eq. (108) S’ denotes the circular area 
O0<r<a’. 


Now, for z = 0, E; is given by eq. (68) inside the area 


0<r<a, (109) 
and it vanishes for a <r <a’. Therefore, eqs. (108) and (68) give 
T; = Ay(ey, hi) + > Ra(en, hd), (110) 
n=1 
where 
(ep, hi) = f i. (e, X h7)-i, dS, (111) 


S being the circular area (109), which corresponds to the waveguide of 
the region z < 0. 


* Here we are only interested in calculating T; and R; for 8; > 0, B; > 0. 
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Equation (79), which was obtained by requiring continuity of H; for 
z = 0, must be satisfied over the area S. By multiplying this equation 
with e,, and integrating over S, we obtain forn + 1 


ARES Tilex ho (ae): (112) 
1 


If the coefficients T; in this relation are expressed in terms of the coef- 
ficients R; using eq. (110), we get forn ¥ 1 


-R, = (Ay + Ry) © (er, hi(en, H+ Rn E (en, hi)? 
1=1 i=l : 


+ ORY (ey Wilen, hi). (113) 


sxn,l1 t=1 


For n = 1, the second sum of the right-hand side should be omitted and, 
furthermore, —R,, should be replaced with A; — Ry. 

We have thus obtained a system of equations in the unknowns Rj, Ro, 
etc. We solve* them in the limiting case where both a’ — a and y’ — y are 
very small, in which case 


(e,,h,)~0O forn ¥s 
R;~0 


; ; 114 
(e,, h’,) - 10 ae) 


and therefore the first two terms of the right-hand side of eq. (113) for 
n ~ 1are respectively equal to 


Ai[(en, hy) + (e1, hy)] 


and R,,. The last term can be neglected. Therefore, eq. (113) gives for 
noel 


1 t t 
R,~ — rs [(e1,h,) + (en, hy)JAr (nm #1). (115) 
Similarly, for n = 1, 
1 ; ; 
Ry 5 [1 — (e1, hy)?]Ai * [1 — (ex, hy)JAj. (116) 


The transmission coefficients can now be determined using eq. (110). 
We find forn # 1 


T, ~ lle, bh Seg hla, md), (117) 


which is a generalization of eq. (84). 


* This derivation is not rigorous, for we neglect to examine the question of convergence 
of the summations in eq. (113). However, the validity of the results appears to be confirmed 
by the experimental results. 
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The coefficients (e;, h,,) can be calculated using eq. (57). If y’ — y and 
a’ — a are very small, we can proceed as follows. The field components 
e, and h,, of the nth mode are considered to be functions of the coordi- 
nates r,@ and of the two waveguide parameters a,y. Therefore, 


oh, oh, 
oy oa” 


where h,,, 0h,/Oy and oh,,/da are evaluated for a’ = a, y’ = y, and oy 
are da denote y’ — y anda’ — a. asimilar relation can be written for e,,. 
It follows from eq. (118) that (e;, h;,) for i ~ n is asum of two terms, 


dh dh 
hi) = (er *) oy + j *) 6a, 11 
(e, hi =(e ay) oY (e 5a) a (119) 


since (e;, h,) = 0. The first term is simply the coefficient (e;, h,,) cal- 
culated for a’ = a; it corresponds to a junction between two waveguides 
of the same radius, but different surface reactance. The second term can 
be interpreted as the coefficient (e;, h,,) relative to a junction between 
two waveguides having the same surface reactance but different radii 
a and a’. Since the term has already been treated in Section V, only the 
latter need be considered. If one sets 


oh, On 
6h, = a 6a, de, = Sy 


and if the ¢ variations of both modes are of the type considered in Section 
I, then, taking into account that e;, = 0 for r = a, using eq. (57) we 
get 


(e;, 6h,,) = 7aoda | -i Be 2 (c 


h, = hy(r, ¢, a’, y’) ~ hy, 





(118) 














da, (120) 





Ohne 4 On nd py, :) 








BPS GN tod ea eee 
. L hi nz 
FP J 3? — Bi B 5 (ee o) ng + °e his) »-o- | (121) 
i n a r=a 


As an application, consider y = ~, in which case (e;, h,,) can be in- 
terpreted as the coefficient (e;, h,,) relative to a junction between two 
smooth waveguides of radii a and a’ = a + 6a, respectively. Assume e; 
isa TE mode and e,, is a TM mode, so that for r = a 
nz 

oa 


where the last term vanishes because h, is a TM mode. Then eq. (121) 
gives 


Ol ng 
(e;, dh,) = rada —— B — B zal —JBn ea en his) o=o0 





Cig = Ciz = =0, (122) 


ot POC Hg 
+ 10: (Shia) gece | (123) 
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Similarly, interchanging n ss 1 in eq. (121) and taking into account that 
forr=a 


0€ jz 





= €ng = Cnz = Anz = 0, 124 
= Cnd = & z ( ) 
we obtain 
(e,, dh;) = 0. (125) 
Now, the two modes are characterized by 
1 
hiz= znd é ui) sing, (126) 
Zo a 
‘ 8 J c ui) 
J ja a 
hig = — 2-1; ——— co, 127 
$ Zo tas : cos (127) 
aye 
and 
r 
Cnz = nd 1 (- Un) cos¢, (128) 
r 
Fi (jun) 
Bra 1 a n 





Ene = Jn sing, (129) 
Un 


—Un 
a 


where the amplitudes 7; and n, are, because of the requirement (72), 
given by 


ee 2 re re ee on ae (130) 
m na av B; Vu? —1|Jy(uj)| * (ka) - 
_4/2Z01 Un I ame, (131) 





me Vira aV Bp |Ji(un)| ¥ (ka) 
From eggs. (123) and (126) to (131), taking into account that J;(u,,) = 0, 
we obtain the final result 
ka ik 6a 
VaBiaBy, Vur —1 a 
Note that in deriving this relation it has been assumed that 6a is suffi- 
ciently small so that 


(e;, dh,) = 2 (132) 


=2 
vB 


~~ 


= 
= 
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If this condition is not satisfied, we should replace 8, with £,, in eq. 
(132). 

Of special interest is the case where e; and e, represent the TE;; mode 
and the TM, mode, respectively. In this case u; = 1.8411 and, letting 
t = 1andn = 2 for these two modes, we get 


ka 6a 
VaBjaBo a- 


From eqs. (115), (117), (125), and (133) we then obtain for the conversion 
coefficients T> and Ro 


(e1, dhe) = 1.2937 (133) 


fia heats = SA, (134) 
apjap2 a 
where |A,|?, |T2|?, and |Ro|? represent the incident power, and the 
powers transmitted and reflected in the TM;; mode. We can verify* that 
Tis smaller by a factor of 2 than the conversion coefficient given in Ref. 
22, which is due to the fact that the assumptions of Ref. 22 imply 


ia ons (e1, h,), (135) 


rather than eq. (117). 
Note that for B2 > 0, we have af) — u2V éa/a , and therefore 


“Us 3/4 
T2* 0.646 \/ 2 (“*) (136) 


apy a 


Note T,2 remains finite even when the TM; mode approaches cutoff 
in the first guide. 


Vil. MODE CONVERSION IN A NONUNIFORM WAVEGUIDE 


Typically, a corrugated feed is made of one or more sections of non- 
uniform waveguide whose surface reactance and radius are functions 
of z. Since a nonuniform waveguide does not in general possess a natural 
mode of propagation, an incident mode will be scattered in forward and 
backward modes. This is true even for a conical waveguide of constant 
surface reactance (except when y = 0 or y = ~). The analysis of Sections 
V and VI gives the differential scattering parameters which allow the 
local coupling into forward and backward modes to be determined at any 
point in a uniform waveguide. We can thus obtain a set of differential 
equations, whose coefficients are given by the above scattering param- 
eters, and which can be solved, at least in principle, for the mode am- 


* In Refs. 22 and 23, the TM); mode was cut off to the left of the junction, and for this 
reason there is poor agreement between those measurements and eq. (134), which is not 
applicable in this case. However, numerical calculations by Masterman and Clarricoats 
agre well with eq. (134) at frequencies well above the cutoff frequency of the TMi; mode, 
as we may verify from Fig. 11 of Ref. 24. 


CORRUGATED FEED PERFORMANCE _ 857 


plitudes. We confine ourselves to a first-order treatment assuming the 
total scattered power is much less than the incident power, since this is 
the most interesting case if the feed is well designed. It is convenient to 
assume for the moment that only y varies with z, in which case the 
waveguide can be approximated by a succession of junctions of the type 
considered in Section V. Let the HE; mode be incident at the input (z 
= z,). We wish to determine the resulting amplitude To(z) of the HE}, 
mode for z = 2. If the variation of y is sufficiently slow, we can neglect 
reflections and determine T'2 assuming the amplitude A, of the HE, 
mode is nearly constant. The transverse field of the fundamental mode 
is then 


Ajeje J P12), (137) 
where 
d® 
or = Bi. (138) 


The effect on the HE,, mode of a small variation dy at z = £ is to produce 
at z ~ 22 a component 


dT» = to(t)6yA ye P18) e —/1Palz2)— P2(8)] (139) 
where 
d® 
7s == Bo, (140) 
z 


and from Eq. (86), 


—1 uUyu2 1 
= en 141 
eat) apy — abo Von | (dy/du;)(dy/duyz)| oe 


Note that both 6; and 82 are functions of z. From eq. (139), integrating 
from 2, to Zo, 





T2(Z2) = i- dT, = Aye —J P22) 
21 


x f- to(z)etslo22)—o1@)] dy dz, (142) 
21 dz 


which assumes a very simple form when ka >> 1, as discussed in the fol- 
lowing section. 


7.1 Conical horn with ka > 1 


Suppose the radius a varies linearly with z, as shown in Fig. 1 for z > 
z,, that the flare angle a is very small, and 
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ka(z) > 1 


for 21 < z < zg. It was shown in Section III that for ka > 1 the properties 
of a mode are entirely determined* by y/ka and, therefore, a mode will 
propagate without variation of its transverse field distribution if y/ka 
= constant. For this reason, it is reasonable to assume that mode con- 
version will be negligible if 


5 cE) = 0), (143) 


Under this assumption, the effect on the amplitude of the HE}, mode 
for z = z2 of asmall variation 6(y/ka) occurring at z = £ can be expressed 
in the form 


dT. = r9(é)6 G) Aye -i#1®¢ -ild2(e2)—#28)], (144) 
a 


where, since eqs. (144) and (139) must agree in the particular case a = 
constant, 





72(f) = kato(é). (145) 
From eq. (96), 
n= —-— 5 (ka > 1). (146) 
ug — U; 


Therefore, 


To(z2) = —Aye~sP2(22) A142 f * pilbie)— P20] a) (147) 
a 


Ug — UY YZ) 


Note that, since ka > 1, 


lu? 
i;Ska--—. 148 
ae A 2ka ee 
Therefore, from eqs. (138) and (140) 
d kug—uj 
— [@ —@® Se a 149 
Now, a varies linearly with z, 
a = (z — Zo) tana, (150) 


* Since now we are dealing with a conical waveguide, each mode is a spherical wave 
centered at the apex A of the waveguide, and the field distribution over a spherical 
wavefront is given to a first approximation (small a) by eqs. (35) and (37). To obtain the 
field distribution over a plane z = constant, we must therefore introdce in eqs. (35) and 
(37) a factor of exp (—jy), ¥ = k(x? + y2)/2R, R being the distance of the plane from the 
apex. 
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where Z9 is the value of z at the apex of the waveguide. Therefore, from 
eqs. (149) and (150) 


P1(z) — So(z) = — : 
1(2) 22) 2 ytana ka 
Of particular interest is the case 
y = constant. (152) 


Then eq. (147) can be readily integrated. Taking into account eq. (151), 
we obtain 


4u?usy? tana? 
2 2 
(uz — u;)4 


2 9 
1-exp [77S = (1-4) || : (153) 


2 y tana kaj, ao 


[Pel? = |Aal? 


x 





where a; = a(z;). Therefore, 

16u?u3y2 tana? 
(ug—uj)4 ” 

where it is recalled that wy = 5.1356 and uy = 2.4048. For a = 4°, which 


is the value chosen in the experiment of Ref. 17, this inequality gives for 
y=l 


|Ts|2 <= |Ay|? (154) 


| T»| 2 -4 
> = 0.6636 10-4, (—41.8 dB), (155) 
|Ai| 
which is a very small value for most aplications. For a = 16°, on the other 
hand, we obtain —29.8 dB, which may no be negligible. 
Note that | 7's|? and |A,|? are, respectively, the powers carried by the 
HE}; mode and the HE;; mode. 


Vill. SUMMARY 


In the feed of Fig. 1, when a TE}; mode is incident at the input, some 
of the incident power is in general reflected. Furthermore, some power 
may be converted to unwanted modes if the corrugated waveguide 
supports more than one mode at the input. Additional mode conversion 
may take place inside the feed if the variation of the radius and of the 
surface reactance is not gradual enough. As a result, a feed will have a 
nonzero input reflection and, at some frequencies, unwanted modes may 
illuminate the aperture of the feed. The consequences of these unwanted 
modes on the radiation characteristics—e.g., enhanced cross polariza- 
tion—are pointed out in Section III and in Ref. 17, where the theory is 
compared with experiment. 

These effects can be evaluated to good accuracy using the expressions 
derived in this article, as highlighted below. For large ka, eq. (35) ex- 


860 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1977 


presses the field shape for all copolarized hybrid modes of various radial 
harmonics with the same ¢ dependence, while eq. (87) corresponds to 
the cross-polarized modes. Equation (36) gives the mode amplitudes 
required to normalize the power carried by the modes. 

A property of corrugated feeds is that the aperture field distribution 
does not remain constant with frequency, as in the case of a feed with 
smooth walls, but varies because of the frequency dependence of the 
surface reactance X,. Thus, although the desired mode has no cross- 
polarized component at the resonant frequency of the corrugations, at 
other frequencies the desired mode does radiate some cross polarization. 
The ratio, C2, between the maximum value of the cross-polarized power 
in the radiation pattern and the maximum value of the copolarized power 
(which occurs on axis) is given by eq. (39). From eq. (39), it follows that 
with large ka and thin disks one can maintain low cross-polarized power 
in the radiation pattern from the desired mode over an octave or 
more. 

At a junction between waveguides of the same diameter but of dif- 
ferent surface reactance, eq. (85) gives the general expression for the 
mode coupling coefficient to modes reflected from the transition, and 
Eq. (86) gives that for modes transmitted forward from the transition. 
Equations (97) and (96) are simplifications of eqs. (85) and (86), re- 
spectively, which apply for ka > 1. When the input waveguide is smooth, 
the mode coupling coefficient is given by eq. (100) for reflected TM 
modes, by eq. (102) for reflected TE modes, and by eq. (103) for hybrid 
modes transmitted forward from the transition. Since the transition from 
smooth to corrugated waveguide is a major source of unwanted modes 
in a corrugated horn, eq. (103) is very useful in determining mode purity. 
Another important formula is eq. (87), which determines the reflection 
coefficient of the dominant mode (return loss) for any transition in X, 
and any ka. 

Another source of generation of undesired modes is the mode con- 
version occurring along the conical taper of a corrugated horn. Equation 
(154) gives the mode-coupling coefficient for the transmitted undesired 
mode due to a conical taper. 

In some cases a step in diameter may be used to match transitions 
between different surface impedances; eq. (134) determines the mode- 
coupling coefficients at a step in diameter. 


APPENDIX A 
Asymptotic Series for u and y in Terms of 1/ka 
We determine the asymptotic series for u and y in terms of 
i) 
os 156 
ha (156) 
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under the assumption y ~ 0. It is convenient to introduce in eq. (20) the 
quantity!? 


_udi(u) _— Jou) _ 





=u (157) 
Jy(u) Jy(u) 
Note that 
dF Jo(u) mos Ji(u) Jo(u) 
—=—— -u- (158) 
du J(u) “ Ti(u) Jy(u) 
But from eq. (157) 
Jo(u) _F if 1. 
Jiu) uw 
Therefore, eq. (158) gives 
a = - 2 FP + (2-1). (159) 
du u u 
It follows that 
Pe eae en (160) 
du 
2 
w+ (14 OF) + 2u = 0, (161) 
a+ ee Satr+2(S) +2=0, (162) 
du? du 
etc. 
In terms of F’, eq. (20) can be rewritten 
4 1\2 
yu2— F— F2+1—-wu? a) = 0, (163) 
ka ka 
or, using eq. (159), 
1 dF TAZ 
ae du ui Gall : ee 
Now assume 
1 1\2 
w= wom | 1+ a 2 + a (-- ) tee]. (165) 
ka ka 
Develop F in a Taylor series about the point u = Ugm, 
dF 
F = F(uom) + Gs) ne ee Ore (166) 
du u=uom 


From egs. (160) to (162), with u = ugm, taking into account that 
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F(uom) = —1, (167) 


we obtain for the derivatives appearing in eq. (166), 


=) 225 

Gs u=U0m is 

d2F 

— = = 3, 168 
ee ( ) 


(=) as He 
du3 u=Udm Uom . 


Substituting eqs. (165), (167), and (168) in eq. (166) one obtains F as a 
series of powers of 1/ka; the coefficients of this series are algebraic ex- 
pressions in Ugm and ay, ae, etc. Similarly, by developing dF/du in a 
Taylor series about the point u = Uom, and then using eqs. (165), (167), 
and (168), we obtain dF/du as a series of powers of 1/ka. Substituting 
eq. (165) and the above series expansions of F and dF/du in eq. (164), 
we can solve for the coefficients a1, ag, etc. We obtain eq. (31). Substi- 
tuting eq. (31) in eq. (166) we obtain an expansion of F' in powers of 1/ka. 
Using these results, from 


etc. 


F 
- =— ee eee (169) 


cos6, u\2 
1 = pate 
Vv i 
we get eq. (32). 


Equations (31) and (32) have been obtained assuming eq. (165), which 
corresponds to the limiting case of eq. (28). If, instead of assuming eq. 
(165), we assume 





u= Un = an | 1+ 6 (;-) +62 (,-) ++ | (170) 


we obtain eqs. (33) and (34). 


APPENDIX B 
Surface Wave Mode 
In addition to the modes considered in Section III, there is a mode for 
which 8 > k. Thus, since for this mode cos 6; > 1, it is convenient to re- 
place 6; with j6, in eqs. (16) and (18). Since 
cos j6 = cosh 6; 

sin JO; = J sinh 0; 
J1Gx) = j11(x), 
Ji Gx) = Iy(x), 
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(171) 


where I(x) is the modified Bessel function of order 1, we obtain from 
eqs. (16) and (17) 





I\(u) ou 
eg 172 
I,(u) cosh 6; i) 
Ses are ar (173) 
ka u Iy(u) ueyu 


where u = ka sinh 6,. We can verify from these two equations that u > 
«© as ka — ~, Now, for large u, 


T,(u) =I,(u) oe (174) 


Therefore eq. (172) gives 
y~o, aska>, 
From eq. (173), for large u and ka, 


y 1 : 
—N--, 175 
ka u a 


Therefore, since u = ka sinh 6, 
1 
sinh 0, ~~ — _ : , (176) 


We now examine the behavior of the field components for ka > . 
Taking into account eq. (171), from eqs. (10) to (15), after replacing 6; 
with j0,, we find for y = —~@ (i.e., for B = 0) that the only nonzero com- 
ponent of the magnetic field is H, and 


1 ; r . 
—~Z Hs = ——— Al ( -) — jBz, 
re sinh 6, a a ous 


Therefore, for kr > y 


-ZoH, ~ A’ Vs cos¢ exp[k(r — a) sinh 6; — j6z], 


where A’ is a constant. This shows that the field is confined to the near 
vicinity of the wall, decaying approximately exponentially from the 
wall. 


We can show that ka > 1.81, the surface wave mode in combination 
with the HE;,, and HE;,, modes comprise the complete set of propa- 
gating modes whose E, azimuthal dependence is cos¢. 


APPENDIX C 
Derivation of dy/du 
From eq. (163), 
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where 
Furthermore, from eq. (164) 


Therefore, 

1 dy 2 1 ud 1 \2 1 
——~ =-—(F-~=)4+-[(—) -1|]/-—F 
ka du a 7 u Ga | u3 


“G-") Ge) ial 07 


For the case y > ~, we now determine 


lim ca 


For y ~ ~ there are two types of modes: TE modes, in which case y ~ 0, 
and T'M modes, in which case y ~ ~. In the former case from eq. (169) 
we have F ~ 0 and, therefore, eq. (178) gives 


ee 1 ka)? —u? 1 
Liar ie errr 
= yur) GG (V0) (179) 
since from eq. (163) for F ~ 0 
ingest 
eee y kau? 


In the latter case (y ~ ©) from eq. (169), we have F ~ ~. More precisely, 
from eq. (163) 


u2 
Fr—y, 
ka = 
and therefore from eq. (178) 
dy_ka,,_ ou se 
— SW — F2N —— y2 ~ o), 180 
ae wy? (ye) (180) 
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APPENDIX D 
Cutoff Frequencies of the Modes of Equations (35) and (37) 
From eqs. (16) and (17) we get 
72+ yo—-1=0, 
where 
yu? 


Qo= =. « 
ka cos6; 


Therefore, either 


Bee eee ge) 
pee aN ee. (181) 
ey) 
or 
—yW — WV a2 
yoeoxet (182) 


respectively, in the two cases of eqs. (35) and (37) (which correspond, 
respectively, to y — 1 and y ~ —1 for ka — ~). At cutoff 8 — 0;1.e., cos#, 
— 0. For cos#; ~ 0, y # 0, we have |w| — © and, therefore, from eq. 
(181) 


_,f9 ify>0 
tad Barer ve 
(184) 
whereas from eq. (182) 
o, ify>0d 
ea ify <0° 


If y = 0, the mode is of the TE type. Now, for a TE mode at cutoff, the 
only nonzero component of the magnetic field is H, and therefore the 
surface reactance X, has no effect on the cutoff frequency. This means 
that the cutoff frequency can be determined by replacing the corrugated 
wall with a smooth wall of radius a, and therefore the cutoff frequency 
is determined by the condition: 


J\(ka) = 0. 


If y = ©, on the other hand, the mode is of the TM type and the only 
nonzero component of the electric field at cutoff is E,. It follows that if 
the disks are very thin (t ~ 0), they can be removed without affecting 
the field. Thus, the cutoff frequency can in this case be determined by 
replacing the corrugated waveguide with a smooth waveguide of radius 
b. It is thus determined by the condition 
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J\(kb) = 0. 


When y = 0, y = £1, and from eq. (16), J;(ka) = 0 for cos0, = 0. Thus, 
when y = 0 both types of mode have the same cutoff frequencies. 
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A corrugated feed with nearly ideal radiation characteristics from 
17 GHz to 29 GHz has been built using a novel fabrication technique. 
The bandwidth of single-mode operation was maximized by properly 
choosing the input parameters of the feed. As a result, only the fun- 
damental mode can propagate at the input from 19 GHz to 28.8 GHz. 
In this frequency range, the far field is essentially polarized in one 
direction. At frequencies higher than 28.8 GHz, there 1s a cross-polar- 
ized component caused by an unwanted mode. An approximate cal- 
culation of the power in this mode ts given. A simple formula for the 
input reflection coefficient is provided. Results are included that show 
how to compute mode conversion in a conical taper, cross polarization 
from a corrugated horn, including contributions from spurious modes, 
and the reflection coefficient from the smooth-guide to corrugated- 
guide transition. Comparison of theory and experiment shows good 
agreement. 


|. INTRODUCTION 


It has been shown by Thomas! that under certain conditions the field 
over the aperture of a corrugated feed is virtually constant over a very 
wide frequency range. For this behaviour to occur, the radius a of the 
aperture must be much larger than a wavelength, i.e., 


ka >1( =<) (1) 


and furthermore, the aperture must be illuminated by a single mode. If 
both conditions are satisfied, we can show that the field distribution is 
essentially independent of the surface reactance of the corrugated horn 
wall, X,, and, as a consequence, it is little affected by the variation of 
X, with frequency. This result is very important, for it implies that it 
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Fig. 1—Dimensions of the corrugated horn. 


is possible to design a feed so that its radiation pattern is circularly 
symmetric and polarized in one direction over a very wide frequency 
range! (an octave or more). The most difficult condition to satisfy to 
obtain such large bandwidth is the requirement that a single mode be 
excited in the feed. This requirement is discussed in a separate article.‘ 
Here, after summarizing certain results of that article, we describe the 
results of an experiment. A long feed, with a flare angle of only 4°, was 
fabricated using a special technique described in Appendix D. At the 
input of the feed, which is shown in Fig. 1, the waveguide dimensions (the 
corrugated depth / and the radius a; see Fig. 1) were optimized so as to 
maximize the bandwidth of single mode propagation. As a consequence, 
unwanted modes were cut off (at the input) over the frequency range 


Wy] <w< 1.68391, (2) 


where w; = 19 GHz and 1.6839w, = 28.8 GHz. The input reflection and 
the far field were then measured from 17 GHz to 35 GHz. The input re- 
flection agrees very well with a simple formula given in eq. (42) and de- 
rived in Ref. 4. Over the frequency range (2), the far field is essentially 
polarized in one direction. At frequencies higher than 28.8 GHz, however, 
a strong cross-polarized component is caused by an unwanted HE}, 
mode, which is excited primarily at the input of the feed. A simple ex- 
pression for the power converted into this mode is given in Ref. 4. 
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ll. PRELIMINARY CONSIDERATIONS 


In the experiment described in Section V, the feed is excited at the 
input by a linearly polarized TE;; mode. Therefore, consideration will 
be restricted to the modes arising for this particular excitation.2-8 

Consider a disk-loaded waveguide centered around the z axis as in Fig. 
1, and let r, ¢, z be cylindrical coordinates defined by x = rcos¢ and y 
= r sing. Assume for the moment that the waveguide parameters h, a, 
b, and ¢ are independent of z. The separation h of the disks, which oc- 
cupy the region a <r < b, is assumed to be much smaller than a wave- 
length A» 


kh «1. (3) 


The region between two consecutive disks forms a radial line whose input 
reactance j/X at r =a is a function of the radial length / = b — a; for ka 
> 1, 

JX ~ jZo tankl, (4) 


where Zo = V 0/e9 . Because of condition (3), the effect of the disks can 
be adequately accounted for by introducing an effective surface reac- 
tance? 


; t 
JX = JX (1 = zs (5) 


where ¢ is the thickness of the disks, and by requiring that the field for 
r <a satisfy the boundary conditions 


forr =a, (6) 





JXs 
where Ey, Hy, HE, are the ¢ and z components of the electric and magnetic 
field. 


2.1 The HE), and HE’ ,,, modes 


The properties of the hybrid modes?-° of a corrugated waveguide are 
determined by the radius a and the surface reactance X, of the wave- 
guide. In general, there is no simple relation? between the propagation 
constant 8 of a mode and the two parameters a and X;. If, however, 
condition (1) is satisfied and furthermore 


y 
ae ee, Ly df 
ha (7) 
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where 


ea (8) 


then we can show?‘ that for all the modes except one, described in Ap- 
pendix B of Ref. 4, the propagation constant is independent of X, and 
is simply given by 

Ba ~ V (ka)? — u?, (9) 


where u is a constant that is either the mth root of Jo(u) = 0, or of Jo(u) 
= 0. Thus, either 


U = Um, Jo(Um) = 0 (10) 
or 
U = Um F2(Um) = 0. (11) 


The modes to which eq. (10) applies will be called HE}, modes and those 
corresponding to eq. (11), HE, modes. The HE}, modes approach 
asymptotically, for large ka, the field distribution?4 


E=Jo é u) e— iB iy, (12) 
a 
whereas the HE}, modes approach the distribution 
E=Jd> (- u) (cos2¢i, + sin2¢i, )e —/8?. (13) 
a 


The absence of a z component in eqs. (12) and (13) is due to the fact that 
this component vanishes like 1/ka, as ka > ~. 

The HEjm modes, which are given by eq. (12), have the important 
property that the electric field is polarized in one direction. Of special 
interest is the fundamental mode (m = 1) characterized by 


u = Uy = 2.4048. (14) 


Note that both eqs. (12) and (13) are frequency independent. 

Thus, over the frequency range in which both conditions (1) and (7) 
are satisfied, the field distribution of an aperture illuminated by the HE, 
mode is essentially frequency independent.!? If, however, only condition 
(1) and not condition (7) is satisfied, then we must add‘ to the right-hand 
side of eq. (12) a component of the type (13), so that 

E=Jy(—u) i, -7 * 

a y+1 

where the factor e—/82 has been omitted. Both y and u are functions of 
ka and X;, and, therefore, they vary with frequency. If one develops jy, 


r ‘ : ; 
Jo G u) (cos2¢i, + sin2¢i,), (15) 
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u in power series of y and 1/ka one finds‘ 


w= tim [t~Sy (--)-| (16) 
y=l1-u?, [eo — Suk + us) (2) --+} (17) 


2.2 Mode conversion in a conical horn4 


Consider a conical horn of constant surface reactance and flare angle 
a, as shown in Fig. 1, from z = 2; toz = 29. If ais sufficiently small, and 
if the input of the horn is excited in the HE; mode, then the resulting 
field inside the horn is very nearly a spherical wave!’ originating from 
the apex of the horn and given by 


Ee ~tkk (18) 


where EF is given by eq. (15) and R is the distance from the apex of the 
horn. Since a is small, 

2 

R~z+—, (19) 

22 
r being the radial distance from the axis and z the axial distance mea- 
sured from the horn apex. From eq. 17, the parameter y which appears 
in eq. (15) is a function of both y and ka; assume |y| # ~, 0. Then, since 
ka increases with z, we have from eq. (17) that y varies with z, and 
therefore the field distribution (15) does not remain constant with z. This 
variation is accompanied by generation of unwanted modes, an effect 
that will be negligible only if a is sufficiently small. 

To determine how small a should be, consider the special case ka(z 1) 
> 1 treated in Ref. 4. Let P, be the total power converted from the de- 
sired mode into the HE,, mode and let Py be the power incident at the 
input. Then, we find that for z = ze 





P. = Po X 3.393 (107) y? tan2a]1 — e/¥| 2, (20) 
where 
ree ir er, 2) 
Therefore, 
P. s Po X 1.857(10~2)y? tan2a, (22) 


where the equality sign is attained for y = (2n + 1)z. In the experiment, 
a = 4°, in which case for y = 1 we find P, < 6.636 X 10-° (—41.8 dB), 
which is negligible for most practical purposes. 
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Ii. RADIATION CHARACTERISTICS2®7 


Let the aperture of the horn be of sufficiently large diameter that the 
far field is simply proportional to the Fourier transform F of the aperture 
field. Also let a be so small that we can neglect the phase variation caused 
over the aperture plane by the variation of R in (18). Then, by taking the 
Fourier transform of eq. (15), we obtain for the field distribution at a 
great distance D from the aperture, for example, using the contour in- 
tegral method of Ref. (4), 





F = No(wv)ix + 1 Nau.) (cos2oi, + sin2¢i,), (23) 
y 
where 
v = ka(zg) sind, (24) 
2 
Ng (uv) = [vIn wi(v) — uw, (vSu(w)], (25) 


and ¢, 6 are spherical coordinates defined by x = D siné cos¢, y = D sin 
0 sind. Equation (23) gives, except for a factor independent of ¢, 0, the 
far-field pattern. For small y/ka we have 


u = uy = 2.4048. 


In this case, from eqs. (23) to (25) the cross-polarization ratio C between 
the maximum value of |£,|, which occurs for v = 3.67, and the maximum 
value of |E,.|, which occurs for v = 0, is given by 


c22 Exim _ @ 962 Ga (26) 
Ops ae +1 
From eq. 17, the asymptotic value of C? for large ka is 
y 2 
C2 0.14 Ga (27) 
ka 


If ka > 1, but y/ka is not small, C2 has the behaviour of Fig. 2. 

Equation (23) assumes that the aperture of the horn is illuminated 
by the HE,; mode. If, in addition to this mode, there is also some HE}, 
mode of power P,, then we have in addition to the component in (23), 
a component? 


vx ev x 4/ Palit) ae N2(u,v)(cos¢ix + singi,), (28) 


where P,, Po are the powers carried by the two modes and y is their dif- 
ference in phase for z = Zo. For ka > 1, 


u ~ uy, = 5.1356. (29) 
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Fig. 2—Cross-polarization peak versus normalized surface susceptance. 


Therefore, from eq. (25) with J(u) = 0, 
a2u 


2 
= Jo(v)si(@)) (30) 


N3 (u,v) = & = 


whose maximum value occurs for v ~ 4.356 and is 


1.271 X 1073 X a4 J?(a) U2. (31) 


Because of the component (28), we have for y = 1 that the ratio C? be- 
tween the maximum value of |E£,|?, which now occurs for v = 4.4 and 
the maximum value of |E,| is 


P 
C2=0.194—. 32 
9 Py (32) 


Note from eq. (23) that the normalized radiation pattern of the HE; 
mode for ka > 1 is simply given by 


Jo(v) 
2 y2 


P(0) = peal (33) 


with u = u, = 2.4048 and v = ka sin 0. If 6; and 62 denote the values of 
6 for which P(6) = 0.5 and P(@) = 0.1, then from eq. (33) 
ka sin 0; = 2.078 (3-dB point) (34) 


and 
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ka sin 0) = 3.597 (10-dB point), (35) 


respectively. 


IV. FEED DESIGN?*9 


The feed of Fig. 1 is now described. Its measured characteristics are 
discussed in Section V. Let {we, w;} be the frequency range over which 
only the dominant mode, the HE; mode, propagates at the input (z ~ 
0). Let wo be the frequency at which the surface reactance is infinite at 
the aperture; then 


be) 20s) at, tian: (36) 


At the input, the ratio between the radii b and a is optimized for maxi- 
mimuM w/w ;. This is shown in Appendix A to require 





BtO), BOTSS SG a5co. (37) 

a(0) 3.8317 7 
Then* 

= a) = 1.6839, (38) 

W] @1/ MAX 


and the surface reactance can be shown to vanish for w = wo, 
y = ~ forz = 0, atw = we. (39) 


In the experiment, the frequency w was chosen equal to 19 GHz, which 
gives wo = 28.8 GHz, a(0) = 0.556 cm, and b(0) = 1.161 cm, as we obtain 
from eqs. (37) and (38) and the condition 


kb = 7.0155 at w = wo, (40) 


shown by Fig. 9 of Appendix A. 

From Fig. 1 the feed consists of two parts joined at z = z;. From z = 
0 to z = 2,, the outer radius is kept constant, so the cutoff frequency of 
the HE}, mode remains constant, as shown in Appendix A, where the 
relations between a, b and the cutoff frequencies of the various modes 
of eqs. (12) and (13) are derived. Note that from z = 0 to z = z, the radius 
a increases, which implies that at any given frequency the surface re- 
actance decreases with z. The importance of this requirement was first 
realized by Bryant®. From z = z, to the aperture, the surface reactance 
remains constant with z. The frequency wo at which it is infinite was 
chosen in the experiment 


* To put this ratio in perspective, the highest frequency of the 6-GHz (TH) common 
carrier band is 1.732 times the lowest frequency of the 4-GHz (TD-2) band; similarly, for 
the 18- and 30-GHz bands, the ratio is 1.695. 
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wo = 1.21 Wy). (41) 


Over the frequency range {w , wo} only the HE; mode is excited at the 
input, other modes being cut off. Therefore, since the variations of a(z) 
and b(z) from the input to the aperture in Fig. 1 are sufficiently gradual 
to ensure negligible generation of unwanted modes, we have that the 
aperture is essentially illuminated by a single mode, the HE); mode, over 
the above frequency range. The radiation characteristics for w; < w < 
wo are, therefore, given accurately by eq. (23) with a = a(z2) and u given 
by eq. (16) for m = 1. 

From eq. (39), the input surface reactance vanishes at wo. Therefore, 
since the corrugated waveguide is connected at the input to a smooth 
waveguide of the same diameter, the input reflection is essentially zero 
at w = wo. For w ¥ wo, however, there is a reflection which, as we shall 
see in Section V, is given accurately by 


|p|? = (Goo ae (42) 
Bi + By 
where {; is the propagation constant of the TE;, mode and 8; is that of 
the HE); mode. Equation (42) is derived in Ref. 4. 
At frequencies higher than we, some of the incident power is converted 
into the HE;,; mode. If P, denotes the converted power and Pp denotes 
the incident power, then 


Pe” Lath —oo| ut 
Po La(B;— Bs)yJ} Ba(u? — 1)’ 
where u = 1.841 and 8; is the propagation constant of the HE, mode. 
Some of the incident power is also converted into the TM; mode of the 
smooth waveguide. If P.. denotes this power, which is reflected by the 
junction, we have4 


(43) 


P. = (Bz = Bo)? Bobi i (44) 
Po (81 + B2)?u® — 1k?” 

where £2 is the propagation constant of the TM,; mode. The above 

coupling equations are derived in Ref. 4 assuming |y| > 1. 

In the experiment described in the following section, the corrugated 
waveguide is connected at the input to a smooth waveguide whose radius 
gradually increases from a relatively small value a’ to the final value a(0). 
The initial value a’ is sufficiently small so that the TM; mode is cut off 
and, as a consequence, the power P, is reflected back towards the junc- 
tion where it is converted into the HE;, mode of the corrugated wave- 
guide. The total power converted into the HE; mode is thus in general 
different from P,. It varies approximately between the two values 


(VP, — VP.) and (VP, + VP;)?, (45) 
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Fig. 3—Return loss of corrugated horn of Fig. 1 vs frequency. 


depending on the phase angle of the TM, mode incident on the junc- 
tion. 


V. EXPERIMENT 


A corrugated feed of the dimensions shown in Fig. 2 was built using 
a technique described in Appendix C. Its measured reflection coefficient 
agrees* very well with eq. (42), as shown in Fig. 3. Its radiation charac- 
teristics are shown in Figs. 4 through 6. From 17.5 GHz to 32 to GHz, the 
radiation patterns are in good agreement with eq. (33), as shown in Fig. 
4, which compares eq. (33) with two measured patterns of |F,.|? in the 
plane ¢ = 45°. Furthermore, from 17.5 GHz to 32 GHz, there is little 
difference between the patterns of |, |? in the two principal planes (¢ 
= 0 and ¢ = 90°) and the pattern for ¢ = 45°. The difference is altogether 
negligible at wo = 23 GHz, which is the frequency for which y = 0 at the 
aperture. Figures 5 and 6 show a few examples of patterns measured in 
the two principal planes. The variations with frequency of the beam- 
widths 26; and 262 (respectively, the 3-dB and 10-dB beamwidths) agree 
very well with eqs. (34) and (35), as shown in Fig. 7, where the measured 
beamwidths in the plane ¢ = 45° are compared with the calculated 
values. 

Finally, the cross-polarized component E, is very small over the fre- 
quency range of eq. (2), as shown in Fig. 8. For w > 28.8 GHz, however, 


* Note that the measured reflection coefficient includes a small reflection due to a 
transition from rectangular to circular waveguide, which was connected at the input of 
the feed during the measurement. 
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Fig. 4—Calculated and measured radiation patterns vs normalized angle at ¢ = 45°. 


the normalized peak C? of Ey increases rapidly with frequency as ex- 
pected because of the HE), mode, which is excited at the input for w > 
28.8 GHz. The solid curves and the dashed curves in Fig. 8 correspond 
to eqs. (27) and (48), respectively. The agreement with the measurements 
is satisfactory, taking into account that eq. (43) is not expected to give 
exactly the total power converted into the HE,, mode, for several reasons. 
In the first place, eq. (43) assumes a very large number of disks per 
wavelength, 


hd, 
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whereas in the experiment, h ~ 0.137) at 30 GHz; the effect of a finite 
number of teeth is briefly discussed in Appendix B. In the second place, 
eq. (43) assumes that only the TE,; mode is incident at the input, 
whereas in practice also the TM); mode is incident, for the reason 
pointed out in Section IV. The total power carried by the TM,,; mode 
is approximately given by eq. (44). 


DECIBELS 








BEAMWIDTH (@) IN DEGREES 


Fig. 5—Measured radiation patterns at ¢ = 0° and ¢ = 90° (principal planes). 


Note, finally, that even if the total power converted into the HE}, mode 
is calculated accurately, to determine the resulting cross-polarization 
peak C2, the difference in phase between the HE;; and HE;, modes must 
be determined. 
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Fig. 6—Measured radiation pattern at ¢ = 0° and ¢ = 90° for 32 GHz (effect of HE,; mode 
on pattern symmetry is evident at upper frequency limit of operation). 


VI. CONCLUSIONS 


Using a novel fabrication technique, described in Appendix C, which 
can be applied at very high frequencies, a corrugated feed of small flare 
angle was fabricated. Its input reflection, found to be given accurately 


by the simple formula 
Bi — Bi 
ial (ew 
lal? = Bi + By 


remained less than —30 dB from about 24 to 32 GHz. Over the frequency 
range w, < w < wo, the far field was essentially polarized in one direction; 
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Fig. 7—Calculated and measured 3-dB and 10-dB beamwidths. 


its pattern is simply given by 


Jo(v) \2 
PO (uP a gi)” 
At frequencies higher than 28.8 GHz, the far field contained a cross- 
polarized component, as predicted by the theory of Ref. 4. 

It was shown that a maximum bandwidth, expressed as the ratio w/w, 
of about 1.68 can be achieved for the type of corrugations considered 
here. By using the corrugations of Ref. 10, which, however, are difficult 
to realize at high frequency, greater values of w/w; may be achieved. 

Both the input reflection and the cross-polarization ratio can be im- 
proved by increasing the thickness t of the disks at the input. Curves of 
|p|? and C? as a function of frequency for different values of t/h are given 
in Fig. 2 of Ref. 4. In the experiment, t/h was kept constant for reasons 
of simplicity, since our main concern was to verify the results of Sections 
II through IV and of Refs. 1 and 4, and to demonstrate the feasibility of 
the fabrication technique described in Appendix D. 
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Fig. 8—Calculated and measured cross polarization; these are maximum values of cross 
polarization (see Section III, eq. 26). 
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APPENDIX A 
Cutoff Frequencies of the HE,,, and HE is Modes 


To determine eqs. (37) to (40) it is necessary to establish the relation 
between a, b and the cutoff frequencies of the various modes of eqs. (12) 
and (13). Assume the thickness of the disks is very small. Consider first 
the cutoff frequencies w., of the HE}, modes, which are characterized 
for large ka by the field distribution of eq. (12). In the vicinity of the 
cutoff frequency wem, we have, for any one of the above modes, 


cos 6, = 0, (46) 


where k cos 6, is the propagation constant in the z direction, 


Un 
COs Om =/1 = aye (47) 
For w~ wem, we can show that the mode degenerates into either a TE 
mode or a TM mode according to the following rule: 
TE mode, if y > 0, | 


48 
TM mode, if y < 0, eo 
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Fig. 9—Region where only the HE; mode propagates. Note: As w varies for a given a,b, 
a point ka,kb moves along a straight line through the origin. 


In the former case, since the only nonzero component of the magnetic 
field of a TE mode near cutoff is H,, we have Hy = 0, and therefore the 
second of the two boundary conditions (6) can be ignored, which implies 
that wm is independent of y. The corrugated waveguide can be replaced, 
therefore, with a smooth waveguide of the same diameter whose cutoff 
frequencies for the TE modes are given by the roots of J;(ka), 


J (ka) = 0 for w = wem (y > 0). (49) 
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Fig. 10—Aluminum and brass disks assembled on mandrel. (a) Before machining. (b) After 
machining exterior surface. (c) With electroformed wall of copper. (d) After machining 
interior surface. (e) After etching away aluminum disks. 


In the latter case, y < 0, the only component of the electric field for 
w ~ Wem is E, and, therefore, the very thin disks can be removed without 
affecting the field. The cutoff frequencies therefore coincide with those 


of the TM modes in a smooth waveguide of radius b, and are given by 
the roots of J (kb), 


J (kb) = 0 for w = Wem, (y < 0). (50) 


Analogous considerations? are valid for the HE,,modes, except that 
now, instead of the rules (48), 


TM mode, if y > 0, 
TE mode, if y < 0,}’ 
and therefore the cutoff frequencies w,, are given by the conditions 
J (ka) = 0 for w = Wem if y <0 


(51) 


: (52) 
J, (kb) = 0 for & = wom, if y > 0. 
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Fig. 11—Modification of assembly to obtain strong attachment to teeth of electroformed 
wall. 


These results are illustrated by the diagram of Fig. 9 showing in the ka,kb 
plane the region* where only the HE,;-mode can propagate. We can see 
that w/w, is maximum when b/a = 7.0155/3.8317 = 1.8309, in which case 
w/w, = 1.6839, as pointed out in Section I before eq. 2. 


APPENDIX B 
Effect of a Finite Number of Teeth Per Wavelength 
In Section IV, we assumed an infinite number of teeth per wave- 
length—i.e., 
h<, 


a condition which seldom holds in practice. For instance, in the experi- 
ment 


h > 0.18716 


for w > 30 GHz. The effect of a finite number of teeth is not difficult to 
evaluate approximately if 


ka> 1, 


In this case, in fact, the field in the vicinity of the corrugations can be 
considered to be locally the field of a plane wave reflected by a corrugated 
plane tangent to the actual corrugated surface. The angle of incidence 


* In Fig. 9 we have not indicated a region for kb < 3.8317 where mode propagation can 
occur, since that region is of no interest to us here. (The HE); mode is cut off for kb < 
3.8317). 
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is 90° — 6,, where 
B 
cos 6, = 


Thus, from Ref. 9, we find that a corrugated waveguide of radius a > 
\ with a finite number of teeth is equivalent to one with h « X, but with 
slightly larger inner radius a’, 


In2 
eae hy. 
Tv 


and the same outer radius b. Although this result is strictly valid only 
if a >> A, measurements!! have shown that it is valid approximately even 
for ka ~ 4. 


APPENDIX C 
Fabrication Technique 
Even if the aperture of a feed is illuminated by a single mode, the far 
field contains a cross-polarized component with amplitude proportional 
to [eqs. (4), (5), and (27)] 
| Ne st Bs! 
ka (1 — t/h) ka tankl| 
To minimize this cross-polarized component, it is important that the 
thickness ¢ of the disks be much smaller than their separation h, 


t<h. (54) 


(53) 


Since h is always appreciably smaller than \/4, and since typically the 
depth / of the grooves is not much different from \/4, condition (54) 
implies 


t <« r/4 (55) 
t<«l. (56) 


Corrugated feeds are difficult to fabricate. When a corrugated feed 
is electroformed, a mandrel of aluminum or other material is first pre- 
pared, and then the corrugated feed is electroformed around the mandril, 
which is then removed with a solvent. However, at high frequencies, say 
at frequencies higher than about 10 GHz, condition (55) demands that 
t be very small (much less than 0.318 cm). Then, taking into account (56), 
the above technique cannot be used. 

Figures 10 and 11 illustrate a technique that can be used at very high 
frequencies, as high as 100 GHz, and which allows very small thicknesses 
t to be realized. First a set of disks of aluminum and brass is assembled, 
as shown in Fig. 10a, to form a single block whose outside surface is then 
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machined as shown in Fig. 10b. The central rod shown in Figs. 10a and 
b defines the axis of the corrugated waveguide. The surface S; corre- 
sponds to the bottom of the grooves, as shown in Fig. 10e. The thickness 
of the brass disks is ¢ and, that of the aluminum disks, h — t. After the 
surface S; is machined, a wall of copper (or other metal) is electro- 
formed, as shown in Fig. 10c. The central rod is then removed, and the 
inside surface Sg, which corresponds to the tops of the grooves (see Fig. 
10e), is machined as shown in Fig. 10d. Finally, the aluminum is removed 
with a solvent and the final product is the corrugated waveguide of Fig. 
10e. An important feature of this technique is that the two surfaces S; 
and Sg can be machined very accurately. 

If the brass disks are too thin, their mechanical adherence to the 
copper wall in Fig. 10e may not be satisfactory. In this case, we may. 
modify the above procedure by adding, after step b, an extra step b’ in 
which the brass disks of Fig. 10b are replaced with disks of somewhat 
larger diameter, as shown in Fig. 11b’. The exact dimensions of the new 
disks are not important. The final result, after the aluminum is removed, 
is shown in Fig. 11e’. The disks are now embedded in the copper wall. 
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Acoustic Properties of Longitudinal 
Displacement in Vocal Cord Vibration 
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We examine the acoustic significance of longitudinal displacement 
in the self-oscillatory behavior of the vocal cords, and inquire into the 
need for representing this detail in speech synthesis. We use computer 
techniques and a previously derived model of the vocal cords to study 
the contribution of longitudinal displacement to the total acoustic 
volume velocity generated at the vocal cords. This volume velocity is 
the effective sound source for production of voiced speech. From com- 
putational results, and from speech sounds synthesized by the pro- 
grammed model, we find that the contribution of longitudinal dis- 
placement is not significant perceptually, and is not essential for 
modeling the dominant acoustic properties of voiced speech. 


I. VOCAL-CORD MODEL 


In earlier work!? we derived an analytical model for the self-oscillatory 
motion of the human vocal cords. We consider the displacing tissue of 
each cord to be approximated by two stiffness-coupled masses (see Fig. 
1). For normal (nonpathological) conditions of phonation, the oscillator 
is bilaterally symmetric, and the mechanical constants of the opposing 
cords are identical. The left-hand mass pair (denoted m, m}) constitutes 
the bulk of the firm cord tissue, while the smaller right-hand mass pair 
(mo, ms) represents the more flaccid mucous membrane covering of the 
firmer tissue. Each mass has associated with it a restoring stiffness and 
a resistive loss. All the stiffnesses and resistances are substantially 
nonlinear,! and in the original work, these elements act to oppose lateral 
motion (x-direction) only. The restriction to lateral motion still permits, 
of course, phase differences in the motion of the coupled masses. Lateral 
displacement of each mass pair determines the cross-sectional area of 
opening at each position. If the length of the cords, or glottal opening, 
is taken as ¢,, then the cross-sectional glottal areas are taken as rec- 
tangular shapes whose areas are Ay; = 2f,x;, 1 = 1, 2, where the factor 


889 


VOCAL TRACT 





Fig. 1—'T'wo-mass model of the vocal cords. Translational displacement is permitted in 
lateral (x) and longitudinal (y) directions. 


2 arises from the bilaterally symmetric cord configuration. These 
cross-sectional areas determine the acoustic properties of the glottal 
volume current U,,, which enters the cord orifice (from the subglottal 
system), and that which leaves it Uz¢ (to pass into the larynx tube). The 
latter volume velocity is the effective sound source for all voiced speech 
sounds. The air pressure just to the left of (beneath) the vocal cords is 
the subglottal pressure P,,, and the pressure just to the right of (above) 
the cords, at the entrance to the vocal tract, is P;. The differential 
pressure (P,, — P;) is the potential that creates the glottal volume cur- 
rents. 

The resulting volume currents depend upon serial acoustic impedances 
dictated by A, and A, and, hence, upon the cord motion, which, in turn, 
is conditioned by the intraglottal pressure distribution in the orifice and 
by the transglottal pressure (Ps, — P;). These serial acoustic impedances 
also are nonlinear (and flow dependent), and represent the mass (in- 
ertance) of air contained within the glottal orifice and the associated 
resistive flow losses.! 

Additionally, there is another potential influence upon the glottal flow, 
namely, the volume of air displaced by the vibrating mass pairs. In 
general, this volume displacement can have components associated with 
lateral and longitudinal motion. In the original work, components of 
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Fig. 2—Equivalent acoustic circuit including lateral and longitudinal displacement-volume 
velocities. 


Table | — Values of impedance components of Fig. 2 




















2 p [Us| _?P ( 1 "r) 
Re a Rio 5 Az, Az, | Ugsl 
Serial 
bgt (i-4) 
Impedances Re D AgsAy 1- |Ug| 
Ro. = 12ué2d,/A3,, Ryo = L2ufedo/A ge 
Lg = Pry ee Lega = pd o/ fAss 
i eat 2 dy - dy 
Longitudinal Uy, = 2€g(d) + do) =, Uye = 28 g(d1 + do) = 
Components d dt 
d 
Ux = dy a (Agi), Ux2 = do = 7 (Aga) 
= 20d; d,,/dt, = 20,do dxo/dt 
Lateral 
Components Cg1 = Agidi/pc? Cgo = Agods/pc? 
G1. = Si) 2 = Ro. Gp= SIV Awo 
2cpp pe 2cpp 
Se ae So = 2(lg + 2xo)de 
p = 1.14 X 10-8gm/cm’, air density 
Constants pw = 1.86 X 10-4dyne-sec/cm2, kinematic-coefficient of viscosity. 
(for vocal ; 
system, c = 3.5 X 104cm/sec, sound velocity 
moist air n = 1.4, adiabatic constant 
ee d = 0.055 X 1073cal/cm-sec-deg, coefficient of heat conduction 
ture)* wo = 27(1000), mid audio range radian frequency 


Cp = 0.24 cal/gm-degree, specific heat at constant pressure. 


* From J. L. Flanagan, Speech Analysis, Synthesis and Perception, second edition, 
New York: Springer Verlag, 1972. 


glottal current corresponding to rate of volume displacement (both 
lateral and longitudinal) were neglected. 


ll. ACOUSTIC CIRCUIT 


Recognizing that the cord dimensions are very small compared to 
sound wavelengths at the frequencies of interest, and that all mechanical 
velocities are small compared to the sound velocity, we derive a one- 
dimensional equivalent circuit for the acoustic quantities involved. Its 
complete form is shown in Fig. 2. The values of all impedance elements 
are given in Table I. 

The serial elements (top branch in Fig. 2) are identical to those of our 
original work!, and relate to time-variation of the acoustic impedance 
of the glottal opening. 
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Fig.3—Simplified equivalent acoustic circuit, including longitudinal displacement 
currents. 


All shunt elements relate to rate of displacement of air volume by the 
moving cord masses. The time variation of all shunt quantities is also 
determined by the motion of the cord masses. 

Lateral motion of the cord masses (normal to the direction of glottal 
flow) displaces air volume at the rate of 


U,; = 20, go cm3/s, i = 1, 2, 
dt 
where x; is the lateral displacement and d; is the depth (thickness) of 
the cord element (mass). Again, the factor 2 arises from the two bilat- 
erally opposing cords. The acoustic compliances, C; and C2, represent 
the compressibility of the small air volumes contained between the op- 
posing cords and the conductances G, and G2 represent the heat-con- 
duction loss at the tissue surfaces of the cords. 

Longitudinal motion of the cord masses is assumed to occur cophas- 
ically and to be translational only. In this regard, consider the y- motion 
of the locked masses to be opposed by a nonlinear spring and loss similar 
to that of k; and r,. The effective surface area exposed to the transglottal 
pressure difference is taken to be the product of cord length and total 
cord thickness, ¢,(d; + dg). No cavity compliances or losses are associ- 
ated with the longitudinal motion, and the longitudinal contribution to 
the total volume velocity is 


d 
Uy; = 2€g(d; + dy) i= 1,2. 


In other words, Uy; and Uy are equal and oppositely poled. 

Notice that in the earlier formulation, the absence of the shunt el- 
ements imposes the constraint U,¢ = Ug, = Ug. The presence of the 
shunt elements (all time-varying with displacements that are determined 
by the equations of motion for the mechanical system which, in turn, is 
forced by the intraglottal and transglottal pressures to close the feedback 
loop of the oscillator) makes the input flow U,; and the output Uge 
typically different. 
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Fig. 4—Computed mechanical behavior of the vocal-cord/vocal-tract model. The vowel 
configuration is //. 


A recent related study? examined the influence of the U,; upon Uge. 
The present study considers separately the effects of the U,;. For this 
purpose, the circuit of Fig. 3 is a simplification of Fig. 2. 


lll. PURPOSE OF PRESENT STUDY 


In the original work,!? we made estimates of the volume displacement 
currents, based upon long-wave assumptions and one-dimensional sound 
propagation, together with what we believed to be reasonable physio- 
logical estimates of cord velocities (compared with volume velocities 
responsive to transglottal pressure). We concluded that displacement 
currents are of second order, and in the original work we chose to neglect 
them in favor of elucidating dominant principles. The original formu- 
lation, therefore, treated only lateral displacement as it affects the serial 
glottal impedances. As a matter of completeness, we more recently have 
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Fig. 5—X-Y trajectories for initiation of oscillation. Trajectories are for pellet positions 
shown in insert. 


returned to a quantitative examination of these assumptions. A first 
study,? now completed, considered the importance of the shunt branches 
that represent the lateral components of volume velocity generated by 
the displacing masses—that is, from the volume current sources U,; and 
U,.2. The results of the study support the original assumptions, and show 
the lateral components to be second order by comparison to the currents 
actuated by the pressure difference acting across the glottal opening. 

The present study examines the contributions of the longitudinal 
displacement to the total glottal volume velocity (specifically, the con- — 
tributions of Uy; and Uy2) and the importance of longitudinal dis- 
placement to the self-oscillatory dynamics of the cord model and to 
sound perception. 

We take the longitudinal restoring stiffness k, typically to be the same 
as the lateral restoring stiffness k,, namely 80 kdynes/cm. The longitu- 
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Fig. 6—X-Y trajectories for steady-state oscillation of the cord model. 


dinal loss (or damping ratio) is also taken to be similar to the lateral one, 
namely ¢, = ¢, = 0.2 These values are based upon clinical observations.? 
We examine these choices subsequently. Further, since the longitudinal 
and lateral motions of the cord masses are considered to be translational 
only, no rotational behavior is included. Still further, while the lateral 
translations of the coupled masses m; and m2 can have large (and 
physiologically natural) phase differences, their longitudinal translations 
are considered to be cophasic, and the internal coupling stiffness is as- 
sumed to act only for lateral motion. The lateral and longitudinal mo- 
tions are, therefore, coupled only through the acoustic variables that 
determine the oscillator forcing functions. In the course of our discussion, 
we will indicate comparisons to actual physiological data to assess the 
realism of these assumptions. 
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Fig. 7—Computed acoustic qualities for the vocal-cord/vocal-tract model. The vowel 


is /o/. 


IV. RESULTS OF COMPUTER SIMULATIONS 


The vocal-cord model, as represented by Fig. 3, was combined with 
a transmission-line formulation of the vocal tract that we have used 
previously in speech synthesis studies.4 The programmed vocal tract 
contains 20 sections which, in addition to the classical acoustic elements, 
represents the yielding soft walls of the tract and sound radiation from 
the yielding walls. This formulation is based upon measurements of 
tissue impedances that we reported earlier.° Also included for the present 
study is a transmission-line representation of the subglottal system. Six 
sections of line represent the trachea, bronchi and lungs, as previously 
described. We implemented the entire system in terms of difference 
equations programmed on a laboratory computer by techniques we have 
described in detail previously.}-? 
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Fig. 8—X-Y trajectories observed from excised larynx of dog (after Baer). 


Most of the data reported here are for the vocal tract configured in 
the shape for the neutral (schwa) vowel /a/. Some data are also included 
for the vowels /i/ and /a/. 

A first step is to ascertain if the cord oscillator, so arranged for longi- 
tudinal motion, performs realistically when compared with observations 
on the human larynx. A second step, then, is to determine the acoustic 
significance of the volume displacement current arising from longitudinal 
motion. 

Throughout these calculations, the laryngeal parameters are set to 
the “standard” values used earlier for phonation by a man’s voice in the 
chest register! (i.e., neutral glottal area Ag. = 0.05 cm?, cord tension 
parameter Q = 0.78, d; = 0.25/Q cm, dz = 0.05/Q cm). Recall that the 
Q parameter scales the values of mass and stiffness and, hence, also the 
values of the d;. Phonation is initiated by raising the lung pressure P, 
smoothly from zero to the standard value of 8 cm H20. The pressure is 
elevated in a 10-ms interval. 


4.1 Mechanical behavior 


As the lung pressure is elevated, the model commences a buildup of 
oscillation. After four or five transient swings, the oscillation settles into 
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Fig. 9—Subglottal pressure variation measured on a human subject (after Sawashima). 


a steady state behavior with a fundamental frequency (pitch) determined 
by the model parameters (the tension parameter Q has the dominant 
effect on pitch!). The initial 80 ms of this synthetic phonation is illus- 
trated in Fig. 4 for the mechanical variables. 

The top two curves show the displacements x; and x2 of mass pair m1 
and mass pair mg, respectively. The first collision of each mass pair is 
indicated by the first flat, negative-going portion of the displacement 
waveforms. For the Ag, = 0.05 cm? value, this occurs for x; = —0.0178 
cm. Note, too, that x; leads x2 in phase by the order of 60°, which is 


/a 
ky = 40 kdyne/cm {x-~y TRAJECTORIES] 


ky = 80 kdyne/cm 


DISPLACEMENT (X1072 cm) 
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0 10 20 30 40 50 60 70 80 90 
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Fig. 10—Effect of longitudinal restoring stiffness k, upon the longitudinal displacement, 
y. Data show oscillation buildup for a lung pressure P, that is raised smoothly from zero 
to 8cm H,0. 
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Fig. 11—Effect of ky upon P;. 


0 


consistent with observations from high-speed motion pictures of the 
human vocal cords. The third trace shows the longitudinal displacement, 
y, which bulges upward as P, is raised. The y motion is roughly sinus- 
oidal. The lower trace shows the net area of glottal opening A, (namely 
the minimum of A,; and Ago). The y-displacement is seen to lead in 
phase the A, wave, again consistent with the upward, rolling motion seen 
in high-speed photography of the real cords. 

An x-y plot of the buildup transient portrays the behavior perhaps 
more graphically. Figure 5 shows the x; vs y and the x2 vs y values with 
time as the parameter. Imagine pellets fixed to the lower and upper inner 
edges of one simulated cord, shown by the inserted anterior-posterior 
view of Fig. 5. The trajectories of the two pellets are plotted for the os- 
cillation buildup. The y-axis is broken and re-originated at y = (d; + 
dy) = [(2.5 + 0.5mm)/Q] = 3.8 mm. The flat portion of the tracks, along 
the vertical midline, reflect collision with the opposing vocal-cord 
mass. 

After several initial swings, the oscillator settles into a steady-state 
behavior. One cycle of this trajectory is shown in Fig. 6. The steady-state 
pitch frequency in this case is Fo = 125 Hz, or a period of T = 8 ms. 


4.2 Acoustical behavior 
The corresponding acoustical parameters, calculated for the same 
buildup period, are shown in Fig. 7. The acoustic pressure at the input 
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Fig. 12(a)—Steady-state X-Y trajectory for k, = 40 kdyne/cm. 


to the vocal tract P; is shown in the top trace. It reflects strongly the 
eigenfrequency structure of the tract, in this case configured for /a/ and 
having formant frequencies of approximately 500 Hz, 1500 Hz, 2500 Hz 
... The transglottal pressure (Ps, — P;), which is the forcing function 
for the y-motion and the pressure potential for the volume flow through 
the glottal opening, exhibits a pronounced pitch-synchronous variation. 
Its peak values, in fact, approach twice the lung pressure value of P, = 
8 cm H20. Recall that P, is the lung pressure input to the simulated 
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Fig. 12(b)—Steady-state X-Y trajectory for ky = 80 kdyne/cm. 


subglottal system, representing trachea and bronchi. But notice that 
the mechanical y-displacement (Fig. 4) does not respond with this detail. 
(Neither do the x; and x2 displacements respond to high-frequency detail 
in their forcing functions—that is, the mechanical system, being mass- 
controlled, filters out this detail.) 

The subglottal pressure P., (the pressure just beneath the vocal cords) 
also exhibits a pitch-synchronous fluctuation, but of somewhat less 
amplitude, namely about +20 to +30 percent of the mean subglottal 
pressure. Its positive peaks correspond to the closing epochs of the glottal 
port. The calculated volume velocity passing the glottal opening U, 
(bottom trace) appears as a traditionally shaped, pulsive waveform. This 
wave is similar to that calculated in previous work (without longitudinal 
motion) but differs in that its values are modified by the effects that U,; 
and U,» couple into the pressure variables. That is, Uy; and Uy can 
influence P,, and P; and, hence, U,¢. The latter three variables, in turn, 
close the oscillator feedback loop by constituting the forcing functions 
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Fig. 12(c)—Steady-state X-Y trajectory for ky = 120 kdyne/cm. 


for the lateral displacement. As was the case in the mechanical variables, 
the U, flow does not reflect a temporal fine structure comparable, say, 
to the (P,, — P:) waveform. The resistive and inertive components of 
the glottal impedance (i.e., the serial components in Fig. 3) act effectively 
as a low-pass filter. It is not unusual, however, to see pronounced tem- 
poral structure that corresponds to the lowest eigenfrequency of the vocal 
tract, especially for low, back vowels (such as /a/), or for tightly articu- 
lated sounds. 

A next question, then, is how do these mechanical and acoustical 
quantities, resulting from the model with longitudinal displacement, 
compare with physiological data. 


4.3 Comparisons to physiological observations 


One qualitative comparison can be made for the mechanical dis- 
placement behavior. Baer’ performed studies on the excised larynx of 
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Fig. 13—Oscillation buildup without subglottal system. 


a dog in which he fixed pellets to the displacing tissue and made optical 
observations under stroboscopic illumination. While his pellet positions 
do not correspond exactly to our mass-pair corners, we can roughly 
compare his observations with our data. Figure 8 shows x—y trajectories 
for one set of conditions for the dog larynx that approximates values used 
in human phonation (namely P, = 8 cm H2O, Ug = 275 cm?/s, and Fo 
= 100Hz). Particles (pellets) 2 and 3 are of interest. While the vibratory 
excursions of the excised dog larynx are larger than those we calculate 
with the model, the qualitative motions are gratifyingly similar. One 
question that arises is how much does the longitudinal (vertical) dis- 
placement depend upon the choice of longitudinal stiffness constant. 
We shall examine this question in more detail subsequently. 

Another comparison can be made in the acoustic domain—namely, 
to the subglottal pressure variation P,, shown previously in Fig. 7. Sa- 
washima® has measured the subglottal pressure during phonation in a 
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Fig. 14—Steady-state oscillation without subglottal system. 


human subject. One of his results is shown in Fig. 9. The qualitative 
correspondence to the model calculation appears relatively good, and 
the acoustic interaction among the simulated cords, vocal tract, and 
subglottal system is realistic. 


4.4 Effect of longitudinal stiffness constant 


In view of uncertainties in the measurement of the stiffness constants 
in physiological preparations, it is important to examine how critical the 
value of k, (the longitudinal restoring stiffness) is to the oscillatory be- 
havior of the model. 

For the bulk of our studies, we have taken ky equal to our typical 
“standard” value of k;, namely 80 kdynes/cm.! We have also used the 
standard value for the damping ratio, {, = {) = 0.2. This choice is based 
upon the physiological measurements on cord tissue conducted by Ka- 
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Fig. 15—Effect of vowel configuration upon longitudinal displacement. 


neko,? who found the stiffness constants for lateral and for longitudinal 
displacement to be similar. Both stiffnesses are also taken to have the 
same nonlinearity, namely, a cubic nonlinearity in the restoring force 
(see Ref. 1). 

To assess the model’s sensitivity to large variations in ky, we calculated 
the buildup of synthetic phonation for k, = 40, 80, and 120 kdynes/cm. 
The resulting y-displacements for these values is shown in Fig. 10, and 
the sound pressure at the entrance to the vocal tract is shown in Fig. 11. 
Also, the x-y trajectories are shown in Fig. 12a, b, and c. 

As would be expected, the greatest influence in this variation is re- 
flected in the y-displacements. The “softer” k, gives larger de dis- 
placement and smaller peak-to-peak vibratory excursions. The P; data 
indicate that the variations in acoustic behavior and glottal excitation 
are very small. The fundamental pitch is sensibly the same for all cases, 
namely 125 Hz. This factor is almost completely dominated by the lateral 
motion. In auditory assessment of the output synthetic sound, the dif- 
ferences are virtually imperceptible, suggesting that the longitudinal 
displacement current is insignificant for speech synthesis. 


4.5 Effect of subglottal system 


A side issue, of some interest in passing, is the effect of acoustic in- 
teraction between the subglottal system and the cord oscillator. If the 
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Fig. 16(a)—X-Y trajectory for the vowel /a/. 


trachea-bronchi-lung system is removed and the lung pressure applied 
directly to the cord input (as a zero-impedance source; i.e., Psz becomes 
a pressure “battery” equal to P,), then the temporal structure previously 
reflected in P,,. is eliminated, the transglottal pressure excursions simply 
equal P,, and the longitudinal component of cord displacement is less- 
ened. This is illustrated by the x-y trajectories for oscillation buildup 
and steady state shown in Figs. 13 and 14. For this case, k, is reset to the 
typical value 80 kdynes/cm. Note the slight lowering of the fundamental 
frequency to 123 Hz. 


4.6 Effect of vowel configuration 


It also is instructive to consider the influence of vowel configuration 
upon the cord model, as presently formulated. Such studies were made 
in detail in the original work.! We therefore compare synthesized results 
for the vowel /a, a and i/. (In this case the longitudinal stiffness constant 
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Fig. 16(b)—X-Y trajectory for the vowel /a/. 


was set to k, = 88 kdynes/cm through an inadvertent keypunch error. 
Because the differences are so small, it did not seem worthwhile to re- 
compute the data for k, = 80 kdynes/cm.) 

The y- displacements are shown in Fig. 15, and the x-y trajectories for 
one period of steady-state oscillation are shown in Fig. 16a, b, and c. The 
longitudinal displacement is not greatly affected by vowel configuration, 
but the constricted articulations /a, i/ clearly lead to slightly greater 
longitudinal peak-to-peak excursions than does the open-pipe (neutral) 
vowel /9/. This typically is owing to the greater acoustic interaction at 
the eigenfrequencies for the configurations with higher acoustic im- 
pedance levels, which in turn leads to greater transglottal pressure dif- 
ferences. This is well reflected in the acoustic variables resulting from 
this calculation. The corresponding acoustic quantities are shown in Figs. 
17 thru 20. In these data, note especially how the tract eigenfrequencies 
are manifest, including in the synthetic output sound pressure from the 
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Fig. 16(c)—X-Y trajectory for the vowel /i/. 


mouth, Pout. Note, too, that the resulting steady-state pitch frequencies 
are about 125 Hz for /a/ and /a/ and about 120 Hz for /i/. 


V. SIGNIFICANCE OF LONGITUDINAL DISPLACEMENT 


Having established that the cord oscillator (with lateral and longitu- 
dinal degrees of freedom) appears to behave realistically, we consider 
the next questions: 

(t) Is the longitudinal motion significant or necessary for proper 
self-oscillatory operation of the model? 

(ti) Is the acoustic volume velocity contributed by the longitudinal 
motion physically or perceptually significant? 

So far as the purposes of speech synthesis are concerned, we answer both 
of these questions in the negative. 

What we wish to do, therefore, is compare the mechanical and 
acoustical behavior with and without longitudinal displacement. Because 
the longitudinal effects are coupled only through U,; and Uy», the lat- 
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Fig. 17—Uye for the vowels /9/, /a/, and /i/. 
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Fig. 18—P,, for the vowels /9/, /a/, and /i/. 
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Fig. 19—P;, for the vowels /a/, /a/, and /i/. 
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Fig. 20—Pout for the vowels /9/, /a/, and /i/. 
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Fig. 21—Computed mechanical behavior of the cord/tract model without y-displace- 
ment. 


eral-motion-only condition (that is, the original formulation of the 
model) is conveniently realized by setting Uy; = Uy: = 0. This condition, 
with no longitudinal displacement flow, yields the mechanical results 
shown in Fig. 21. These results are essentially identical (at least to three 
significant figures) to the corresponding quantities in Fig. 4. In particular, 
note that the crucial A, waveforms are virtually identical. 

We can now examine pertinent acoustic quantities. Calculation of the 
conditions with y- displacement yields the result shown in Fig. 22. The 
figure shows one cycle of the steady-state oscillation. Recall that Ugg is 
the total volume velocity at the larynx tube entry to the vocal tract. U, 
is the flow component through the actual opening of the glottis. Note 
that U,¢ is non-zero during the time the cords are actually closed, cor- 
responding to an upward (vertical, longitudinal) displacement of air 
volume that adds positively to Ug. Similarly, later in the cycle, downward 
longitudinal displacement subtracts from U,. The difference between 
these volume velocities is 


(Uge - U, ) = Uy = Uy 


by virtue of the assumption of cophasic longitudinal motion. This dif- 
ference is plotted on an X10 enlarged scale in Fig. 22. The peak value 
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Fig. 22—Glottal volume velocities calculated with y-displacement (refer to Fig. 3). 


of the difference for this condition is on the order of 10 cm?/s, or about 
9 of the peak value of the total Uz. This result is not just peculiar to 
this range of volume velocity, but rather it scales comparably at louder 
and softer phonation. For example, if the lung pressure P, is doubled 
say to 16 cm H20, the longitudinal displacement current increases be- 
cause the transglottal pressure and the longitudinal displacement in- 
crease. But the U, flow also increases and remains far and away the 
dominant quantity. 

The amplitude spectra of these quantities provide convenient corre- 
lation with auditory percepts. The spectra for Ugz and U, are shown in 
Figs. 23a and b. A close comparison shows the differences to be less than 
2 dB, an amount that is not significant perceptually. The more relevant © 
comparison is obtained when the effect of y-motion is eliminated (by 
removing Uy; and Uy). The corresponding glottal waveform for no y- 
displacement is illustrated in Fig. 24. It is denoted U,?. Also reproduced 
is the U,¢ with y-displacement. Further, the difference between the 
longitudinal displacement and lateral-only conditions (U,¢ — Ug?) is 
shown on an-X10 enlarged scale. During the glottis-closed time, this 
difference is identical to the (U,¢ — Uz) difference of Fig. 22, because 
U, = 0. During the glottis-open time, the (Ug~ — Ug?) difference differs 
from the (U,~¢ — Ux) difference. In other words, U, differs from Ug? 
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Fig. 23(a)—Amplitude spectrum of U,¢, which includes the effects of y- displacement. 


essentially by the influence that U,; and Uys have upon the transglottal 
pressure difference (Ps, — P:). 

Again, the more perceptually relevant comparison is to the amplitude 
spectrum. The spectrum of U,, is given in Fig. 25. A close comparison 
to the Uz¢ spectrum of Fig. 23 shows the differences to be less than about 
2 dB. Auditory assessment of the output synthetic vowels shows them 
to be indistinguishable even in close comparison. 


VI. CONCLUSION 


In view of these results, we conclude that realistic acoustic behavior 
(which is needed in speech synthesis) can be obtained in the cord model 
without the additional complexity of longitudinal displacement. Lon- 
gitudinal displacement is not necessary for realistic self-oscillation of 
the model. The important vertical phase differences in the two-mass 
motion are adequately duplicated by lateral displacement only, as is the 
significant acoustic interaction between vocal tract and vocal cords. 
Further, the rate of volume displacement owing to longitudinal motion 
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Fig. 23(b)—Amplitude spectrum of U, calculated with y-displacement. 


is clearly perceptually not significant and need not be represented with 
added detail. 

These conclusions about the mechanical and acoustic behavior have 
a corollary in a companion study on the rate of displacement of air vol- 
ume owing to lateral motion only.? This contribution was examined by 
making use of the shunt branches in Fig. 2 that include U,;, Uy. Cal- 
culations and computer simulations showed that the contribution to 
glottal volume velocity of the air extruded from the glottal port by lateral 
tissue displacement is barely discriminable in a differential auditory 
comparison. In fact, the perceptual effect for the lateral volume dis- 
placement is just slightly larger than for the longitudinal displacement. 
Both are quite second-order in importance. 

We have found in the present study that proper acoustic and oscilla- 
tory behavior of the model does not depend significantly upon longitu- 
dinal displacement. The longitudinal motion is relatively insensitive to 
acoustic loading and to changes in longitudinal stiffness. The longitu- 
dinal motion influences fundamental frequency only slightly. What, 
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Fig. 24—Comparison of waveforms for U,¢, which includes y-displacement current, 
and U,%, which is calculated without y-displacement. 


then, are the critical and sensitive parameters of the cord model? In other 
words, what parameters are most influential upon the perceptual at- 
tributes of U,2, since the end product—the output sound—depends 
directly upon U,¢? The results of our earlier work can be combined with 
the insights obtained here to consider this question. 

The original study showed that the intra-glottal pressure distribution, 
and the fluid flow laws used to deduce it, are quite important to proper 
oscillatory behavior, to proper generation of the Ug, flow, and to realistic 
acoustic interaction between the vocal tract and vocal cords. To a large 
extent this pressure distribution determines how the pitch frequency 
varies with subglottal pressure and with articulatory configuration. The 
mass-stiffness product (i.e., the natural frequency of the mechanical 
system) is quite dominant in determining pitch range. Subglottal pres- 
sure, assuming it to be above an initiation threshold of several cm H2O, 
is primarily correlated with sound intensity, a relatively noncritical factor 
for voiced sounds. Mechanical parameters such as cord thickness, 
damping ratio, and nonlinearity are relatively noncritical except as they 
influence duty factor and “flow chopping” at collision (which yields a 
broad-spectrum U,¢ function). None of the mechanical variables, lateral 
or longitudinal, reflects the temporal fine structure of the acoustic 
variables, but both must and do reflect the open-close cycles of the vi- 
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Fig. 25—Amplitude spectrum of U,?, the glottal volume velocity without y-displace- 
ment. 


brating cords. The realistic phase differences in motion of the upper and 
lower edges of the cords (mg and m, in the model) allow phonation 
smoothly over a wide range of input impedances to the vocal tract (both 
inductive and capacitive), and this behavior can be obtained satisfac- 
torily by permitting lateral displacement only of the stiffness-coupled 
masses. The computational complexity of anything more detailed does 
not seem necessary from the standpoint of duplicating realistic acoustic 
behavior, which is the objective in speech synthesis. 

On the other hand, if the objective were a detailed study of tissue de- 
formation (as might be the case in simulations for clinical diagnosis or 
for representing pathological conditions) then the computational com- 
plexity of longitudinal displacement might be considered. In such a case, 
the vocal-cord model should be treated as a more distributed system. 
For the representation and synthesis of normal speech, however, these 
details do not appear perceptually significant and are not needed to 
represent the dominant properties of vocal-cord vibration. 
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Faulty-Trunk Detection Algorithms Using 
EADAS/ICUR Traffic Data 
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A class of algorithms for detecting abnormally short-holding-time 
trunks has been developed that utilizes individual trunk data available 
in EADAS/ICUR (Engineering and Administrative Data Acquisition 
System/Individual Circuit Usage Recorder). This data consists of a 
two-dimensional statistic that compresses the raw trunk measure- 
ments—the state of the trunk (busy or idle) sampled every 100 or 200 
seconds—into a manageable form. Because this data is essentially a 
sufficient statistic for the stochastic process used to model the (unob- 
servable) trunk state measurements, one of the algorithms developed 
is Wald’s sequential probability ratio test. Two of the algorithms de- 
veloped have been implemented in ICAN (Individual Circuit Analysis 
Program), and are currently being used to test trunks associated with 
the No. 1 crossbar, No. 5 crossbar, crossbar tandem (1XB,5XB, XBT), 
and step-by-step switching machines. The focus in this paper, however, 
is on the modeling and analysis aspects of the problem, and only slight 
attention is paid to the various trade-offs and real-world constraints 
encountered in implementing the algorithms. 


I. INTRODUCTION 


A message trunk, the basic connecting link in the switched telephone 
network, provides the communication path between switching machines 
as well as certain call setup capabilities, such as supervision, signaling, 
and ringing. For an important class of trunk faults that cause call failure, 
the trunk is released by the switching system upon customer abandon- 
ment and is again available to fail another call. As a result, a single un- 
detected faulty trunk of this type can fail a significant fraction of the 
offered attempts to a group and will characteristically have an abnor- 
mally short holding time. 

Because of their potential service impact, significant efforts have been 
made to understand and quantify the impact that such abnormally 
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short-holding-time trunks have on central office and network service.!~4 
It is now widely understood as a result of these studies that this generic 
trunk fault results in a fraction of service attempts “killed,” which is out 
of all proportion to their number in the trunk population. Consequently, 
traffic data available from new and existing traffic data-acquisition 
systems has been viewed in the light of increasing interest in trunk-fault 
detection. In particular, with the advent of the Bell System EADAS/ICUR 
(Engineering and Administrative Data Acquistion System/Individual 
Circuit Usage Recorder) system,” it was natural to ask whether the new 
individual trunk data available could be used to detect such “killer” 
trunks.* 

This paper discusses the theoretical aspects of a class of killer-trunk 
detection algorithms that utilize the individual trunk traffic data 
available in EADAS/ICUR. These algorithms were designed for, and 
practical versions of them are presently implemented in, the ICANT 
portion of the EADAS/ICUR system. We focus here, however, on the 
problem formulation, modeling, and analysis aspects of the algorithms 
without bringing in many of the diverse factors and trade-offs encoun- 
tered in the actual implementation. 

Because the holding time of a trunk affects the statistical properties 
of the trunk data in EADAS/ICUR, it is natural to formulate the killer- 
trunk detection problem as a problem in the testing of statistical hy- 
potheses. In this context our modeling effort is basically an attempt to 
precisely define the state of a trunk (normal or killer) and expose the 
relevant underlying distributions. Well-known aspects of the theory of 
hypothesis testing are then applicable and immediately suggest a number 
of different tests. Sequential tests are naturally considered since the 
EADAS/ICUR data evolve sequentially in time. Questions about the ro- 
bustness of the models, and the structure and performance of statistical 
tests, are addressed using standard analytic tools. 

The material in this paper has been organized into six major sections, 
whose content we briefly describe. After considering the data available 
in EADAS/ICUR (Section II), we proceed to model a trunk (Section 3.1), 
motivate an appropriate set of statistical hypotheses suitable to our 
problem (Section 3.2), and briefly review several classical tests for de- 
ciding between statistical hypotheses (Section 3.3). With these prelim- 
inaries out of the way, we develop individual trunk algorithms based 
solely on individual trunk data. Proceeding in a heuristic manner, we 
use the individual trunk data to “derive” an ad hoc killer-trunk-detection 


* The term killer trunk has been widely adopted in referring to a faulty switching- 
machine-accessible trunk in a group whose average holding time is substantially smaller 
than the average group holding time. 

¥ Individual circuit analysis program—a software program that analyzes much of the 
EADAS/ICUR traffic data and maintains the system data base. 
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algorithm (Section 4.1). Although the insight gained in proceeding in 
a heuristic fashion is significant, we shift our emphasis in Section 4.2 and 
rigorously derive an optimal test statistic. It is interesting to find that 
the ad hoc statistic is essentially one of two symmetric statistics which 
comprises the optimal test statistic. The relationship between these 
individual trunk statistics is further explored in Section 4.3. 

In Section V we factor grouping information (which essentially 
identifies all trunks common to a trunk group) into the picture, and 
develop detection algorithms tailored to trunks associated with the 
No. 5 crossbar switching machine. This development necessitates 
modeling the 5XB trunk-group selection procedure, and several results 
due to Forys and Messerli? are utilized. In Section VI we shift our dis- 
cussion to the performance of the 5XB group algorithms, deriving ap- 
proximate expressions for the mean statistic update and mean detection 
time in Sections 6.1 and 6.2, respectively. The paper concludes in Section 
6.3 with an approximate analysis of the false-alarm probability of the 
5XB group algorithms. 


ll. EADAS/ICUR DATA 


The structure of a killer-trunk detection algorithm is largely de- 
pendent on the type of individual trunk measurements available.* In 
EADAS/ICUR, the raw (unobservable) data consists of the state of each 
trunk (busy or idle) every 100 or 200 seconds. Fortunately, the data ac- 
cumulations available essentially summarize all the relevent information 
in the raw data. 


2.1 Switch-count and transition data 


The data available from the EADAS/ICUR system, which can be used 
to distinguish between normal and killer trunks, is obtained by sampling 
individual trunks at 100 or 200 second intervals. This data consists of 
periodic accumulations (typically hourly, two-hourly, or three-hourly) 
of both the Busy states, and the State transitions. The busy state accu- 
mulation is usually referred to as the switch count. For the 200-second 
sampling option with a one-hour accumulation period, the switch-count 
is an integer between 0 and 18. A state transition occurs whenever the 
state of a trunk (busy or idle) is different at two successive scans. For the 
200-second sampling option with a one-hour accumulation period, the 
number of state transitions is an integer between 0 and 17. 

If we denote the ith scan during an accumulation period in which m 
scans occur by x;, and let 0 and | correspond to trunk idle and trunk 


* Until very recently, almost all trunk-traffic measurements were obtained on a group 
rather than on an individual trunk basis. 
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Fig. 1—1ICUR data. 


busy, respectively, the available data may be written 


(i) n@n) = > x; (switch count) 
i=1 


m 
(ii) t(m) = © |x; — x;-1| (state transitions). 


t= 
Thus, the raw (unavailable) data in the form of a binary sequence, 
Xm = (x4,° ss “Xa ds 


is compressed into the two statistics n(m) and t(m). 

Because the holding time of a killer trunk is generally on the order 
of a few tens of seconds, it should have substantially more state transi- 
tions than a normal trunk, for a given switch count. Figure 1 illustrates 
the sampling process on both a normal and killer trunk.* For the pur- 
poses of this figure individual calls are represented by rectangles, call 
durations correspond to the width of the rectangles, and a half-hour 
accumulation period with the 200-second sampling option is used. 

We note in passing that for the 200-second sampling option, very little 
information is lost by “compressing” the raw data Xm = (X1,*++,Xm) into 
the two statistics n(m) and t(m). Thus, normal conversation lengths tend 
to be in the vicinity of 3 to 4 minutes and, hence, with the 200-second 
sampling option, we expect that only adjacent samples are significantly 


* The realizations shown in Fig. 1 are more or less typical for a 5XB trunk group witha 
mean group holding time of approximately 4 minutes operating at about 40-percent oc- 
cupancy, and having a killer trunk with a mean holding time of approximately 1 min- 
ute. 
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correlated. If {x;}”, is Markovian, then [n(m),t(m)] is almost a sufficient 
statistic ([n(m),t(m),x1,Xm]| is sufficient) for x,, (see Section 4.2). Note 
also that for a trunk in the killer state, succesive samples should be es- 
sentially independent (for both sampling options). 


2.2 Grouping information 


In addition to the switch-count and transition data available from the 
EADAS/ICUR system, we are also able to utilize a system map to identify 
(t) all trunks common to a trunk group, and (ii) the trunk-selection 
procedure* associated with the trunk group. It turns out that using this 
grouping information,’ in addition to the switch-count and transition 
data, enhances the detection potential considerably. 

Thus, we divide the class of algorithms into two types according to 
whether or not grouping information is utilized. The first type, which 
uses only switch-count and transition data is referred to as an individual 
trunk algorithm. These individual trunk algorithms are applicable to 
all trunks—including two-way trunks—independent of the type of 
switching machine they are associated with. They do however assume 
knowledge of the trunks nominal holding time. The second type of al- 
gorithm uses the grouping information in addition to the switch-count 
and transition data and is referred to as a group algorithm. Group al- 
gorithms are “tailored” to a specific kind of trunk-selection procedure 
and, hence, apply to trunk groups associated with specific switching 
machines. For the purposes of this paper, the trunk-selection procedure 
considered is random selection of idle trunks. This procedure models 
the selection procedure of trunk groups associated with the 5XB 
switching machine. Group algorithms generally apply only to one-way 
trunk groups. 


Ii. PRELIMINARY CONSIDERATIONS 

In attempting to quantify the intuitive notion that a killer trunk will 
exhibit more transitions than a normal trunk (see Fig. 1), for a given 
switch count, it is natural to consider the transition probabilities: 


Pio = Plait, = O/x, = YY 
and 


Pos = Pixiy, = 1/x; = 0}, 


* The map in EADAS/ICUR indicates the type of switching machine that the trunk group 
is associated with, and this allows us to model the trunk-selection procedures (see 
Ref. 2). 

tT We will consistently use “grouping information” to refer to both the identification 
of all trunks common to the group and the trunk-selection procedure associated with the 
group. 
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where x; denotes the state of a trunk at epoch t and 7 denotes the sam- 
pling interval. Of course, to evaluate these transition probabilities, we 
must be concrete about how we model a trunk. 

Before tying ourselves down to any specific model, however, it is useful 
to view these conditional probabilities in a canonical form. Thus, suppose 
we begin by assuming only that the binary valued process x; is stationary. 
We have then the following simple result: 


Lemma 1: Let x; be a binary valued, stationary random process and let 
p and R(-) denote its mean and covariance function, respectively. 


Then, 





(1) 


() Paolo) =») {1- Z | 


R(0) 
(tt) pPy,0(p,7) = (1 — p)Po,1(p,7). (2) 
Proof: Part (1) is a consequence of the definition of R(-). That is, 
R(t) = E(x,xX147) — p* 
= P(x, = 1x4, = 1) - p?, 
where 
p = E(x:) = P(x; = 1). 
Part (ii) follows from the two identities: 
p = pP11(p,7) + (1 — p)Po,1(0,7) 
and 
1 = Py 1(p,7) + Pi 0(p,7). 


A consequence of this result is that uncorrelatedness and independence 
are equivalent: 


Corollary 1: For the process in lemma 1, x;,x¢4, are independent if and 
only if R(r) = 0. 


Note that the dependence of R(-) on p has been suppressed for conve- 
nience. 
3.1 Modeling an individual trunk 


A particularly simple way to model a trunk is as the server in a single 
server loss* system with a Poisson arrival process and an exponential 
service time distribution. This model is commonly denoted by 
M/M/1-loss.® Let x, denote the state of the server: 


* In a loss system, customers who are blocked depart without waiting. 
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acts 1 if server is busy at epoch t 
‘ |O. if server is idle at epoch t 
and R(-) the covariance function of x;. It is easily shown’ that for the 
M/M/1-loss system, 


R(t) = R(O) exp {—(A + w)7}, 
where \ and uw are the mean arrival and service rates, respectively. Thus, 
P10(p,7) may be written as: 


Pyo(p,7) = (1 —p) | 1— exp =k (3) 


where the trunk occupancy p is equal to A/(A + »). Throughout this paper 
we will be concerned with 7 = 100 or 200 seconds and a nominal holding 
time 1/u in the vicinity of 3 minutes. The mean holding time of a killer 
trunk 1/y* will always be expressed as 1/ru with r typically in the range 
5 to 15. Thus, if we denote uz by S, we may write the transition proba- 
bility P{x,4, = 0/x, = 1} for a trunk with mean occupancy p as 


Prolor) = (1p) [1- exp (=) ], (4) 
Ra 
where r = 1 corresponds to a normal trunk. (Since P19 = 1 — p implies 
that x,,%9;,.-+ are independent, we will assume independence for r 
sufficiently large in subsequent sections.) 

Figure 2 is a plot of P19 vs p corresponding to S = 10/9 (200-second 
sampling and a 3-minute mean holding time) for several values of r. P1,o 
is essentially equal to 1 — p for r = 5. Figure 3 is a similar plot of P19 vs 
p corresponding to 100-second sampling and a 3-minute mean holding 
time (S = 5/9). In this figure P19 is essentially equal to 1 — p for r= 
yess 

Before putting too much emphasis on the transition probabilities in 
Figs. 2 and 3, it is prudent to consider the effect of factoring more realistic 
assumptions into the single server loss model. Thus, while the Poisson 
arrival process assumption is probably a reasonable assumption for a 
trunk in a 5XB trunk group (random selection of idle servers), it poorly 
models the overflow nature of the traffic offered to trunks in a 1XB/XBT 
trunk group.? In the latter case, it is more appropriate to model the input 
stream to a trunk as a peaked process.® Figure 4 is a plot of Pi vs p 
parameterized by the peakedness (z) of the input stream. This figure 
is based on an expression for Po derived for a GI/M/1-loss model! with 


1 The trunk-selection procedure for 1XB and XBT trunk groups is essentially a two- 
sided ordered hunt.? 

1 GI/M/1-loss denotes a single server loss system with a renewal process input stream 
(GI) and an exponential (M) service time distribution. 
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M/M/1—LOSS MODEL 
p71 = 180 SECONDS 
T = 200 SECONDS 


->7 KILLER TRUNKS 


2 NORMAL TRUNK 
a (r= 1) 





0 20 40 60 80 100 
MEAN TRUNK OCCUPANCY IN PERCENT 


Fig. 2— 1 — 0 transition probability for the 200-second sampling option. 


M/M/1 — LOSS MODEL 
pt-| = 180 SECONDS 
T = 100 SECONDS 


~~ >>KILLER TRUNKS 


NORMAL TRUNK 
(r= 1) 





0 20 40 60 80 100 
MEAN TRUNK OCCUPANCY IN PERCENT 


Fig. 3— 1 -> 0 transition probability for the 100-second sampling option. 
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——KILLER TRUNKS (r= 5) GUM AS TOsn MOREL 


Z = 1,2,5 ~' = 180 SECONDS 
T= 200 SECONDS 
Z= PEAKEDNESS 


a Z=2, 5 
(INDISTINGUISHABLE 
ON THIS SCALE) 


a 


NORMAL TRUNKS 





0 20 40 60 80 100 
MEAN TRUNK OCCUPANCY IN PERCENT 


Fig. 4— 1 > 0 transition probability for the GI/M/1-loss model with aswitched Poisson 
arrival process—2U0-second sampling option. 


a switched Poisson input stream (commonly used to model overflow 
traffic). Appendix A contains several details on the model and deriva- 
tion. It is clear that the effect of peaked traffic on the transition proba- 
bility is very small (z = 1 corresponds to a Poisson stream). 

Recent data? indicates that the service time distribution of a normal 
trunk has a coefficient of variation significantly greater than 1 (the ex- 
ponential case). Thus, in Appendix A we derive the covariance function 
of the server process x; for an M/G/1-loss! model with a mixed exponential 
type of service distribution. Figure 5 is a plot of P1,9 vs p parameterized 
by the coefficient of variation of the mixed exponential service distri- 
bution. We see that increasing the coefficient of variation has a noticeable 
effect on the transition probabilities, but the effect is to increase the 
discrimination between the normal and killer-trunk transition proba- 
bilities. 

Thus, it would appear that the transition probabilities based on the 
M/M/1-loss model are reasonably robust to perturbations in the trunk 
model. In addition, one suspects that using these transition probabilities 
in a detection scheme, which exploits the basic differences between killer 
and normal trunk transitions, might lead to a conservative design. 


t M/G/1-loss denotes a single server loss system with a Poisson (M) input stream and 
a general (G) service time distribution. 
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20 40 60 80 100 
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Fig. 5— 1— Otransition probability for the M/G/1-loss model with a mixed exponential 
service time distribution—200-second sampling option. 


3.2 Testing statistical hypotheses—a basic idea 


Suppose a trunk has constant mean occupancy p (Ex; = p) and we 
observe it for h seconds during which n switch counts and to 1 > 0 state 
transitions accumulate. We may write 


and 


where 


Hence, we have 


and 


E(ty0) = (m — 1)E([xp- — X(k-1)r]7) = (: = 1) pP 0. 


-‘m 
n= > Xr: 
pel 


m 
tio= Dd [xer = eA 


k=2 
h 0 z20 
m=-andz7 = : 
a 5 z<0 
h 
E(n) =p (~) 
: 


(5a) 


(5b) 
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Thus, for example, a trunk at 20-percent occupancy (sampled at the 
100-second rate) would accumulate 108 switch counts, on the average, 
in 15 hours. The mean number of 1 — 0 state transitions in this time 
interval corresponding to a normal trunk is 43, but the corresponding 
mean number of 1 — 0 state transitions for a killer trunk (during the 
same time interval) is 86. (Referring to Fig. 3, we see that P1,o(p = 0.2) 
is 0.40 and 0.80 for a normal and killer (r = 7.5) trunk, respectively.) 

With this example in mind, it is natural to consider a detection scheme 
of the following type: 


(1) Wait until we accumulate no switch counts on a trunk. 

(ii) When n = no, compare the accumulated 1 — 0 state transitions 
t49 to some threshold T'. 

(iit) Ifty9 > Ty, decide that the trunk is a killer, otherwise decide that 
the trunk is normal. [t¢19 is not directly available (see Section 2.1), but 
it may be estimated by t/2.| 


If the trunk occupancy were known and fixed, this scheme would ap- 
pear to be very reasonable. The analogy to the usual scheme suggested 
for deciding between a fair and a biased coin is clear: np is the number 
of (hopefully) independent experiments (analogous to the number of 
coin tosses), with each experiment having just two possible outcomes: 
the scan which follows the switch count is either 0 or 1. Thus, each switch 
count is associated with a 1 — 0 state transition (“heads”) ora 1— 1 state 
transition (‘‘tails”’). 

From the point of view of statistical hypothesis testing, we are thinking 
of two underlying states: 


Null hypothesis Ho: P;,o(p) = P trunk normal 


Alternate hypothesis H;: P1,9(o) = P* trunk killer 
Thus, our intuition suggests that a threshold test of the type sketched 
above is natural for distinguishing between Hp and H,. We will see 


(Section 4.1) that a (nonoptimal) test of this form arises naturally from 
pursuing the coin tossing analogy further. 


3.3 Problem formulation—sharpening the focus 
To simplify matters, assume to begin with that 


(i) The nominal mean holding time 1/y is known. 

(it) The trunk occupancy p is known and constant. 

(iti) The switch-count and transition data accumulations n(m) and 
t(m) are continuously available (scan by scan). 


With these assumptions, it is an easy matter to conceptually describe 
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the “optimum” scheme for deciding between the two simple hypothe- 
Ses, 


Ho: Trunk normal (mean holding time 1/p) 
H: Trunk killer (mean holding time 1/ry), 


with it understood that “trunk” refers to one of the specific models de- 
scribed in Section 38.1 (for concreteness assume the M/M/1-loss 
model). 

Let Xm = (%1,°-+%m) be the sequence of trunk states up to and in- 
cluding the mth scan, and let the available data be (as before) 


m m 
n(m) = > x;,t(m) = Y |x; — x;-4]. 
=) i=2 


Let 
Pim(n,t) = P(n(m) = n,t(m) = t/H;)i = O0orl 
and let 
P,,,(n,t) 
m(nt) =— 
ae SEE) 


The joint probability distributions P;,(n,t), 1 = 0,1, are well defined 
for any specific trunk model, but they may be nontrivial to derive. 2,,(.,.) 
viewed as a function of the vector random variable [n(m),t(m)] is re- 
ferred to as the likelihood ratio statistic and plays a central role in the 
theory of statistical hypothesis testing. More specifically, the optimum 
test (in a variety of senses) for deciding between two simple hypotheses 
involves suitably comparing ¢,, to a threshold (or thresholds) in order 
to make a decision. 

We briefly review two optimum tests, the Neyman-Pearson (fixed 
sample) test!° and Wald’s sequential probability ratio test (SPRT), using 
notation appropriate to our (discrete) problem. 


3.3.1 The Neyman-Pearson test 

Suppose a and @ denote the type 1 and 2 errors* of the test, 
Choose H, if €,, = T 
Choose Ho if €m < T, 


and suppose a’ and @’ denote the type 1 and 2 errors of any other test 
(requiring m samples) for deciding between Hp and H,. Neyman and 


* The type 1 and 2 errors, « and 8, are often referred to as the probability of false alarm 
and the probability of miss, respectively (a = probability of choosing H; given Ho is the 
true state, 8 = probability of choosing Ho, given H; is the true state.) 
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Pearson’s classical result is: if a’ < a, then §’ = 8. Thus, of all tests re- 
quiring m samples and having a false-alarm probability not exceeding 
a, the likelihood ratio test achieves the minimum probability of miss 
(maximum probability of detection). Since a = P(£ > T/Ho), choosing. 
asample size m and threshold T to achieve a S ao requires knowledge 
of the (conditional) distribution of @,,. Similarly, having chosen m and 
T, calculating B = P(£,, < T/H,) requires the distribution of £,,. (con- 
ditioned on H,). Note also that with such a fixed sample test, we decide 
in advance to accumulate exactly m samples before making a decision. 
In many contexts, data accumulates sequentially in time, and rigidly 
requiring m samples—independent of the particular realization that is 
unfolding—is not an optimal strategy. 


3.3.2 Wald’s sequential probability ratio test 


Using Wald’s SPRT, 1011 we continue to update ¢,, k = 1,2,--+ and 
defer a decision as long as ¢;€(T,T1). We make a decision the first time 
€, falls outside the interval (7'o,7';). Thus, 


if To <p, <T, R=1,2,-4m—-1 
and €m ¢ (T0,T)), 
then choose AH, if 2, = T1 
and choose Ho if €, = To. 


Clearly the stopping time m of the SPRT is a random variable, and the 
mean of m (given either hypothesis) is a measure of the time it takes to 
reach a decision. (Under a wide variety of circumstances, the SPRT ter- 
minates with probability 1.) Let E;(m) (i = 0 or 1) denote the mean 
stopping time, given that hypothesis 1 is in effect. Given a SPRT with type 
1 and 2 errors a and 8, and with mean stopping times E'9(m) and E;(m), 
consider any other test (sequential or not) with type 1 and 2 errors a’ 
and #’, and with mean stopping times Eo(m) and E,(m). The SPRT has 
the following optimal character! 


if a’ < aand &’ < B, 
then E,(m) = Eo(m) and E(m) = E,(m). 


Thus a SPRT is superior to a fixed sample test, if both tests have the same 
type 1 and 2 errors, in the sense that on the average it reaches a decision 
more quickly (under either hypotheses). 

In sharp distinction to the fixed sample test, the thresholds T'y and 
T required to approximately achieve specified type 1 and 2 errors are 
trivially determined.!! On the other hand, even determining the mean 
and variance of the stopping time is often a difficult chore. 
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In Section 4.2, we explicitly calculate the SPRT* for the simple hy- 
pothesis testing problem described at the beginning of this section. 
Before looking at this optimum test, however, we describe an ad hoc 
algorithm which is very robust and consequently attractive from a 
practical point of view. 


IV. INDIVIDUAL TRUNK ALGORITHMS 


A basic underlying assumption in this section is that the normal mean 
holding time of a trunk is known. Thus, if the algorithms in this section 
are designed relative to a normal mean holding time of 3 minutes, they 
will not discriminate between normal trunks having a mean holding time 
in the vicinity of 40 seconds,* and an actual killer trunk with the same 
mean holding time—both of these trunks will be detected as killer 
trunks. 

The rationale for studying this type of detection problem is two-fold: 
from the practical point of view the simplicity of implementation and 
general applicability! of these algorithms is attractive, and EADAS/ICUR 
can flag trunk groups which should not be studied by the killer trunk- 
detection algorithms (thus preventing false alarms on normal short- 
holding-time trunks). From the theoretical point of view, it was natural 
to consider this problem before factoring group information into the 
picture. 

Another modeling assumption used in this section (as well as in sub- 
sequent ones) is that the arrival process is stationary within data accu- 
mulation intervals, but the mean arrival rate may change arbitrarily from 
one accumulation period to another. Since we use equilibrium analysis 
(e.g., in calculating Po) we assume, in effect, that equilibrium is achieved | 
instantaneously. 


4.1 Anad hoc algorithm 


The essential idea of the test suggested in Section 3.2, is to decide on 
the state of a trunk by comparing the number of 1 — 0 state transitions 
(t19) to some threshold Ty, conditional on having accumulated a fixed 
number of switch-counts. We heuristically* proceed to derive such a test, 
using a standard likelihood ratio formulation, and explicitly take into 
account the time-variability of traffic. 

Let Xm = (%1,°++,Xm) correspond to the (unobservable) binary se- 


* Based on the M/M/1-loss model for a trunk. 

t Trunks in special-purpose trunk groups (credit checking, weather, etc.) will typically 
have mean holding times in the vicinity of 40 seconds. 

! The individual trunk algorithms can be used to test any trunk—regardless of the type 
of switching machine the trunk is associated with. 

* The distributional assumptions made in this section are intuitively motivated, but 
cannot be rigorously justified. We examine these assumptions carefully in Section 4.3. 
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quence of trunk states during an accumulation period in which m scans 
occur. Let t19(m) and n(m) be the number of 1 — 0 state transitions and 
switch counts associated with x,,. Denote the conditional probability 


P(ti9(m) = t/n(m) = n) 


for a normal and killer trunk by P(t/n) and P*(t/n), respectively. These 
conditional distributions depend, of course, on the trunk’s occupancy 
and on the particular trunk model we have in mind. [P*(t/n) also de- 
pends on the killer parameter r.] However, for the purposes of the heu- 
ristic development of this section, we do not precisely define which trunk 
model we have in mind. 

Since each switch count is associated with eitheral—-Ooral—>1 
state transition with probabilities P;9 and P;,; = 1 — P), respectively 
for a normal trunk, and since we expect successive transition events on 
a trunk to be essentially independent,* it seems reasonable to assume 
that P[t19(m) = t/n(m) = n] for anormal trunk is binomially distributed 
with parameters n and P). This same argument applies to a killer trunk. 
Denote the binomial distribution with parameters n and p by b(k;n,p) 
k = 0,-++,n, where 


b(k;n,p) = ) pri=p)er* 
Thus, we may think of a trunk with occupancy p during an accumulation 
period as having a conditional distribution 
P[tio(m) = t/n(m) = n] = b[t3n,P10(0,r)], (6) 


with r = 1 andr = ro corresponding to the normal and killer states of the 
trunk. (Recall that P;o(p,r) is essentially independent of r for r = ro with 
ro = 7.5 and 5.0 for 100- and 200-second sampling, respectively.) With 
these assumptions, we may think of testing the two simple hypothe- 
ses: 


Ho: P(t/n) = b(t3n,P1,0) Pio = Prole,1) 
Hy: P(t/n) = b(t3n,Pio) Pio = P1,0(p,ro)- 


If the 1 — 0 transition and switch-count accumulations for two successive 
and contiguous accumulation periods are (t1,n1) and (to,n2) respectively, 
we assume that 


P(t1,t2/ny,n2) = P(ty/n1)P(t2/n2). 


The idea here is that the only dependence between the two successive 


* The idea is that if significant correlation extends only one or two scans back, then 
successive transition events (events “triggered” by switch counts) should be essentially 
independent. 


TRUNK-DETECTION ALGORITHMS 933 


bit streams x}, = (x1,+++Xm) and x, = (xm44,°++X2m) is essentially due 
to the dependence between x,, and X41. 

Thus, if we denote the transition and switch-court accumulations for 
the ith accumulation period (in which m; scans occur) by [¢i9(m;),n(m,)] 
during which the trunk has occupancy p;, we have 


P{tyo(m4) = t1,*++tio(mp) = te|n (my) = ny,-++n(me) = ni] 


oa 


= [ bltisni-Pioleir)], (7) 


l 


where p; i = 1,---,k are the occupancies for the k accumulation periods. 

If tp = (t1,--+,t,) and mn, = (n4,-++,n;) consider the likelihood ratio: 

k bltisni,P1 (01,70) 

€(t,/n,) = . : 
MU Ble P vont) 

Denote the log likelihood ratiot log ¢(t,/n,) by €(t,/n,) and note 

that 





(8) 


A(tn/ me) = & blei/na, 


where 


b[ti;ni,P1,0(pi,ro)| 
b[ti3ni,P1,0(p:,1)] 
The expression ?(t;/n;) can be written as 


E(t;/n;) = a(p;)ti — a(p;)ni 


&(t;/n;) = log 








with 
_ 1 — Pi olp,1) 
a(p) = log (= ber) (9a) 
and 
_ P10(p,ro) 
a(p) = a(p) + log Pier (9b) 
Thus, we have 
. k 
e(t,/m,) = xu [a(p;)t; — a(pi)ni). (9c) 


Unfortunately, the occupancy in the ith accumulation period (p;) is 
unknown and hence equation (9c) cannot be used as a test statistic. One 
obvious “fix” is to estimate p; by p; = n;/m;, where n; and m; are the 


t1>T iff g(¢) > g(T) if g is monotone increasing, so the tests é > T and g(¢) > g(T)are 
equivalent. 
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switch count and the number of scans, respectively, during the ith ac- 
cumulation period. In stationary traffic, the estimate 


i i 
pi = DN; / uM; 
J=1 j=l 
would be used [o(A;) = 1/Vi o((;) if mj = m for all J]. 
Corresponding to the sequence of accumulations, (¢1,;,m;) i = 1,2,--- 
we define r; and R;, 1 = 1,2,--+ by 
ri = ac(pj)ti — a(pi)ni, pi = nj/m; (10a) 
and 
R; = R;-1 +r; with Ro = 0. (10b) 
Thus, we arrive at the sequential test: 


(1) Compute R;, 1 = 1,2,--- and defer making a decision as long as T'9 
<R; <T}. 

(ii) If i = k corresponds to the first accumulation period for which 
R; ¢ (To,T;), then 


Ry, = To = Trunk normal 
R, = T; = Trunk killer. 


If we ignore the fact that we are estimating p; by #;, and by assuming that 
the various assumptions made are valid (see Section 4.3), we identify the 
above test as Wald’s SPRT and as such 79 and 7; can be calculated as 
follows:!! to approximately achieve type 1 and 2 errors, a and 8, re- 
spectively, a + B < 1, choose 


To = log (=) (10c) 


and 





T, = log (- — =); (10d) 


a 


Throughout this section, we have assumed that the 1 — 0 transitions 
(tio) are available when, in fact, only the total transitions (t) are available. 
It should be clear that t19 can differ from t/2 by at most +44. To be pre- 
cise, let t19(™m), t9,(m) be the number of 1 — 0 and 0 — 1 state transitions 
corresponding to a bit stream x,, = (%1,++-+, Xm). If n(m) is the switch 
count corresponding to x,,, then we have 

n(m) = tyo(m) + tiy(m) + xm (11a) 
and 


n(m) = toi(m) + t11(m) + x1, (11b) 
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where t;(m) is the number of 1 — 1 state transitions. Therefore 
tyo(m) + Xm = toi(m) + x4, 


which together with t(m) = to,;(m) + ti9(m) yields 


(12a) 


tro(m) = 5 ttn) a —— 


ys 


and 
1 = 
toi(m) = — t(m) — Goa, (12b) 
2 2 
Thus, we can write the statistic update (eq. 10a) as 
t — 
a-—an+t+ farsa a 
2 2 


[It is easy to show that E[(x; — xm)a(p)] = 0.] 
We conclude this section with an interpretation of the statistic update. 
Rewriting the statistic update as 


r= (a—a)tyj9 — a(n — t49) 
and using eq. (11a), we obtain 
r= (a—a)ty9 — ati, — aX. (13) 


Now, from eqs. (9a) and (9b), it is clear that a > a > 0. Thus, each 1 > 
0 transition is weighted positively (evidence of a killer) while each 1 > 
1 transition is weighted negatively (evidence of a normal trunk). This 
is an intuitive explanation of the fact that the random walk (eq. 


(10b)) 
Ry = Re-1 + re 


has a positive drift if the trunk is a killer and a negative drift if the trunk 
is normal. 

The fact that the update assigns a negative weight (—a) whenever the 
last bit (x) is 1 uncovers a modeling deficiency. Recall that in eq. (6) 
we assumed 


P(tyo(m) = t/n(m) = n) = b(t3n,P}o), 


even though x,, = 1 can not contribute to an observable 1 — 0 transition. 
In this way we effectively modeled in a bias towards making “trunk 
normal” decisions. We can easily correct eq. (6) by conditioning on 
whether x,, = 0 or 1, obtaining: 


P(tyo9(m) = t/n(m) = n) = (1 — p)b(t3n,Pio) + pb(t3n — 1,P0). 


Now, proceeding as before in formulating the log likelihood ratio yields 
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a statistic R,, where 
R, = Ry-1 + fr 
and 
Fy =1r +t Qe, 


where r, is defined by eq. (10a), 


Pk tr l 
Vg) a) eee) 
AE = prJX\ net I= Prolbero)/ 


Pk tr i 
1+ (775,) (1-3) GP oth) 
1 — pp np/ \1 — Py o(6p,1) 


and pp, ty, and n» are the trunks occupancy estimate, 1 — 0 state tran- 
sitions, and switch count, respectively, during the kth accumulation 
period. 

Thus, we obtain our original test statistic with the correction term gq, 
added on. Note that gq; = 0, qx ~ 0 as py > 0, and qx > aas pp > 1, 
which is just the type of behavior expected, to offset the bias term in 
rp. 

Having heuristically developed an ad hoc sequential algorithm that 
is intuitively appealing and easily implementable, it is natural to ask: 
how does it compare to the optimum sequential algorithm? In the fol- 
lowing section, we rigorously develop an optimum sequential test. 


gr = log (14) 


4.2 An optimal algorithm 


Consider the two simple hypotheses: 
Ho: Trunk normal (mean occupancy p, mean holding time 1/z) 
Hy: Trunk killer (mean occupancy p, mean holding time 1 /roy). 
The optimum test for deciding between the two hypotheses—in the sense 
of minimizing the mean decision time—for given type 1 and 2 errors, is 
Wald’s SPRT (see Section 3.3), and it is based on the likelihood ratio 
statistic ,(t,n) given by 
ok = = 
ee P*(t(m) = t,n(m) = n) (15) 
P(t(m) = t,n(m) =n) 

Thus, it is clear that the ad hoc test described in Section 4.1 is not opti- 
mal, based as it is on an assumed conditional distribution, 


P(tyo(m) = t/n(m) = n). 


Before proceeding to study eq. (15), we must define the trunk model 
precisely. In the developments that follow, we model a trunk as the server 
in an M/M/1-loss system (see Section 3.1). The model implies that the 
sequence of trunk states x;, t = kr, k = 1,2,---is Markovian. Note that 
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although this appears to be a reasonable model for a normal trunk with 
200-second sampling, it ignores the conditional dependence “‘2 samples 
back,” which is more important for 100-second sampling—e.g., x; given 
x,-; is independent of x;—2, for the M/M/1-loss model. Taking this de- 
pendence into account in a trunk model would not be useful however, 
since the data needed to implement dependence “two scans back,” is not 
available. 

Since we are modeling the sequence of trunk states as a binary valued 
Markov process xx,, k = 1,2,-+ -,in equilibrium,it is clear that this process 
is characterized by @ = (Pi,0,Po,1), where P),9 and Po, are the transition 
probabilities 


P(xe47 = O/x; = 1) and P(xi4, = 1/x¢ = 0), 


respectively. (In general, a binary valued Markov process x;,,, k = 1,2,-++ 
in equilibrium, can be characterized by any two of the three quantities 
p, P1,0,Po,1. For our special Markov process (based on the M/M/1-loss 
model), both P;,9 and Po, and hence the process itself is determined by 
p alone.) Now having observed any m-tuple of the samples, which we 
denote by xm = (x1,°++,Xm), it is trivial to show that the statistic 


T(Xm) = (t(m),n(m),x1,xm) 


is a sufficient statistic for @. Thus, except for the initial and terminal 
states (x; and x), the transition and switch-count accumulations 
summarize all the “relevant information” in X;. 

Our hypothesis-testing problem can now be formulated as follows: 
X1,X%9,* ++ is a binary-valued Markov chain in equilibrium with parameter 
6 = (Po,1,P1,0) or 6* = (Po1,P 3,0). That is, our two states are 


Ho: {x;} Markovian, characterized by 6 = (Po,1,P1,0) 
Hy: {x;} Markovian, characterized by 6* = (Po1,Pi0). 


Now, because (t(m),n(m),x1,xm) is a sufficient statistic for 6, we know 
that the likelihood-test statistic based on the raw (unobservable) data 
Xm = (X1,*++,Xm) will be expressible in terms of t(m), n(m), x; and x», 
only. Thus, instead of studying eq. (15), we proceed (for simplicity) to 
study the likelihood-ratio statistic: 


P*(Xm) 


Pa. ue 


?(Xm) = log 
In Appendix B we study bmn (t,n) and find that it differs from ?(x;) only 
in an end-effect term. In @(x,,) this term depends on x; and x,,, whereas 
in £,(t,n) the corresponding term is a function of t and n. 
Since 
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P(Xm) = P(x) I P(xj/xi-1), 


we may write 


P(Xm) = P(x) Pi? Pit! Pe Pim 











So 
P*(x1) Po 1 — Pio 
?(xm) = lo tp lee eal G ) 
(Xm) °8 P(x ) 10 BD. 11 10g LPs 
Py 1- Po 
a ea | ( ). 17 
01 OB oo log 1—Po (17) 


Note that a trunk with mean occupancy p is busy and idle with proba- 
bility p and 1 — p respectively, independently of the state it is in (normal 
or killer). Thus, 


log —— ey) = 0 and eq. (17) can be written 
P(x) 
&(xm) = [(a — a)tyo — aty1] + [(8 — b)tor — btool, (18) 
where the parameters 6 and # are defined by 
1 > Poi 
b = log (————— 19 

- ; 7 B) ( a) 
B=b + log (5), (19b) 

Po. 


and the parameters a and a are defined as in Section 4.1 (eqs. (9a) and 
(9b)). Po and P;9 correspond to Po,i(p,r) and P1,0(p,r) with r = ro. 

Before discussing the symmetric structure of the optimum statistic 
[eq. (18)], we examine the Po, characteristics for the M/M/1-loss model. 
Using eqs. (2) and (3), we can obtain Po; vs mean-trunk occupancy p for 
a normal (r = 1) and killer (r = ro) trunk. Figures 6 and 7 are plots for 
the 200- and 100-second sampling option, respectively, with a mean 
holding time of 180 seconds. It is clear from Fig. 6 that a0 — 1 transition 
is just marginally more likely to occur on a killer trunk than on a normal 
trunk with a 200-second sampling rate. Although, the difference in the 
0 — 1 transition probabilities between a normal and killer trunk in- 
creases substantially with the 100-second sampling rate, it is clear that 
these differences are still quite small—compared to the spread between 
the P; 9 and Pio plots (see Figs. 2 and 3). Note that eqs. (2) and (4) show 
that 








P 1 — exp — 
—O) ae lforrg>1 (19c) 
Po,1 —s 
1 — exp 
Lp 
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M/M/1 — LOSS MODEL 
p>! = 180 SECONDS 
T = 200 SECONDS 


KILLER TRUNK 


r25 


NORMAL TRUNK 
{r= 1) 





fe) 20 40 60 80 100 
MEAN TRUNK OCCUPANCY IN PERCENT 


Fig. 6— 0 — 1 transition probability for the 200-second sampling option. 


M/M/1 — LOSS MODEL 
p-' = 180 SECONDS 
T= 100 SECONDS 


KILLER TRUNK 


NORMAL TRUNK 
(r= 1) 





0 20 40 60 80 100 
MEAN TRUNK OCCUPANCY IN PERCENT 


Fig. 7— 0 — 1 transition probability for the 100-second sampling option. 
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and, hence, we have G > b > 0. Using egs. (4) and (19c), we see that 
Po3/Po,1 = P}o/P10 and, therefore, 


B-b=a-a (19d) 
Equation (18) shows that the optimum statistic is the sum of two 


symmetric statistics: 


(i) The statistic [(a@ — a)t1q — at,,], which is essentially the ad hoc 
statistic (see eq. (13) and related discussion). 

(ii) An additional statistic [(G — b)to1 — btoo], which weights 0 — 1 
transitions positively (evidence of a killer) and 0 — 0 transitions nega- 
tively (evidence of a normal trunk). 


Note that by interchanging the role of 0 and 1 in either of these two 
statistics, we obtain the other—b is obtained from a and £ is obtained 
from a by replacing P19 with Po). 

By using eq. (11a) and the analogous equation 


m —n(m) = too(m) + toi(m) + xi, (Xn = 1 — Xm) (20) 
in eq. (18), the optimum statistic can be written 


&(xm) = [atio(m) — an(m)] + [Bto1(m) — b(m — n)] teil(xm), (20) 


where the end-effect term e;(x,,) is given by 
€1(Xm) = AXm + OxF. 


To implement ?(x,,) with only t(m) and n(m) available, necessitates 
estimating both t19(m) and to;(m) by t(m)/2. That is, using eqs. (12a) 
and (12b) in eq. (21) yields. 


f(x.) = | « au an(m) Fe [aS —b[m — n(m)}} rary re aay 





Z 2 
(22a) 
where 
e(X1,X%m) = (a — B) ea) + €1(xm) 
= (¢ =) (xy +xm) + bt (22b) 
or 
a at+ Bp 
é(x») = re) t(m) — (a — b)n(m) — bm + e(x1,%m). (23) 


t Recall that a — 8 =a — b [eq. (19d)]. 
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In the development of the ad hoc algorithm, we assumed that the 
statistics corresponding to successive accumulation periods are inde- 
pendent. We conclude this section by examining the independence as- 
sumption and with some remarks on implementation. 

Thus, turning to multiple accumulation periods, suppose x‘, 1 = 
1,2,---,k are the (unobservable) bit streams for k successive (and con- 
tiguous) accumulation periods, where x!, = iG- m+" *%Xim). ASs- 
suming stationary traffic, and noting that {x;}24 is Markovian, we can 
write 


LP(Xin44i/Xim) 


1 = 
P(xi,+++,xh) I Plxi) xT PGi) (24) 
and, therefore, 
“a k P* (xim+1/Xim)/P* (xim+1) 
O(xh,0++Xm) = + > log |e me (25 
(not m) = Ye E(x) = Ot Bene Pleund bo 





where P(-) and P*(-) denote the distribution under Ho (trunk normal) 
and H, (trunk killer), respectively. But, as we have seen, 
PF (Xinta/Xim) | —ros | 
eS |S exp (26) 
P* (Xim+1) lip 
and hence xjm+1 and xjm are essentially independent for ro sufficiently 
large. Therefore, 


P* (Sinai) Nin) = PE aa) 


and, hence, eq. (25) may be written as 
a k k A 6 k-1 
E(Xms Xm) = Do e(Xm) — LL (XimjXim4+1), 
i=1 i=l 


where 


PGineilsim) 
POipavr) 


is recognized as the mutual information random variable, which plays 
a central role in information theory. !? It is well known!” that (under Ho) 
E\I(xim3Xim+1)} is non-negative, and hence to ignore the end-effect 
term 


Lig tine) =] (27) 


k-1 
DX [(XimjXim+1) 


1w=1 


by implementing the statistic 


(Xin) 


II iM > 


i 
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would tend to make a normal trunk look more like a killer trunk on the 
average. In Appendix C, however, we show that the mean end effect 
E\I (xim3Xim+1)} is negligible compared to the mean statistic update 
E{é(xi,}. 

The “optimal” sequential algorithm is implemented in the same 
manner as the ad hoc algorithm (see eqs. (10a) and (10b)) except that 
now, corresponding to the sequence of accumulations (¢j,n;,m;)i = 
1,2,-- +, we define r; by 


ri = (<<) t; -(a- b)n; — bm;. 
The term ¢; denotes the total number of transitions in the ith accumu- 
lation interval. As is the case for the ad hoc algorithm (which corresponds 
to b = B = 0), the weights are functions of the trunk occupancy estimate 
pi = n;/m,. In a practical nonstationary environment, no claims of op- 
timality are made or implied. The term “optimal” is applicable only in 
the context of the equilibrium (e.g., stationary) model with known trunk 
occupancy. 


4.3 The ad hoc algorithm reviewed 


The assumption that t19(m) conditioned on the switch count n(m) 
is binomially distributed, is the basic assumption in the development 
of the ad hoc statistic. Although this assumption is incorrect (as we will 
soon see), the ad hoc statistic is essentially (except for an end-effect term) 
one of two symmetric statistics whose sum is the optimum statistic. Our 
purpose in this section is to examine the binomial assumption and to 
explain the relationship found between the ad hoc and optimal statis- 
tics. 

Since the optimal statistic was developed for a trunk modeled as a 
server in an M/M/1-loss system, it is natural to examine the binomial 
assumption (eq. (6)): 


Plt1o(m) = t/n(m) = n] = bf[t3n,P1,0(0,r)], (28) 


where Pj o(p,r) is given by eq. (4) in this context. Consider a killer trunk 
with r sufficiently large and suppose Xm = (X1,* + +,Xm) is the bit stream 
for a killer trunk during some accumulation period. Then, for all practical 
purposes [see eq. (26)], the trunk states x; i = 1,2,---+,m are independent 
and identically distributed Bernoulli random variables: 


p ifx =1 
l-p ifx=0. 


Thus, it is clear that the switch-count distribution on a killer trunk is 
the binomial: 


P(x; =x) = 
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P(n(m) = n) = b(n;m,p). (29) 


Now suppose A,,, (t,n) denotes the number of binary m-tuples having 
exactly t 1-0 transitions and n ones. If Xm = (x1,+++,Xm) is a se- 
quence of trunk states for a killer trunk with n(m) = n, it is clear that 
each such sequence has probability 


P(Xm) = p™(1 — p)™—" 
and therefore 
P(tio(m) = t,n(m) = n) = Am(t,n)p"(1 — p)™-" (30) 


for a killer trunk (with r sufficiently large). Equations (29) and (30) show 
that for a killer trunk, 


Pinions = (31) 


where we have used the fact that 


satan (7°) (0) @ 


It is interesting to note that while our assumed distribution for a killer 
trunk (28) differs from the correct distribution (31)—note that (31) is 
independent of p—there are some interesting similarities. For example, 
the assumed distribution peaks in the vicinity of (n + 1)(1 — p) and has 
mean equal to n(1 — p)t whereas the true distribution peaks in the vi- 
cinity of (n + 1)[1 — (n/m)] and has mean equal to n[1 — n/m]. Note that 
for “typical” realizations (x1,-++,x), we have 


a 

and, hence, the two distributions have the same general location and 
scale. [In fact, expression (31) is a hypergeometric distribution, which 
converges to (28) asm — © if n = pm (Ref. 13).] Thus, although incor- 
rect, the binomial distribution approximates the true distribution of the 
killer trunk. 

The following result helps to put the relationship between the ad hoc 
and the optimal statistics in perspective. 


Lemma 2: If {x;} ts a binary state stationary Markov chain with tran- 
sition probabilities Po, and P;, and if Xm = (X1,°++,X%m), then we 


+ Py olp,r) ~1-pasr—>o. 
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have 


P(xm) = 6(t10;n,P1,0) X b(to,14;m—n,Po,1) X q, (33a) 
where 
= P(x1)(1 — Pi0)7*"(1 — Powe (33b) 
Oars 
Cio to1 
Proof: 


P(Xm) = P(x) UL P(xi/xj-1) 


= P(e) Pe Pye Par Poe 
= P! (1 — Py 9) Pf! (1 — Po.3)! X P(x). 


Using eqs. (11a) and (18) to express f1; and tog in terms of tj and ft, 
respectively, yields the result. 

Thus, given a binary state stationary Markov chain {x;}, it is clear from 
the above lemma that the log likelihood ratio 


rs P*(xXm) 
£(Xm) = log —-— 
Vt) Ee 
formulated for the two hypotheses 
Ho: {x;} Markovian, characterized by (Po,1,P1,0) 
H;: {x;} Markovian, characterized by (Po1,P 10) 
is the sum of three terms: 
b(t10:n,Pio) b(to1;3m — n,Po,1) Af (“) 
og 
b(t1037,P 1,0) b(to1;m — n,Po,1) 








O(Xm) = log 


The first term is the ad hoc statistic (atj9 — an), the second term is the 
additional statistic [Gto, — b(m — n)], and the third term is an end-effect 
term (ax + bx‘,). 


** P* 
log =) =AXm+tbxot+ jog = axm + bxG, 
q P(x1) 
since, P*(x,) = P(x1) = p. 
The ad hoc algorithm, although based on the approximate binomial 
distribution, is very attractive for a number of practical reasons: 


(i) For the 200-second sampling option, it is essentially optimum in 
a practical sense, since the Po, characteristics for a normal and killer 
trunk are not far enough apart to exploit (see Fig. 6). 

(it) When we exploit the grouping information for the 5XB trunk 
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group in Section V, it will become obvious that the additional part of the 
optimum statistic is quite sensitive to the trunk occupancy and, hence, 
to the trunk-selection procedure modeled. (The sensitivity of the ‘‘ad- 
ditional” statistic to trunk occupancy stems from the fact that Po, is 
“almost” proportional to p.) On the other hand, we will see that the ad 
hoc 5XB algorithm is relatively robust to minor perturbations in the 
trunk occupancy (and therefore to the trunk-selection procedure) and 
hence might be expected to perform well in a real 5XB environment. 


V. THE 5XB TRUNK-GROUP ALGORITHMS 


In addition to utilizing individual trunk switch-count and transition 
accumulations the 5XB group algorithms exploit the following: 


(1) The identity of all trunks common to a group. 
(it) The trunk-selection procedure. 


The resulting 5XB group algorithms typically are faster* than their 
individual trunk counterparts and are also less sensitive to the groups 
nominal holding time. 


5.1 The 5XB trunk-group model 


For the purposes of this paper, we model a 5XB trunk group (with all 
trunks normal) as an M/M/N-loss model with random selection of idle 
trunks.? The same assumptions apply if the group contains one or more 
killer trunks, but in this case we assume that killer trunks have a mean 
holding time equal to 1/r that of the normal mean holding time. In ad- 
dition to being convenient theoretically, this idealized model has also 
been very useful in developing the 5XB group algorithms presently 
implemented in ICAN. 

If all N trunks are normal, the random selection rule implies that all 
trunks have the same mean occupancy. In Ref. 2, the birth and death 
equations for the above model with a single killer trunk were solved in 
closed form, and in Ref. 14 this was generalized to an arbitrary number 
of killer trunks. These analytic results turn out to be quite useful, and 
in what follows we will need the following results derived in Ref. 2. 


Theorem 1: For the above 5XB trunk group model having a single killer 
trunk with parameter r and an offered load of a erlangs, the blocking 
probability B(N,a,r) and the mean occupancy p*(N,a,r) of the killer 
trunk are given by: 


* For given type 1 and 2 errors, the group algorithms typically have a considerably 
smaller mean decision time than their individual trunk algorithm counterparts. 
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NB(N,a) 
Nr — (r — 1)a{1l — B(N,a)] 


0 —————— (34b) 


1 qi — B(N — 1,a)] 
a 


(i) B(N,a,r) = (34a) 


where B(N,a) is the usual erlang B blocking associated with an M/ 
M/N-loss system with all trunks normal and an offered load of a 
erlangs. 


It is easy to see that the occupancy p, of each of the N — 1 normal 
trunks must satisfy the conservation equation: 


rp, + (N — 1)p, = a[1 — B], (34c) 

and the trunk-group occupancy ¢; is defined by: 
PM = pr + (N — 1)pr 
r N ., 
(For a 5XB trunk group having a killer trunk with parameter r: p,; and 
p, denote the mean occupancy for a killer and normal trunk and ¢, de- 
notes the mean group occupancy.) 

Although eqs. (34a) through (34d) define an implicit relationship 
between p,(N,a,r) and ¢;(N,a,r), it will be very useful to have a simple 


explicit relationship. If the blocking term in eq. (34b) is ignored and if 
we “associate” ¢,; with a/N, an approximation suggested is: 


‘ or 
Sa os 35a 
Pr r= (p= 1)é ( ) 
This approximation, although quite good for large N, is rendered ob- 
solete by the following exact result: 


(34d) 


Theorem 2: Consider a 5XB trunk-group model with all trunks normal 
and mean-group occupancy @. (We will let ¢ denote the (mean) group 
occupancy for a5XB trunk group with all trunks normal.) If one of the 
trunks ts replaced by a killer with parameter r, then 

pr = p(¢,r), 


where, 


p(o,r) = ———* (35b) 


(r—1)o° 
Of course 
_ all — B(N,a)| 
ag N 
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where B(-, -) is the usual erlang B blocking expression. 

This surprising result, which follows easily from eq. (34b), is proved 
in Appendix D. As a consequence of this theorem, “5XB group occu- 
pancy”’ will be used to denote the occupancy of a 5XB group with all 
trunks normal. 

The mean occupancy of the normal trunks in a 5XB trunk group 
model having a single killer trunk no longer is given exactly by ¢. But 
the following result, derived in Appendix D, shows that ¢ is a good ap- 
proximation. 


Theorem 3: Consider a 5XB trunk group model with N trunks having 
a single killer trunk with parameter r = 1. The mean occupancy (p,) 
of the N — 1 normal trunks satisfy 


N-1 
r—(r—-1)¢ 


ii (<5) (r—1)¢ 
re OC (36) 


s 
IIV 
> 

IIV 


where ¢ is the mean-group occupancy with all trunks normal. 


Theorems 2 and 3 are proved in Appendix D, where an exact expression 
for p, is also derived. These results are special cases of general results 
obtained for the random selection model.!4 


5.2 Exploiting the 5XB Grouping Information 


To simplify matters, we assume that a trunk in a 5XB group* with 
mean-group occupancy ¢ has mean occupancy p(¢,r) given by 


d 

p(¢,r) Gee 
where r = 1 corresponds to a normal trunk. Thus, if a group has no killer 
trunks, all normal trunks satisfy p = ¢ and eq. (37a) with rr = 1 yields the 
correct occupancy. If, however, the group has a killer trunk, then all 
normal trunks satisfy inequality (36) so eq. (37a) with r = 1 is an ap- 
proximation that increases in accuracy with the size of the group. Of 
course eq. (37a) is exact for a (single) killer trunk in a 5XB trunk 
group. 

It is clear from eq. (37a) that 


’ (37a) 


7 p(¢,1) 
BT aa Tag) aD 


i We use “5XB group” and our idealized model of a 5XB trunk group interchange- 
ably. 
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and hence p(¢,r) is typically much smaller than p(¢,1). (For ¢ = 0.50 and 
r = 10, p(¢,r) = 2/11 p(¢,1). Thus, it would appear that considering the 
1 — 0 transition probability as a function of the 5XB group occupancy 
would effectively “spread” the Po characteristics in Figs. 2 and 3 further 
apart. That is, for a given ¢, we propose comparing 


Ple(¢,1)] and Pj o[o(¢,r)] 


[rather than P;9(p) and P}o(), as in Section IV]. 
Denoting the composition P; o[o(¢,r)] by P1,0(¢,r), we have 


Py,o(¢,r) = (1 — p(¢,r)) (1 — exp [—)) (38) 
1 — p(¢,r) | 
which is plotted in Figs. 8 and 9 for the 200- and 100-second sampling 
options, respectively. The normal holding time used in these figures is 
180 seconds, and the killer-trunk characteristics are drawn for r = 5, 10, 
and 15. 

The increased “spread” between normal and killer P; 9 characteristics 
obtained in this way is simply a consequence of exploiting the distinctly 
different occupancies of a normal and killer trunk in a 5XB trunk group. 
Figure 10 is a three-dimensional sketch of the composition of P;,9 and 
p. Because all normal trunks in a 5XB group have the same mean occu- 
pancy, we see that a single P19 vs ¢ characteristic suffices to describe 
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Fig. 8— 1 — 0 transition probability for the 200-second sampling option. 
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Fig. 9— 1 — 0 transition probability for the 100-second sampling option. 






P, 9/P PLANE:—~— 


M/M/1 — LOSS 
MODEL 


< 74 —— ~p,PPLANE: 5XB GROUP 
MODEL 


p(d.r) 


Fig. 10—Sketch of the composition of Pi 9(e) with p(¢,r). 


all normal trunks. This fact allows us to translate the individual trunk 
algorithm’s development to this 5XB context with essentially only no- 
tational changes. 
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Fig. 11— 1 — O transition probability for the200-second sampling option (mean group 
holding time = 45 seconds). 


Figures 11 and 12 are plots of eq. (38) drawn for a normal group mean 
holding time of 45 seconds.* The normal trunk characteristic corre- 
sponding to 180 seconds is shown in dashed lines. We see that with 5XB 
grouping information factored into the picture, considerable discrimi- 
nation exists between both normal trunk characteristics as well as be- 
tween the normal trunk having a holding time of 45 seconds and the killer 
trunks. The discrimination that exists between the normal trunks per- 
mits us to make the 5XB group algorithm adaptive to the group mean 
holding time. [Although we will not pursue this topic, the basic idea is 
that 2jt104)/2jn(V) (sums are over all trunks in the group) is an estimate 
of p10 and can be used to decide which (of several) normal pj charac- 
teristics constitutes H.| 


5.3 The ad hoc and the optimal 5XB group algorithms 


We assume that the mean group holding time is known and consider 
formulating a hypothesis-testing problem similar to that in Section 4.1. 
Thus, we denote P(ti9(m) = t/n(m) = n) by P(t/n) and consider the two 
hypotheses: 


* For normal holding times in the vicinity of 45 seconds, a killer parameter r in the range 
3 to 5 probably is typical. An r of 10 or 15 in this context is unrealistic. 
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Fig. 12— 1 — 0 transition probability for the 100-second sampling option (mean grouping 
holding time = 45 seconds). 


Ho: P(t/n) = b[t;n,P1,0(¢,1)] 
Hy: P(t/n) = b[t3n,Pi0(¢,r)], r= ro. (39) 


There are two differences between this formulation and the one in 
Section 4.1: 


(1) The trunk occupancy p in Section 4.1 is replaced by the group 
occupancy @¢. 

(it) The alternate hypothesis H, is composite since P;,o(¢,r) for r = 
ro are distinct. 


The approach taken in dealing with (ii) is a natural one often 
adopted;!! since P; o9(¢,r) is monotone increasing in r {this follows from 
eq. (38) upon noting that r/[1 — p(¢,r)] = r+ ¢/(1 — ¢)}, then testing 
between Ho and the simple alternate hypothesis 


Hy: P(t/n) = b[t;n,P1,0(¢,ro)], 


say with type 1 and 2 errors a and 8, respectively, implies that if the true 
state of nature is H; with r = ro > ro the resulting type 2 error will not 
exceed 8. With this approach, we can simply translate the ad hoc algo- 
rithm results developed in Section 4.1 to this 5XB group context by 
making the appropriate changes in notation. 
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Thus, the ad hoc 5XB group algorithm can be described as follows: 
corresponding to the sequence of data accumulations (¢/,n}) i = 1,2,--- 
for the jth trunk in a group of N trunks, where t/ and ni are the 1 — 0 
transition and switch-count accumulations, respectively, during the ith 
accumulation period in which m; scans are made, define s/ and S/i = 
1,2,- . 3 J = 1,--+; N by 

s! = a(¢;)ti — a(g;)nt (40a) 

and 
Si = Si_, +s! with Sj, = 0, (40b) 
where @; is the group’s occupancy during the 1th accumulation period. 


The sequential test for the jth trunk in the group, j = 1,---,N is defined 
by 


(1) Compute Si i= 1,2,---and defer making a decision as long as To 
<S1< Ty. 
(it) If i = k corresponds to the first accumulation period for which 


S! € (ToT), 
then 
Si < To = trunk j is normal 
Si = T= trunk j is a killer. 
The weights a(¢) and a(¢) are defined by 


5 LS Pi 0(¢,1) 
BG) 108 1 — P10(¢,ro) oa 
and 
ad) = a(4) + log Loe") (41) 


Pi 0(¢, 1) 


where P; 0(¢,r) is defined by eq. (38). 

Just as in the individual trunk algorithm, the actual occupancy re- 
quired to choose the weights a and a is unknown and must be estimated. 
Thus, the group occupancy ¢; during the ith accumulation period is 
estimated by 


b= AY dln (42) 
pi = Na i 


Note that 4; is a “better” estimator than pj; = n;/m; (the estimator used 
in the individual trunk algorithm) in the following sense: given a 5XB 
group with all trunks normal and mean-group occupancy ¢ (in equilib- 
rium), we have 
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(i) E@) = p=¢ 

(ii) E(b) = @ 

(iit) var(@) < var(p). 

In addition to the better occupancy estimate available on a group basis, 
the fact that the group P 9 characteristics are “flatter and broader” then 
the individual trunk P; 9 characteristic implies that the group algorithm 
more faithfully tracks the required weights, than does the individual 
trunk algorithm. 

The ad hoc 5XB group algorithm has the same pleasant intuitive in- 
terpretation that the ad hoc individual trunk algorithm had (see eq. (13) 
and related discussion). It is also easy to show how the optimal individual 
trunk algorithm development of Section 4.2 carries over to the 5XB 
group context. 

Thus, consider the two states of a trunk to be described by: 


Hg: {x;} Markovian, characterized by @ = (Po,1,P1,0) 
Hy: {x;} Markovian, characterized by 6* = (Po1,P 10), 
where P19 = P1,0(¢,1), Pio = P1,0(¢,ro), and [see eq. (2)] 


p(o,r) 
1 — p(¢,r) 


with p(¢,r) defined by eq. (37a). The assumptions that lead to a con- 
sideration of these two statistical hypotheses as a model of the normal 
and killer states of a trunk can be found in Section 4.2. 

Proceeding as in Section 4.2 leads us to the optimum statistic ?(x;) 
for distinguishing between the two simple hypotheses under consider- 
ation: 


Po1(9¢,r) = P30(¢,r) (43) 


&(xm) = [(a — a)t1o — ati] + [Btoo — (B — b)toi] + ce (44) 
P(x) 


where the parameters a and a are defined by eqs. (41a) and (41b) and 
the parameters b and @ are defined by 


- 1 — Poi(¢,ro) 
b = log ———— 45 
ean a Po,1(¢,1) on 
and 
7 j- Po 1(¢,1) 
= 6b + log ———_.. 45b) 
aces eC ET 


As in eq. (39), the alternate hypothesis is really composite since 
[Po,1(¢,r), Pi,0(¢,r)] for r = ro are distinct. Since Po,1(¢,r) is monotone 
decreasing in r, the approach discussed earlier of treating H; as asimple 
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hypothesis with r = ro is followed. The parameters 6 and f are defined 
a bit differently than were the parameters b and @ for the optimal indi- 
vidual trunk algorithm (see eq. (20a) and (20b)) in order to obtain non- 
negative weights. Thus, in the individual trunk algorithm context, we 
had Po; > Po,1 (eq. (19c)), but in the present 5XB group context we have 
Po, > Po1. The reason for this “flip-flop” is easy to see: for a given group 
occupancy ¢, we are now contrasting the 0 — 1 transition probability 
for two trunks which differ, not only in their hang-up rates but also in 
their occupancies as well. Thus, the difference in the two occupancies 
dominates the effect that the hang-up rates alone have. Roughly 
speaking, the 0 — 1 transition probability of a trunk is approximately 
equal to its occupancy (conditioning on the last scan has little effect) 
and hence since p(¢,ro) « p(¢,1) it is clear that we should have Po,1(¢,r) 
< Po,1(¢,1). Figures 13 and 14 are plots of Po,1(¢,r) for the 200- and 
100-second sampling options, respectively. In both figures, the killer 
trunk characteristics have been plotted for r = 5, 10, and 15 and a normal 
mean holding time of 180 seconds is assumed. These figures are very 
insensitive to the assumed normal mean holding time, since they es- 
sentially reflect eq. (37a), which is independent of the mean group 
holding time. . 
The “additional” statistic which appears in eq. (44), 


btoo — (8 — b)toy 
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Fig. 13— 0 — 1 transition probability for the 200-second sampling option. 
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Fig. 14—0 — 1 transition probability for the 100-second sampling option. 


shows that 0 — 0 transitions are weighted positively (evidence of a killer 
trunk*) and 0 — 1 transitions are weighted negatively (evidence of a 
normal trunk). This additional statistic is strongly influenced by the 
occupancy of a trunk, and only slightly by its hang-up rate. 

Note also that the term log P*(x )/P(x 1) is nonzero in the 5XB context 
since 


1 
P*(x 1) _|r —(r-1)¢ 
P(x) r 
r—(r-—1)¢ 


if x,=1 
if x, = 0. 


VI. PERFORMANCE OF THE 5XB GROUP ALGORITHMS 

In common with all sequential detection algorithms, the time required 
by the killer-trunk detection algorithms to reach a decision (trunk nor- 
mal or killer) is a random variable. In this section, we obtain an ap- 
proximate formula for the mean time required by the 5XB group algo- 
rithms to reach a decision. This result is used to contrast the performance 


* Killer trunks in the 5XB group model have very low occupancy, and hence 0 — 0 
transitions are likely. 
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of the ad hoc and optimal algorithms as well as to point out the consid- 
erable effect that the sampling rate has on each. In addition, we also 
obtain an approximate expression for the false-alarm probability of the 
5XB group algorithms. The analysis for the individual trunk algorithms, 
although differing in several respects from the group algorithms, involves 
the same sort of considerations and is omitted. 

The analysis in this section assumes a server system in equilibrium 
and, therefore, the mean trunk-group occupancy ¢ is assumed constant. 
In addition, to simplify the analysis, we assume that ¢ is known; an ap- 
proximate analysis which does not require this assumption is sketched 
in Section 6.1. A consequence of this assumption is that the algorithm 
weights are treated as constants rather than random variables. This 
assumption is not unreasonable because for multiple-hour accumulation 
periods var(¢) is quite small (¢ is the switch-count estimate of ¢). (Var(#) 
has been derived for an M/M/N-loss system.!°) 


6.1 Mean statistic update 

Corresponding to a sequence of trunk states x1,X9,+ + +X in an accu- 
mulation period with m scans, define a sequence of transition updates 
22,°* 42m by 


b if (Sn=1%n) = (0,0) 
wae (B-6)if — (xn-1,%n) = (0,1) AG) 
(a—a) if (Xn-1,4n) = (1,0) 


—a if (tp 21,%0) = (1,1) 
The optimum 5XB statistic (eq. (44)) may therefore be written: 


P*(xm) 


P(Xm) 


In practice the end-effect term cannot be implemented and all the 
transitions in eq. (44) must be estimated in terms of t(m), n(m), and m. 
Thus, if we denote the implementable version of eq. (44) by Sm(¢,ro),? 
use eqs. (11), (12), and (20) in eq. (44), and drop all end-effect terms we 
obtain 


A(xm) = Yo 2; + log (47) 
1=2 


Sm(¢,ro) = [ io An a(¢)n(m) 


(48) 


+ | Bort Gls) sgl 


t ro is the value of the killer parameter used in defining the alternate hypothesis and 
hence the algorithm weights (see eqs. (41) and (45)). 
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Now using eqs. (11), (12), and (20) once again, we see that S,,(¢,ro) may 
be written 


Sm(b,ro) = 3 2; + (bx%, — axm) + (a + B) aoa! (49) 
i=2 2 


Since the 5XB group is in equilibrium, 29, - -,zm are identically dis- 
tributed and it is easily verified that their common mean is given by 


E(z) = (a P19 — a)p + (6 — Po,18)(1 — p). (50) 


We recall that p (eq. (37a)) as well as the transition probabilities (eqs. 
(38) and (43)) are functions of the group occupancy ¢ and the state r of 
the trunk. The weights a, a, 6 and # are only functions of ¢ for a specified 
choice of the parameter (ro), which characterizes the alternate hypoth- 
esis. Thus, the mean statistic update for the optimum 5XB statistic and 
its “implementable version” is given by: 


» _ o P(x) 
E{0(xm)} = (m — 1)E(z) +E log ee (51a) 
and 
E\Sm(,ro)} = (m — 1)E(z) + 6(1 — p) — ap. (51b) 
Note that it is easily shown that 
P¥(x)) Site: r= 6) . 
E }log P(c,) | = —log [ro — (ro — 1)¢] + i ET log ro 





which is negative for r = 1 and positive for r 2 ro. 

Although the increments (transition updates) defined in eq. (46) are 
identically distributed, they are not independent. In fact, since the state 
sequence {x;} has been modeled as a Markov chain (see Section 4.2), it 
is easy to see that the sequence {z;} defined by (46) is also a Markov chain. 
Relabeling the four natural states, 


b, —(8 — 6), (a — a), —a, 
of the chain {z;} by 
0, 1, 2; 3, 


respectively, it is easily seen that the one-step transition matrix 7 for 
this chain is given by 


1-Poi Poi 0 0 
» 0 0 Pio 1- Pro 
ae oe ee 0 
0 0 Pio 1—Pio 


958 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1977 





(a) 


OPTIMAL STATISTIC 
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Fig. 15—-State diagram for the ad hoc and optimal 5XB group statistics. (a) Optimal 
statistic. (b) Ad hoc statistic. 


where 7;,; = P(Zp = J/2n—1 = i). Of course, the stationary distribution 
P satisfying Pr = Pis 


P=[(1— p)(1 — Pox), 1 — p)Po1, pP1,0, p11 — P10). 


The state diagram corresponding to the Markov chain {z;} is shown in 
Fig. 15a. If 6 = 6 = 0, we obtain the ad hoc algorithm, for which the se- 
quence {z;} is a three-state Markov chain, with natural states 0, a — a 
and —a. The state diagram for this chain is shown in Fig. 15b. 

Figures 16 and 17 are plots of the mean statistic update vs group oc- 
cupancy for the implementable version of the 5XB group algorithms. 
(Log base 10 is used in this paper.) These figures are drawn for a killer 
trunk with parameter r = 10. Also shown is the corresponding plot for 
a normal trunk (r = 1). 

We close this section by indicating an approximate analysis of E(z) 
which doesn’t assume that ¢ is known. Thus, in practice, ¢ is unknown 
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Fig. 16—Mean statistic update for the 5XB group algorithms—200-second sampling 
option. 
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Fig. 17—-Mean statistic update for the 5XB group algorithms—100-second sampling 
option. 
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and is estimated by ¢. Hence eq. (46) should read: 
= B(¢o) if (xn-1xn) = (0,0) and $= do. 


Zn 


Conditioning on ¢, we have 
E(z) = E{E(z/)}, 
where 
E(z/$) = 6(6)P(xn—1 = 0,%n = 0/6) +++. 


Now by assuming that (x,—1,x,) and ¢ are independent,* we get eq. 
(50) with the constant weights a(¢),--- replaced by the mean values 
E{a(@)},- - -. The mean values can be approximated in either one of two 
ways: 


(i) E}a(¢)} = e[E(¢)] = a(¢), which obviously amounts to assuming 
¢ is known. 
var(d) d2at 


Gry Ela OE seO) Si aa lay 





which factors the available variance of ¢ into the picture. 


6.2 Mean time to detection 


The basic structure of all the detection algorithms in this paper are 
the same: a statistic s; is evaluated at the end of the :th accumulation 
period, i = 1,2,--- and a decision is made the first time that the sum s1 
+ s9+--- falls outside an interval (To,71). Presumably, the random walk 
type statistic S; has a negative drift under Ho (trunk normal) and a 
positive drift under H, (trunk killer). Wald’s SPRT always has the ap- 
propriate drift: if Hp and H, correspond to the probability distributions 
Po(w) and P,(w), respectively, and if @(w) is defined by 


Pi(w) 
Po(w)’ 


then E{?(w)} < 0 under Hp and E{(w)} > 0 under Hj. The proof is im- 
mediate by using the inequality! 





o(w) = log 


— dpi log pi < — D pi log qi, 
L L 


* For reasonable-size trunk groups, we expect very little dependence between the 
sampled state process of an individual trunk (x1,---) and the group process ¢. 
t This is a Taylor expansion to second order. 
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where {p;} and {q;} are distinct probability distributions. In Appendix 
E, we show that both the ad hoc (6 = B = 0) and the additional (a = a 
= 0) parts of E(z) (eq. 50) have the appropriate drift. This in turn shows 
that both the ad hoc and additional parts of the 5XB group statistic (eq. 
(51a)) have the appropriate drift. 

Suppose Yj, Yo,--- are 1.i.d. random variables with common mean yu 
and consider a random walk 


So= VY; n=1,2. 
i=1 


with absorbing barriers at 7) and 7}. If u is small compared to Tp and 
T,, then the mean stopping time (mean number of steps to absorption) 
E(n) is approximately given by 


PoTo + PiT, 
m , 
where Po and P, are the probabilities of absorption at To and 7}, re- 
spectively. This follows from Wald’s identity!! 


E(S,) = wE(n) 


if we approximate the mean value of the random walk at absorption by 
PoTo + P,T}. 

In our detection theory context, Po and P; correspond to 8 and 1 — 
G, respectively, if the trunk is a killer (8 = probability of miss) and 1 — 
a and a, respectively, if the trunk is normal (a = probability of false 
alarm). If we denote the mean number of accumulation periods needed 
to reach a decision under Hy and H, by E(Tyn) and E(T;), respectively, 
and assume (1) successive statistic updates are independent and (ii) the 
mean statistic update is small compared to T'9 and T), we obtain 


(1 = a)To a aT; 


E(TxN) =-——_—_ . 52 
(EN) = "Fs (bral pee 
and 
BT9>+ 1 —- B)T, 
A = 52b 
oS ral Cen 


where the mean statistic update E{S,, (¢,ro)} is evaluated for r = 1 in (52a) 
and for r = rg in (52b). 

The assumption that successive statistic updates are independent is 
not strictly true if the successive statistic updates are contiguous. 
However, one expects that the slight (end-effect) dependence will not 
give rise to very much error. 

Let T;,, be the time required to decide (incorrectly) that a killer trunk 
is normal. Similarly, let T;,, be the time required to decide (correctly) 
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that a killer trunk is a killer. Then we have 
E(T;,) = BE(T xo) + A — BE(T» ,), 
which suggests that 
E(T,,) = E(T;,) if B<«1. (53) 


The moments of the conditional stopping times T),, and T;,, can be ob- 
tained by using a well known technique of Wald’s,* and one finds that 
approximation (53) is reasonable if JT; and (—To) are sufficiently 
large. 

In our faulty-trunk detection context, a type 1 error (‘‘false alarm”’) 
may result in the misuse of craft resources (e.g., testing a perfectly good 
trunk). A type 2 error (“miss”) on the other hand, will result in an in- 
creased time to detection. Assuming type 1 and 2 errors of 10~® and 1072, 
respectively (realistic implementation values), implies approximate 
thresholds To = —2 and T, = 6 (formulae 4.5c and d, log base 10). With 
these parameter values, the mean detection time E(T;,) can be ap- 
proximated by 


T 
ENS (¢,ro)} 


For a normal trunk, the mean statistic update is comparable to T) and 
hence expression (52a) isn’t applicable, nor is it needed since the 5XB 
group algorithm reaches a decision on a normal trunk after one or two 
updates. 

Figures 18 and 19 are plots of the mean-detection time vs group oc- 
cupancy for the implementable versions of the 5XB group algorithms. 
[The dashed line portion of Figs. 18 and 19 indicates where approxi- 
mation (54) involves considerable error (e.g., the region in which 
E{Sm(@,ro)} is a significant fraction of T;).] It is apparent from these 
figures that 


E(Tp,) = (54) 


(1) The mean detection time for both the ad hoc and optimal algo- 
rithms is enhanced by using the 100-second sampling option. This en- 
hancement is far more pronounced for the ad hoc algorithm. 

(it) The optimal algorithm is “faster” than the ad hoc algorithm. This 
contrast is greater for the 200-second sampling option. 


6.3 False alarm probability of the 5XB group algorithms 


If B(a,8) and A(a,@) are the test thresholds that result in type 1 and 
2 errors a and £, respectively, for a SPRT, Wald showed!! that using 


* See eq. 158 and 159 in Appendix A.5.2 of Ref. 11. The method ignores (as usual) the 
“excess over the boundaries” and hence yields approximate results. 
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Fig. 18—Mean time to detection for the 5XB group algorithms—200-second sampling 
option. 


thresholds T7'9 and T; defined by* 
B 
To =1 ( 
o = 10g 1 


— a 
yield type 1 and 2 errors a’ and ’ which satisfy 
and a’ s——. 


ar 
Bs l-a ~ [-£s 
The proof of this result is trivial and depends only on the assumption 
that the SPRT terminates with probability 1. This assumption is satisfied 
by a wide class of SPRTs,!® including the case of interest to us, where the 
underlying distribution is “Markovian”. Thus, the probability of false 
alarm (a’) for a SPRT (an example of which is the “optimal” 5XB group 
algorithm) satisfies a’ < 1077). 

Because the ad hoc algorithm is not a SPRT (the true underlying dis- 
tribution is not binomial, see Section 4.3), we may wish to study the 
consequences of using the above thresholds appropriate for a SPRT. T’o 
do this, we can proceed in (at least) two distinct ways: 





) and TT, = log (—) (assuming a + 6 < 1) 








(i) consider the statistic S »,(¢,ro) to be the basic update in the algo- 
rithm. 


* This assumes that the log likelihood ratio is used in defining the SPRT. 
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KILLER TRUNK PARAMETER, r= 10 
SINGLE~HOUR ACCUMULATION, m = 36 
UPPER TEST THRESHOLD, T, = 6.0 
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Fig. 19—Mean time to detection for the 5XB group algorithms—100-second sampling 
option. 


(it) Consider the transition update Z;(x;-1,x;) to be the basic update 
in the algorithm. 


In either case, we study a process of the form S; = Sl-, x;,j = 1,2,--- 
with absorbing barriers T') < 0 < T;. The advantage of proceeding as in 
(t) above is that the increments x; may be assumed to be i.i.d. We briefly 
sketch this approach. 

Consider a random walk S;, 7 = 1,2,--+ with i.i.d. increments {x;}. If 
n isa stopping time associated with S;, j = 1,2,---, Wald’s fundamental 
identity!!!” is given by 


EfeSty(t)-"} = 1 (all ¢ satisfying | x(¢)| 2 0), (55) 
where 
x(t) = Efe*4 


is the moment generating function corresponding to the common dis- 
tribution of the increments. If Pr, and Pr, denote the probabilities of 
absorption at J and T), respectively, then rewriting (55) in a standard 
way yields 


Pr EleSn'y(t)-"/Sy = Tol + Pr EleS='x(t)-"/Sn = Ti} = 1. (56) 


Now if the x; take on both positive as well as negative values with nonzero 
probability and have a nonzero mean, then the equation 


x(t) = 1 
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has only one nonzero root ho with the property that ho and E(x) have 
opposite sign.* That is, E(x) > 0 = ho < Oand E(x) <0 =hpo > 0. As- 
suming that the excess of S,, over the boundaries is small, eq. (56) yields 
the standard approximation!!}!7 


ehoT1 — 1 
Pry = ehoT — ehoTo == ehoTo (57a) 
and 
1- hoTo 
Pr, = - (57b) 


Zz ehoT1 = ehoTo 


Note that the probability of false alarm (type 1 error) corresponds to Pr, 
if the random walk increment is that of a normal trunk [Sp (¢,ro9), r = 
1]. Similarly, the probability of miss (type 2 error) corresponds to Pr, 
if the random walk increment is that of a killer trunk [S,,(¢,ro), r = 
rol. 

To use these approximations, we must compute the moment gener- 
ating function 


X(u) = Efe le(y)t/2-a(4)n]u} (58) 


by using the joint distribution p(t,n) derived in Appendix B. [Note that 
our discussion applies equally well to the additional and optimal sta- 
tistics. In general, we need E{exp [S(¢,ro)u]}, where Sm (¢,ro) is given 
by eq. (48).] Choosing the test thresholds 7’ and T, for the ad hoc al- 
gorithm according to the Wald SPRT formulae [eqs. (10c) and (10d)] 
typically results in Py, < B and Pr, <a. 


Vil. SUMMARY 


A class of killer-trunk detection algorithms has been developed that 
use the individual trunk usage and transition accumulations available 
in EADAS/ICUR. Because this data is essentially a sufficient statistic for 
the Markov chain used to model the (unobservable) sampled data, one 
of the algorithms developed is Wald’s celebrated SPRT. 

The detection algorithms developed can be partitioned in two natural 
ways: 


(t) By sampling rate (100 or 200 seconds). 
(it) According to whether grouping information is used. 


The algorithms which do not use grouping information are applicable 
to all trunks (one way or two way) independent of the type of switching 
machine used. A version of one of these individual trunk algorithms is 
currently in use in ICAN, testing trunks on 1XB, XBT and step-by-step 


* See Appendix A.2.1 of Wald’s original treatise (Ref. 11). 


966 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1977 


switching machines. The algorithms that exploit grouping information 
detect killers more quickly but are tailored to a specific switching ma- 
chine. A “group” algorithm of this type is currently being used to test 
trunks associated with 5XB switching machines. 

In the course of this study several problems of independent interest 
were studied. These include: 


(i) The server covariance in a M/G/1-loss and GI/M/1-loss system. 

(it) The structure of the likelihood statistic that arises in testing 
simple hypotheses characterized by a binary valued Markov chain. 

(iit) The occupancy of nonidentical trunks in a random-selection 
(Markovian) loss system. 


The major conclusion in this study is, of course, that accumulated 
switch-count and state-transition data on individual trunks (based on 
sampling intervals on the order of a normal holding time) can be used 
to reliably detect abnormally short holding time trunks. Moreover the 
(near) optimal sequential detection algorithms using this accumulated 
data are easily exhibited, simple in structure, and intuitively appeal- 
ing. 
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APPENDIX A 
Sensitivity of the Transition Probabilities to Modeling Assumptions 


To get an idea about the sensitivity of the algorithms to some of the 
modeling assumptions, the transition probability P;9 = P(x:+, = O/xt 
= 1) was studied for the following two cases: 


(t) M/G/1-loss, where the service distribution function F(-) is the 
mixed exponential given by 


F(t) = 1 — dye~*!* — doe ~*2é t=0 


(ii) GI/M/1-loss, where the arrival process is the switched-Poisson 
process’ commonly used to model overflow traffic. 


Because the methods used to obtain Po for these two models differ, 
we discuss these models separately. 
A.1 The M/G/1-loss model 


An observer viewing the server in an M/G/1-loss system sees an al- 
ternating sequence of busy and idle intervals. The busy intervals are 
distributed according to some distribution F(-) and are independent. 
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The idle intervals are exponentially distributed with mean \~! (\ = mean 
arrival rate) and are independent. Thus, the sequence of alternating busy 
and idle intervals constitutes an alternating renewal process. In this 
context, the conditional probability 


P10 = P(xi4, = 0/x; = 1), (59) 
where 
mae 1 : server , busy at epoch t 
0 ifserver is idle at epoch t 
has already been studied,!® and we have the following result: 
Theorem 4: Consider an M/G/1-loss system in equilibrium with a ser- 


vice time distribution F(t) having Laplace transform f*(s). If P},o(s) 
denotes the Laplace transform of P;,o(r), then we have 


u[l — f*(s)] 
s{s + d[1 — f*(s)]} 
where X and ware the mean arrival and service rates, respectively. 
Proof: See Section 7.4 of Ref. 18 [f5(s) = \/s + and fy(s) = f*(s)]. 


For the mixed exponential service time distribution mentioned above 
let 


Pi o(s) = (60) 


x1 = 2ud 
and 
xo=2Qu(1-d), Osdsl. 
If T is distributed according to this two-parameter (u,d) family of dis- 
tributions, then 
(1) E(T) = pol. 
20 





(ii) var(T) = p-2 ( ), where 5 = 4d(1—d). 


tiene (- = a 


BD 4% 


Thus, the mean is fixed at u~! and the coefficient of variation satisfies 
c = 1, with equality occurring when 6 = 1, which corresponds to the 
M/M/1-loss system. 
Equation (60) is easily inverted for this family of mixed exponential 
distributions, obtaining the following result: 
ryt 6 rot 6 


Pippa =p) + ae t+ ee. * (61) 
ry(r1 — re) ro(r2 — 11) 


where 
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r= - ss es E Sata) (5) |" 


2 
oo) -[o-oases QT” 


p44), oe 
1—p 





and p = \/(\ +p) is the mean occupancy of the server. 
In Fig. 5, P1.0(p,7) [eq. (61)] is plotted vs p for several different values 
of c (c = 1, 1.5, and 2) assuming 


(1) w~! = 180 seconds 
(iz) + = 200 seconds. 


Also shown is a plot of Pj,o vs p for a killer trunk with r = 10 [u in (61) 
is replaced by ru]. The normal trunk P, 9 characteristic with ¢ = 1 cor- 
responds to the M/M/1-loss system. 


A.2 The GI/M/1-loss model 


The covariance function R(-) for the GI/M/1-loss model with a 
switched Poisson arrival process has the form: 


R(r) = c, e¥!7 + cye¥27 + cge7@tY)7, (62) 


where the coefficients c; and the exponents w; are messy expressions 
involving the three switch parameters w, , and A (Ref. 8) and the mean 
service rate u. The derivation of this covariance function is straightfor- 
ward but tedious and is therefore omitted. (For the switched Poisson 
arrival prccess, the Markovian state equations can be solved for 
P(x;,x1+,), where x; = state of server at epoch t.) Our purpose here is 
to explain how eq. (62) was used in generating Fig. 4. 

If p is the mean occupancy of the server (p = E(x;)) and a is the offered 
load in a GI/M/1-loss system, then it is easy to show that the peakedness® 
z satisfies 





z=a(—*). (63) 


For GI/M/1-loss system the call congestion is ¢() and is related to the 
time congestion (p) by a(1 — ¢(u)) = p. So using z(u) = 1/[1 — ¢(u)] - 
a yields the result.? (#(-) is the L.S. transform of the interarrival time 
distribution.) Therefore, specifying p and z uniquely determines a. 
Hence, with a and z known, we obtain the equivalent random parameters 
and use the three-moment match to obtain the switch parameters (see 
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Ref. 8). Using this procedure, we obtain R(r) vs p parameterized by z. 
Equation (1) then yields P1. 


APPENDIX B 
The Likelihood Statistic Based on the Observable Data 


For ease of derivation, the likelihood statistic derived in Section 4.2 
was based on the raw data Xm = (X1,*°*,Xm) rather than on the observable 
data [¢(m), n(m)]; t(m) and n(m) are defined in Section 2.1. We will now 
study @,,,(t,n) and verify that the two statistics differ only in their end- 
effect terms. We will also examine the end-effect term based on t(m) and 
n(m) and show how it “tracks” the end-effect term based on x; and 
os! 

We begin by expressing the probability of x, = (x1,-++,Xm) in terms 
of t(m), n(m), x1, and xm: 


Lemma 3: If {x;} ts a binary state stationary Markov chain with tran- 
sition probabilities Po; and Py,9 and if Xm = (X1,°++,Xm), then 


ENTE Seneca 
PHP — Posy™-"-@/2), (64) 


where 
Q(x1,Xm) = Piyrtm?2 Py fritem)/2 Parr? Piytenii2+1 (65) 


and 


if x,=1 


a ae if x,=0 


Proof: This result is obtained from lemma 2 using eqs. (12a) and 
(12b). 
For convenience, we introduce the following notation 


(i) f GC n,P o) =P Pi/2(1 — Pi0)"—“/®. If t is even then of course 
t 
b(5 P ) 
9 n,f°1,0 


t 
pisP ) = 
f(g a n 


iva) 
(wi) g (—im = n,Po,) = P PY20 — Poi)™-"-"/2, 


(iii) Py y(t,n) = P(t(m) = t,n(m) = n,x1 = x,4m = y). 
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(iv) S,(t,n) = number of binary m-tuples satisfying t(m) = t, n(m) 
=n,xy =X, and xm = y. 


Thus, we can write 


1 


Mr 


P(t(m) = t,n(m) =n) = 


x 


Py y(t,n) (66) 


0 y=0 


and 


t t 
Pry(tin) = (n,Pro) & (Gm — nsPo1) POe)Sxy(tnQ(x,9). (67) 
Equations (66) and (67) imply that 


P(t, n(m) =t,n(m) =n) 


=f (SmPio) ¢ (5m ~ Pon)| SY Plx)Szy(tn)Q(xy)| (68) 


x=0 y=0 


and therefore it is easily seen that: 


P*(t(m) = t,n(m) =n) 
P(t(m) = t,n(m) =n)’ 


=| of -an|+[65-bin-m) | + 2m, (69) 


én (t,n) = log 


where 


x=0 y=0 


XY P¥(x)Syy (tn) Q*(x,y) 
LC) = le (70) 


= PRIS. En)OGY) 
x=0 y=0 


Comparing eqs. (70) and (22), we see that ?,,(t,n) and @(x») differ 
only in their respective end-effect terms. The following result is perhaps 
a bit surprising: 


Lemma 4: &(t,n) = (xm) if t is odd. 
Proof: t odd = x1 ¥ Xm = e(0,1) = e(1,0) = a + b/2 (see eq. (22b)). todd 
= Soo(t,n) = S11(t,n) = 0so E(t,n) for t odd may be written: 
(1 — p)So,1Q* (0,1) + pSi1,0Q* (1,0) 
(1 — p)S0,1Q(0,1) + oS1,0Q(1,0) © 


Now using (4) E(t,n) for t odd can be manipulated into the following 
form 





E(t,n) = log 


(1 — p)So1Q* (0,1) + pS1,0Q* (1,0) 


E(t,n) =1 ) 
(t,n) = log (1 — p)So,1Q(0,1) + pS1,0Q(1,0) 
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but Po,1/P1,0 = Po,1/P 1.9 = p/1 — p which completes the proof. 
If t is even then x; = Xm, 


e(0,0) = b 
and 
e(1,1) =a. 


For even t, eq. (70) can be written as 


E(t,n) = log 





Lp p 
Sic Gat S 
G - a) oa ey PY eas 
Thus, for even t, E(t,n) is a complicated function* of t and n. Note 
however that E(t,n) < 0 and 


a if Soo(t,n) =0 
E(t,n) = 
(in) f if S14(t,n) =0 
It can be shown that 
Pa: 
Sii(én) _ 2 
So,o(t,n) t 
m-n- . 


if (t,n) is such that S;,; and Soo are nonzero. 


APPENDIX C 
The End Effect E\I( Xim3Xim+1)} 
The end effect 
PT Otinn inesi)| = H(i, a9) ze A xineea/ Sim); 
where 
A (Xim41/Xim) = (1 — p)H (xim41/0) + pH (Xim+1/1) 
and 
A(xim+1), A (xim+1/0) and H (xXim4+1/1) 


are the binary entropy function A(x) evaluated at p, Po; and P;,, re- 
spectively (A(x) = —x logx —(1—x) log (1—x),0 Sx <1). 

Table I exhibits E{(xim3Xim+1)} and E{?(x,)} as a function of the 
trunks occupancy p—for a normal trunk, with a single-hour accumula- 
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Table | — End effect as a function of occupancy 


(7 = 100 seconds) (7 = 200 seconds) 

p E\I(xim+1;Xim)} E\2(xm)| E\I (Xim+1;Xim)} E\2(xm)} 
0.10 0.04 —1.36 0.01 —0.22 
0.20 0.05 —1.60 0.01 —0.20 
0.30 0.04 —1.50 0.01 —0.15 
0.40 0.03 —1.21 0.01 —0.09 
0.50 0.02 —0.84 0.00 —-0.04 
0.60 0.01 —0.47 0.00 —0.01 


tion period (m = 36 and 18 for the 100- and 200-second sampling options, 
respectively). The M/M/1-loss model is used with 1/u = 180 seconds. It 
is clear that the mean end effect is negligible compared to the mean 
statistic update. 


APPENDIX D 
Occupancy Formulae for a Random-Selection Loss System 


Theorem: Consider an N server Markovian loss system with random 
selection of idle servers and 


(1) N — 1 servers with mean service rate wu. 
(it) 1 server with mean service rate rp (r > 0). 
(itt) Mean arrival rate X. 


Let p, and p,; denote the mean occupancy of the servers with mean 
service rates u and ru, respectively, and let B denote the blocking (call 
congestion). Also let ¢ and @~ denote the carried load per server in an 
N server and N — 1 server Markovian loss system, respectively, (as- 
suming all servers have rate ), given an offered lead a = \/u. That is, 
¢ =a[1 — B(N,a)|/N and ¢~ = a[1 — BIN — 1,a)|/(N — 1),where B¢-, -) 
is the Erlang blocking formula. Then 


p 


(1) p; = ae (71) 
(ii) B(N,a,r) = ae [B(-, +) is the Erlang (72) 
blocking formula]. 
cha —(r-1 — 
(it) 62 p= ee errmeee (73) 
; _ gr —(r — 1)¢7 
Proof, ©) Pr = (ears Jee ) $. (74) 
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(1) This part of the theorem follows immediately from eq. (34b) 


‘ 1 
pr = 


1+ i r{1 — BIN — 1,a)] 
a 
by using the well known Erlang B recursion® 
B(N,a) a 
————— = — BIN - 1, 
cea No 
and the expression for the mean group occupancy 
a[l — B(N,a)] _ ’ 
N 


(tt) This part of the theorem follows immediately from eq. (34a) by 
using the above expression for mean-group occupancy. 


(tit) The lower bound part of eq. (73) follows from eq. (34c) by noting 
that B(N,a,r) < B(N,a)i and using eq. (71) of this theorem. Thus, we 
obtain 


re, 
r—(r-l1)¢ 
which can be arranged to yield the lower bound in eq. (73). The upper 
bound in (73) is an immediate consequence of (74) since 6 = ¢. Thus, 
it remains to prove (74). 
(iv) We prove this part as follows: Using eqs. (71) and (72) in the 
conservation eq. (34c), 


rp; + (N — 1)p, = a(1 — B) (75) 


(N — 1)p, 2 a[1 — B(N,a)] = No 


yields 
_a(r-1)A-¢@)+N—- ro 


= 76 
or (N= Dlr - = De] “ 
Therefore, eq. (74) holds if and only if 
si _ alr — I — ¢) + (N=r)g 
[r-(r-1)¢-] X= (N — 1) : (77) 
The right-hand side of eq. (77) can be rewritten as: 
rhs = ——— {r(N - 1 —a) ta —- B[N — (a + 1)r +.a]} (78) 


N(N — 1) 


t This follows eq. (34a) by noting that a[1 — B(N,a)] <= N. 
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using 
» - gl = BW.) 
N 


Similarly, using ¢~ = a[1 — B(N — 1,a)|/(N — 1), the left-hand side of 
eq. (77) can be written as: 


lhs = {r — (r — 1)(a/N — 1)[1 — BIN — 1,a)]}(a/N)[1 — BWN,a)]. - (79) 


Now using the Erlang B recursion formula in eq. (79) and rearranging 
terms yields eq. (78). Thus, eq. (77) and consequently eq. (74), of the 
theorem is proved. 


APPENDIX E 
The Sign of the Mean Statistic Update 
The mean statistic update of the optimal 5XB group algorithm is given 
by eq. (51a): 
» ~  s P*(X4) 
Et e(xm)} = (m — 1){(aP1,0 — a)p + (6 — BPo1)(. — p)} + DEBE 


The mean update of the ad hoc and the additional statistics correspond 
tob= 8 = Oanda = a = 0, respectively. The following lemma is needed 
to study the sign of E{@(x)}: 


Lemma 5: If 0<q <p < 1 define 





1- 
a = log ( S) and a=at log”, 
I—p q 
then 


a 
o> <p: 
a 


Proof: Consider the log likelihood ratio statistic @(m) for this simple 
hypothesis testing context: the observed process is n(m) = 272, x;, where 
x; are iid. Bernoulli random variables with P(x; = 1) given by q and p 
under Ho and Hy, respectively. Therefore, 


e(m) = og ———— = an —am 


b(n;m,q) 


and hence the mean of @ (m) is m(aq ae a) under Hp and m(ap — a) under 
H,. But as we noted in Section 6.2, E'{} is negative under Ho and positive 
under H;, in a general discrete setting.* Thus, we must have 


aqgq-a<0 
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and 
ap—-a>Q0O, 


which completes the proof. 
An immediate consequence of this lemma is that 


a(¢,ro) 


P 
(1) Py,0(¢,1) <———~ alera 


< P1,0(¢,ro) 


and 


NCCE sere < Poald,)), 


where a and « are defined by eqs. (41a) and (41b) and 6 and B are defined 
by eqs. (45a) and (45b). Since r = rp implies that P; o(¢,r) = P1,0(¢,ro) 
and Po 1(¢,r) = Po.1(¢,ro), we see that 


sgn(a(¢,ro)P1,0(¢,r) — a($,ro)) = sgn(b(¢,r0) — Po,1(¢,r)B(,ro)) 
positive ifrZ=ro 
negative ifr=1. 
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Cancellation of Polarization Rotation in an Offset 
Paraboloid by a Polarization Grid 


By T. S. CHU 
(Manuscript received July 20, 1976) 


The polarization rotation properties of the field radiated from a 
polarization grid have been found to be similar to those in the aperture 
of an offset parabolic antenna. This observation suggests broadband 
cancellation of the polarization rotation in a large offset reflector by 
the opposite prerotation of the incident feed radiation via a polariza- 
tion grid. Reduction of cross polarization from —24 dB without the grid 
to —39 dB with the grid wires is predicted in a numerical example. A 
previously unexplained polarization rotation measured using an offset 
parabolic grid is shown to be in good agreement with calculation. 


|. INTRODUCTION 


Current interest in frequency reuse through orthogonal polarizations 
creates a strong incentive for improving polarization properties of an- 
tennas. Cross-polarized radiation from an offset reflector! is often re- 
garded as a blemish on an otherwise excellent antenna, which offers both 
low sidelobe level and good impedance matching.” Although the cross 
polarization can be minimized using a large effective f/D ratio, the cor- 
responding requirements of small offset angle and large feed aperture 
are not always convenient in applications. Recently a trimode feed horn® 
has been discussed as a means of reducing the cross-polarized radiation 
of offset reflectors; however, a multimode arrangement is inherently 
narrow in bandwidth. Our purpose here is to propose a broadband 
method for reducing the cross-polarized radiation of the offset reflector 
antenna. This scheme is based upon cancellation of polarization rotation 
due to reflector curvature! by opposite rotation generated by a polar- 
ization grid.4 

In Section II we examine the similarity between cross polarization in 
the aperture of an offset reflector and that of a polarization grid. The 
condition of cancelling the first-order cross-polarization terms is then 
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Fig. 1—Geometry of offset reflector. 


deduced. In Section III numerical calculations of practical examples are 
given. 

In the appendix the previously unexplained measured polarization 
rotation* of a parabolic grid® is calculated. This example is given here 
as experimental evidence of agreement with the predicted rotation of 
radiation from a polarization grid. 


il. CANCELLATION OF POLARIZATION ROTATION 


Let us first briefly review the salient properties of the cross-polarized 
field in the aperture of an offset paraboloid as shown in Fig. 1. For a 
balanced feed radiation, 


cos ¢’ sin ¢’ _: 
Ey=FO',0)| v= / ee (1) 
sin ¢’ cos ¢’ 


where (9, 6’, ¢’) are spherical coordinates with respect to the z’ axis and 
(6’, 6’) are the corresponding unit vectors. Since the reflected field from 
the paraboloid is 


E,; = —E,; + 2n (Ky -n), 
* The design was initially proposed by Comsat for the Comstar satellite. 
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where /7i is a unit vector normal to the reflector surface. The principal 
and cross-polarized field components in the reflector aperture can be 
written respectively:! 


hp _ F(0’, #’) 


M=E,-. [sin 0’ sin 09 cos ¢’ — sin? ¢’ 
Yp tp 
- (cos 09 + cos 6’) — cos? ¢’(1 + cos 99 cos 6’)} (2) 
y F(0’,’) .. . ; 
N= E, -”? = mn) [sin 6’ sin 89 sin ¢’ 
Xp tp 


— sin ¢’ cos ¢’(1 — cos 6’)(1 — cos 49)], (3) 


where t = 1 + cos & cos 4) — sin 0’ sin 99 cos ¢’, M2 + N? = F?/p2, and N 
vanishes when 6 = 0. The offset angle 9 is between the feed axis and the 
reflector axis. The sign combination in eqs. (2) and (3) indicates that the 
rotation of the polarization vector due to offset in a paraboloidal aperture 
has the same magnitude and is in the same direction as illustrated in Fig. 
1 for any orientation of the incident linear polarization. The projection 
of the intersection of a circular cone (with vertex at the focus) and the 
offset paraboloid onto the x,y, plane is a circular aperture with cen- 
ter 


a 2f sin 0 
“cos 49 + cos 4 


(4) 


and radius 
2f sin 0, 
cos 09 + cos 0.” 


(5) 


where @. is the half angle of the cone. Equations (4) and (5) will be used 
later to obtain the relations in eqs. (12) and (13). 

Radiation from transmitting and reflecting wire grids can be obtained 
by magnetic and electric current sheet models, respectively. The prin- 
cipal and cross-polarized components are? 


P = —C[I1 — cos? ¢’ (1 — cos 6’) — sin 0’ cos ¢’ tan 6] (6) 
X =+C[sin ¢’ cos ¢’(1 — cos 6’) + sin 0’ sin ¢’ tan 6], (7) 


where C is a proportionality constant, #’ and ¢’ are the spherical coor- 
dinates of the feed (z’) axis, and 6 is the angle between the conducting 
wires and the x’y’ plane as shown in Fig. 2; the expressions inside the 
brackets are identical to those of eqs. (9) and (10) in Ref. 4, provided one 
makes the following substitutions: ¢’ = ¢ — 90° and 6 = 90° — y. The 
changes of notation are made for the purpose of comparison with eqs. 
(2) and (3). The upper and lower signs in eq. (7) correspond to the 
transmitting and reflecting cases. The orientations for the transmitting 
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Fig. 2—Configuration for cancellation between polarization rotations of an offset par- 
aboloid and a polarization grid. 


and reflecting polarizations together with a given grid geometry are 
shown in Fig. 2 where the conducting wires are parallel to the plane of 
the figure. 

One notes that the leading terms, which contain first power of 6’, in 
eqs. (3) and (7) have the same sinusoidal dependence on @’ and ¢’. Fur- 
thermore, the sign combination in eqs. (6) and (7) indicates that the 
transmitting and reflecting orthogonal polarizations rotate in the same 
direction, opposite to the rotation in the aperture of an offset parabo- 
loid. 

Let us take the first-order approximation—i.e., cos 0’ ~ 1—in eqs. (2), 
(3), (6) and (7), but sin 6” is retained. Then the cross polarization in eq. 
(3) normalized with respect to eq. (2) cancels that in eq. (7) normalized 
with respect to eq. (6): 


6 
= Sy jiras 


9 
P M o. (8) 


Polarization rotation of the radiation from a wire grid, as predicted 
by eqs. (6) and (7), also explains data measured using a cylindrical re- 
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flector made of a curved wire grid. The report on this experiment® mis- 
interpreted the polarization rotation as a consequence of diffraction and 
the offset geometry. A comparison between the calculated and measured 
polarization rotations of this case is given in the Appendix. 


Ill. NUMERICAL EXAMPLES 


Since the cancellation of polarization rotation discussed in the pre- 
ceding section only eliminates the leading terms, it is of interest to de- 
termine the residual cross polarization. Assuming that the offset reflector 
is located in the far zone of the radiation from a wire grid, as shown in 
Fig. 2, the principal and cross-polarized components in the reflector 
aperture can be written 


P = F(#’)[1 — cos? $p (1 — cos 8,) + sin 0, cos dp tan «| (9) 

X = +F(6’) [sin dp cos dp (1 — cos 0p) — sin Op sin op tan e], (10) 
where F(0) is the feed-radiation pattern and 

0’ = cos! [cos 6, cos 99 + sin Op sin 99 cos dp). (11) 


The above equations are simply a decomposition of the grid radiation 
into the two orthogonal components of a balanced feed whose axis 
coincides with the paraboloidal axis. The expressions inside the brackets 
of eqs. (9) and (10) are of the same form as those of eqs. (6) and (7); but 
0, and ¢p are the spherical coordinates with respect to the paraboloidal 
(z,) axis instead of the feed axis, and « = (69 — 6) is the angle between 
the conducting wires and the xpy, plane, as shown in Fig. 2. 

To relate eqs. (9) and (10) to the normalized aperture coordinates r 
= (p,/a) and gq, the following expressions can be obtained with the aid 
of eqs. (4) and (5): 


Biogen E (r sin 0. sin ¢q)* + (sin 09 +r sin 8, cos ¢,)? | 
e cos 8. + cos 4 
(12) 
bp = tan=! | rsin be sin og i (13) 
sin 09 + rsin @, cos dq 


Numerical examples of several combinations of parameters (9, 0. and 
¢) have been calculated for the principal and residual cross-polarization 
components from eqs. (9) and (10). The feed pattern has a gaussian shape 
with 10-dB taper at the edge of the reflector. The principal polarization 
is close to unity (0 dB) around the center of the reflector aperture. The 
maximum, calculated, residual cross polarization is given in Table I for 
a number of examples. Fig. 3 shows a plot of both principal and cross 
polarizations for the case 4) = 50°, 0. = 20°, and « = 25°. Only half of the 
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Table | — Cross Polarization in the aperture of an offset reflector 




















Max. Residual Cross Max. Cross Pol.® 
Oo A. € Pol. With Grid Balanced Feed 
(Deg) (Deg) (Deg) (dB) Without Grid (dB) 
50 20 25 —38.6 —24.0 
50 20 23 ~—36.4 —24.0 
50 20 27 —36.1 —24.0 
60 20 30 —38.0 —22.5 
60 30 30 —30.9 —18.0 
90 20 45 —34.3 —17.5 
90 14 45 —40.5 —20.0 








aperture needs to be shown for each polarization because of the sym- 
metry. The maximum residual cross polarization in this case is —38.6 
dB, reduced from —24 dB for the same reflector aperture illuminated 
by a balanced feed without a polarizer. Keeping the same set of param- 
eters, Ip = 50° and 0, = 20°, the calculated residual cross polarization 


PRINCIPAL POLARIZATION (dB) «—— | ——®* CROSS POLARIZATION (dB) 














Fig. 3—Relative amplitude levels in the aperture of an offset paraboloid. #9 = 50°, 0 
= 20°, « = 25°. 


982 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1977 


becomes —36.4 dB for « = 23° and —36.1 dB for « = 27°. These results 
indicate that the residual cross polarization is not overly sensitive to a 
slight departure from the optimum orientation of ¢ = 69/2. 

The examples for 09 = 60° and « = 30° show residual cross polariza- 
tions of —38.0 dB and —30.9 dB for #, = 20° and 30°, respectively. The 
second-order terms are not quite negligible at 0. = 30°; however, the 
cancellation of cross polarization appears to be significant. 


IV. DISCUSSION 


In view of the residual second-order (1 — cos 6’) terms, the half-cone 
angle of the reflector subtended at the focus should not exceed about 
20° in order to take full advantage of the cancellation. If the reflectors, 
such as in an offset cassegrain configuration, do not cause significant 
cross polarization, the conducting direction of a polarization diplexing 
grid should be oriented to avoid introducing any polarization rotation, 
as discussed in Ref. 4. 

When 0p — @, is less than about 30° and ¢ = 6 = 6/2, the feed horn 
associated with the polarization reflected from the grid will introduce 
some blockage, as shown in Fig. 2. The blocking problem can be eased 
by using a smaller value of «. This practical difficulty may prevent op- 
timum orientation of the grid wires for reflectors of small offset angle, 
and hence reduce the effectiveness of the cancellation. The practical 
application of this scheme should indeed lie in reflectors with large offset 
angle. 

The expression € = 60/2 implies that the grid wires are approximately 

_parallel to the tangent plane at the center of the offset reflector. One 
notes the similarity between this case and a symmetrical small-cone- 
angle paraboloid illuminated by a grid-covered feed. 
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Appendix 
Polarization Rotation of an Offset Parabolic Grid 


An offset cylindrical-reflector system fed by two line sources of pillbox 
type (as shown in Fig. 4) was proposed? as a dual-polarized antenna with 
an elliptically shaped coverage pattern. The reflecting system consists 
of a vertically polarized grid attached to the surface of a parabolic cyl- 
inder, the front surface of the grid having the same curvature as the cy- 
lindrical reflector. The measured data® showed excellent orthogonality 
(within 1°) between vertical and horizontal polarizations over the whole 


PARABOLOID ANTENNA 983 


BEAM DIRECTION 








_- OFFSET 
-~~ PARABOLIC 
CYLINDER 


_ — OFFSET VERTICALLY 
POLARIZED GRID 


GR __—GRID 
ti LINE > “VERTEX 
Ae 
/ 
VERTICALLY — —— >; fk ~~ PARABOLIC 
POLARIZED / XN CYLINDER 
~PARABOLIC VERTEX 

LINE SOURCE Cv LINGER 

FOCAL LINE 


| 
HORIZONTALLY 
POLARIZED 
LINE SOURCE 


Fig. 4—Configuration of a dual-polarized cylindrical reflector antenna (Ref. 5). 


6.8° X 3.4° (3 dB) elliptical beam. However, significant polarization 
rotations, in the same direction for both polarizations, were observed. 
Maximum rotations of about 2° occur on the major axis (xz plane) of the 
half-power ellipse of the beam. No adequate explanation was given for 
this measured polarization rotation. Here we explain the rotation using 
electric and magnetic current sheet models for the grid. 

The first-order approximation of polarization rotation by a wire grid 
is simply sin @ cos ¢ cot y (eq. 7 in Section II or eq. (10) in Ref. 4), where 
6 is the angle off the beam axis, cos ¢ is unity in the xz plane of maximum 
rotation, and y the angle between the conducting wires and the beam 
axis. Since the wires of the curved grid have a variable direction, the 
rotation A in the plane of maximum rotation can be obtained by aver- 
aging cot y over the parabolic curve, y2 = 4fz. 


sin 0 ses oY (dy)? + (dz)? 
EVE 


_sin 6 


Sr f Ecoty ds = 


(14) 


where £ is the aperture field distribution. Let u = y/2f, eq. (14) be 
transformed into 
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Fig. 5—Measured angles of nominal vertical and horizontal polarizations around —3-dB 
pattern contours of 6.8° X 3.4° for a dual-polarized cylindrical antenna with a parabolic 
grid; + indicates locations of measurement. 


1 
sin 0 EV1+u2udu 
0.1763 
SS SS, (15) 


1 
J, EvV1+u2du 
0.1763 


where the upper and lower limits of integration are obtained from y = 
2f tan y/2 with wy = 90° and 20°, respectively. A is not sensitive to the 
aperture distribution. Numerical calculation gives A = 0.61 sin 6 for both 
cases when FE is assumed to be uniform and when E has a 10-dB edge 
taper—i.e., when 


y 27 
E = -0.396 + 4.746 ip 4.034 (5 ;) (16) 
The quadratic form has been chosen to perform the integration in closed 
form. When 0 = 3.4°, A = 0.0362 rad = 2.07°; agreement between this 
calculated value and Wilkinson’s measured rotation shown in Fig. 5 is 
indeed very good. 

Furthermore, eqs. (6) and (7) in Section II indicate that the trans- 
mitted and reflected orthogonally polarized fields from the grid rotate 
in the same direction, just as in the measured rotations. The planes of 
measured maximum and null rotation also agree with the predictions. 
Since the two polarizations have the same sense of rotation, orthogonality 
is preserved. However, if this rotation occurs in the feed radiation, 
nonorthogonal elliptically polarized radiations from the reflector will 
result. 
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Ten Years of Power Aging of the Same Group 
of Submarine Cable Semiconductor Devices 


By A. J. WAHL 


(Manuscript received January 18, 1977) 


The active devices in the first Bell System transistorized submarine 
cable system (SF) are unpassivated, diffused-base, germanium tran- 
sistors and oxide-passivated, silicon diodes. At the date of this writing 
these devices have accumulated over 550 million device hours of pow- 
ered operation in service and on the aging rack without failure or im- 
paired device performance for a demonstrated failure rate of less than 
0.00045 percent per thousand hours (4.5 FITs) with 90- percent confi- 
dence. This paper reports details of the behavior of 500 of these devices 
that have reached ten years of controlled power aging. The overall re- 
sults indicate that initial reliability objectives are being achieved and 
that the semiconductor devices should not be expected to limit the 
desired life span of SF submarine cables. 


I. INTRODUCTION 


In the initial stages of the project, the decision to use germanium 
transistors in the first Bell System transistorized submarine cable was 
a rather bold gamble to take advantage of their then superior frequency 
capability while also attempting to mitigate any reliability inferiority 
to silicon transistors. When the decision was made, in the time period 
of 1960-61, silicon transistors were clearly the wave of the future and, 
although still inferior to germanium transistors in frequency capability, 
had already indicated their reliability superiority. Nevertheless, the 
results reported here, together with the results of regular aging and actual 
service experience thus far, offer rather convincing evidence that the 
desired reliability objectives are being achieved. The first submarine 
cable using these transistors, a relatively short cable with 136 repeaters, 
has been in service over eight years. A transatlantic cable with 363 re- 
peaters has been in service nearly seven years. No semiconductor device 
failure nor degradation toward failure has yet been observed in any of 
the systems that have no redundancy and in which the failure of even 
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Fig. i—Collector-base breakdown voltage in volts at 100 vA for L2287 transistors— 
before and after 10 years of power aging at 75 mW. 


one device could disable the system. At the time of this writing the de- 
vice-hours of service in these cables totals more than 250 million for both 
transistors and diodes. On the aging rack these devices have accumulated 
over 300 million device-hours of power aging—also without failure. On 
the basis of these results the demonstrated failure rate, with 90 percent 
confidence, is less than 0.00045 percent per thousand hours (4.5 
FITs). 

The method of providing semiconductor devices for the SF submarine 
cable system has been described previously.! As a further part of this 
program, a hundred devices of each of the five basic device codes of 
transistors and pn diodes have been retained on extended power aging. 
A milestone was observed in this program when the aging time on these 
devices reached ten years. The primary purpose of this paper is to show 
in a condensed way the behavior of these devices across ten years of 
power aging. 

The transistors are of diffused-base germanium with no passivation 
or junction protection on the active chip. The diodes are of diffused 
silicon with grown oxide passivation over the junctions. All devices are 
encapsulated in hermetic packages. 
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Fig. 2—Collector-base reverse leakage current in nA at 15 V for L2287 transistors— 
before and after 10 years of power aging at 75 mW. 


ll. DATA COLLECTION 


All transistor parameters are measured before the devices are placed 
on long-term power aging, which is a period of about six months for 
candidate devices. From previous work with devices of this kind, we 
found that device current gain was the only parameter expected to show 
any consistent aging trend, and it is also the most critical parameter for 
the repeater amplifier circuit. For practical testing reasons during 
long-term aging, the de current gain hyg was measured. From the 
standpoint of the repeater circuit, however, the critical parameter is the 
low-frequency ac current gain hy, rather than the de current gain hrp. 
The two are not the same, but they bear a relationship to each other, as 
described in any standard text on transistor theory. If one shifts, the 
other shifts in the same direction. 
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Fig. 3—Emitter-base breakdown voltage in volts at 100 wA for L2287 transistors—before 
and after 10 years of power aging at 75 mW. 


Previous work with transistors of this kind had also indicated that any 
interruption of the bias on the device would cause the current gain to 
suffer a disturbance that would tend to obscure the real aging trend. For 
this reason, during long-term power aging, only the dc current gain was 
measured, and the measurement was made very precisely without dis- 
turbing the bias on the device. From the standpoint of device parameters 
directly relevant to circuit performance, therefore, the complete picture 
of the aged transistor cannot be known until the device is removed from 
long-term aging and all device parameters have again been measured. 
For such transistors, the aging time interval is ten years rather than six 
months. 

Since the aging trends of diode parameters were not disturbed by 
change of electrical bias for testing purposes, most of the sensitive diode 
parameters were tested during aging. The diodes, therefore, were not 
removed from the aging rack at the end of the ten-year aging time. In- 
terim aging information on these diodes has already been published.? 
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Fig. 4—Emitter-base forward voltage in volts at 15 mA for L2287 transistors—before 
and after 10 years of power aging at 75 mW. 


lll. PRESENTATION OF THE AGING RESULTS 


A rather simple and clear picture of device aging behavior can be 
shown in the scatter plot in which the initial value of a device parameter 
is plotted vs the final value across the aging period for each device. For 
the transistors in this case, the information obtained at times between 
the initial and final values is not available for reasons discussed pre- 
viously. For the diodes the information is available, but not significant 
since general aging trends, where they exist at all, are monotonic. For 
simplicity and consistency, therefore, all the ten-year aging results are 
shown here in the form of scatter plots. 


IV. COMMENTS ON THE RESULTS 
4.1 General 
In most cases where any detectable parameter shifts have occurred 


across ten years of aging, the magnitudes of the shifts have been small 
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Fig. 5—Low-frequency common-base current gain for L2287 transistors—before and 
after 10 years of power aging at 75 mW. 


and very uniform among all the devices, almost as if the devices had 
shifted in unison. Such behavior immediately raises the suspicion that 
the apparent shift may actually be an offset in testing rather than a real 
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Fig. 6—Typical behavior of dc common-base current gain for diffused-base germanium 
transistors on long-term aging. 
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Fig. 7—Noise figure in dB at 5.5 MHz for L2287 transistors—before and after power 
aging at 75 mW. 


shift. Unless the testing offset can be definitely confirmed, however, the 
shift must be accepted at face value. In the latter case, the uniformity 
of the shift makes the extrapolation very simple for more devices and 
longer time. The extrapolated result is then quite obvious and com- 
fortable. Possibly the only deviation from this situation is the behavior 
of the leakage current of the L2320 diodes. Further discussion of this case 
follows under the heading of that particular diode. 


4.2 L2287 transistors 


L2287 transistors are split into 38A, 38B, and 38D codes for use in the 
first two stages of the amplifier and in the oscillator of the SF repeater. 
The two amplifier stages operate at power levels of 25 and 75 mW, re- 
spectively. The oscillator operates at a lower power level. During aging, 
all of these devices are held at a power level of 75 mW, composed of 5 V 
and 15 mA. The extended aging program has included one hundred of 
these transistors. 

As shown in Fig. 1, the collector-base breakdown voltage of these 
transistors showed almost no discernible change after ten years of power 
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Fig. 8—Collector-base breakdown voltage in volts at 100 uA for L2288 transistors— 
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before and after 10 years of power aging at 0.99 W. 
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Fig. 9—Collector-base reverse-leakage current in nA at 15 V for L2288 transistors— 
before and after power aging at 0.99 W. 
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Fig. 10—Emitter-base breakdown voltage in volts at 1 mA for L2288 transistors—before 
and after 10 years of power aging at 0.99 W. 


aging, with the exception of one device which showed a decrease of about 
1.5 V. This device showed an increase in collector-base leakage current 
of about 0.15 wA, as seen in the lower part of Fig. 2, and a decrease in the 
emitter base forward voltage of about 5 mV, as seen in the lower part of 
Fig. 4. These changes, however, would have no effect whatever on the 
operational reliability of the device. No attempt has been made to de- 
termine the cause of the changes since all parameters of the device are 
still well within initial specification limits. 

Figure 2 shows that the lower values of collector-base reverse-leakage 
current show no discernible change while those above about 0.5 nA show 
a very slight upward shift. Such a small shift, however, is of no practical 
significance. 

The emitter-base breakdown voltages of most of the devices show 
remarkable stability, as depicted in Fig. 3. Only those above about 
1.6 V show a small amount of scatter in the behavior. With the one ex- 
ception, previously noted, the emitter-base forward voltages showed 
essentially no discernible change after ten years of power aging, as seen 
in Fig. 4. 

Perhaps the most gratifying result of this extended aging program is 
the behavior of the low-frequency common-base current gain hyp, as 
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Fig. 11—Emitter-base forward voltage in volts at 90 mA for L-2288 transistors—before 
and after ten years of power aging at 0.99 W. 


shown in Fig. 5. From the beginning, this parameter was expected to 
show a significant aging trend. Even under the conservative assumption 
that the aging of this parameter is linear with time, Fig. 5 indicates that 
in ten years of aging the average of the parameter has shifted only about 
a quarter of the way to end-of-life for the group. Actually, the aging of 
this parameter is known to occur at a decreasing rate with increasing 
time, as shown in Fig. 6 for a typical case of the aging of the dc com- 
mon-base current gain. With similar behavior of the ac common-base 
current gain, the end-of-life for the group of devices would appear to be 
at some indeterminate time far beyond forty years. 

Only one-third of the L2287 transistors (those for the first stage of the 
amplifier) are required to meet the initial noise figure specification limit 
of 4 dB at 5.5 MHz. Figure 7 shows that in fact nearly all of these devices 
still meet this limit after ten years of power aging. In the evolution and 
refinement of techniques for measuring the noise figure over the ten-year 
period of extended device aging, it was more important to provide the 
best and most realistic noise figure data than to maintain reference to 
the past. Other experience in aging indicates that if the other device 
parameters remain well-behaved, the noise figure does not change. 
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Fig. 12—Low-frequency common-base current gain for L2288 transistors—before and 
after power aging at 0.99 W. 


4.3 L2288 transistors 


The L2288 transistors, coded as 38C and operating at a power level 
of 0.99 W, are used in the output stage of the amplifier in the SF repeater. 
During aging, all of these devices are held at this same level of power, 
composed of 11 V and 90 mA. The extended aging program has included 
one hundred of these devices. 

Figure 8 shows that the collector-base breakdown voltage of these 
transistors has very uniformly shifted downward by about a quarter volt 
after ten years of power aging. This is the kind of result which arouses 
suspicion that the apparent shift may actually be due to an offset be- 
tween tests conducted over ten years apart. However, since the fact of 
testing offset could not be definitely established, the shift must be ac- 
cepted at face value. Even if it is real, such a shift is of no practical sig- 
nificance. Under a simple linear extrapolation the voltage would be de- 
creasing at the rate of 1 V in forty years, and a decrease of at least 30 V 
would be required before even the first device would begin to impair 
system performance. 
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Fig. 13—Second-harmonic distortion in dB below fundamental for L.2288 transistors— 
before and after power aging at 0.99 W 


In Fig. 9 the lower values of collector-base leakage current (below 
about 4 uA) show no discernible change after ten years of power aging. 
Those above about 4 nA show, if anything, a very slight decrease. 

The emitter-base breakdown voltages shown in Fig. 10 indicate very 
consistent behavior, with the possible exception of one device whose 
voltage decreased by about 0.1 V and another whose voltage increased 
by about 0.05 V, again with no significant effect on circuit or system 
reliability since the emitters are always forward biased. 

The emitter-base forward voltages shown in Fig. 11 indicate a small 
but consistent decrease of about 1 mV after ten years of power aging. 
Such a decrease is not significant enough to argue about whether it is 
real or an offset in testing. 

The behavior of the low-frequency common-base current gain shown 
in Fig. 12 is again very gratifying. An increase in this parameter was 
expected, but the magnitude of increase has turned out to be smaller than 
expected. If the first ten years of aging serve as any indication, the end 
of life limit as an average for all devices will probably never be reached 
for all practical purposes. 

Harmonic distortion for these devices was determined by measuring 
the level of the second and third harmonic at the output with a sine wave 
input at a frequency of 1.0 MHz at a power level of 5.0 dBm. Measure- 
ment precision is not in the same class with that of the dc and low-fre- 
quency parameters, especially for the initial data. During the ten-year 
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Fig. 14—Third-harmonic distortion in dB below fundamental for L2288 transistors— 
before and after power aging at 0.99 W. 


period of extended aging, the precision of measurement was improved 
considerably. Nevertheless, the aging results shown in Figs. 13 and 14 
should be considered merely indicative that the distortion parameters 
are quite well behaved and that after ten years of power aging, they still 
meet all of the device initial requirements. 


4.4 L2317 diodes 


The L2317 diodes, coded as 467A, are large area pn junction protective 
devices used in the power path of the SF repeater. During normal system 
operation, each device merely draws reverse-leakage current. In the event 
of a cable fault, however, it would conduct a large transient of reverse 
current to protect the repeater circuit from excess voltage. These devices 
are aged with 13-V reverse bias, which is a few volts below breakdown. 
Although the ten-year aging results are again presented in the form of 
initial and final values for each parameter, a large number of interme- 
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Fig. 15—Breakdown voltage in volts at 1.0 mA for L2317 diodes—before and after 10 
years of power aging at 13 V reverse bias. 
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Fig. 16é—Reverse-leakage current in nA at 13 V for L2317 diodes—before and after 10 
years of aging at 13 V reverse bias. 
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Fig. 17—Forward voltage in volts at 500 mA for L2317 diodes—before and after 10 years 
of power aging at 13 V reverse bias. 
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Fig. 18—Forward voltage in volts at 200 mA for L2318 diodes—before and after 10 years 
of aging at 10 mA rms and 60 Hz. 
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Fig. 19—Breakdown voltage in volts at 2 mA for L2320 diodes—before and after 10 years 
of aging at 2 mA in breakdown region. 


diate data points taken during the ten years show no significant excur- 
sions away from the end points shown. 

Figure 15 shows the behavior of the breakdown voltage of one hundred 
L.2317 diodes across ten years of power aging. A very consistent upward 
shift of about 0.01 V is readily apparent but is too small to be of any 
practical significance. The behavior of the reverse leakage current shown 
in Fig. 16 also exhibits a small but very consistent upward shift of about 
1 nA, while Fig. 17 shows a consistent downward shift of about 2 mV in 
the forward voltage. 


4.5 L2318 diodes 


The L2318 diodes, coded as 468A and 468B, are used for surge pro- 
tection in the oscillator circuit of the SF repeater. Each device consists 
of two silicon pn junction diode chips mounted on a common metallic 
heat sink and oppositely poled so that only a forward characteristic is 
observed with an applied voltage of either polarity. During long-term 
power aging, a 60-Hz voltage is applied to make the diode conduct 10 mA 
rms in each direction. The forward voltage drop is measured at two 
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Fig. 20—Reverse leakage current in nA at 8 V for L2320 diodes—before and after 10 
years of aging at 2 mA in breakdown region. 


current levels in each direction during long-term power aging. 

Figure 18 shows the behavior of the forward voltage of one hundred 
L2318 diodes in one direction and at one current level across ten years 
of power aging. No discernible aging can be seen. The behavior of the 
other forward voltages is not shown because they look no different. 


4.6 L2320 diodes 


The L2320 diodes, coded as 467B, are silicon pn junction limiter de- 
vices used in the third stage of the amplifier of the SF repeater. Operation 
is at 2.0 mA in the breakdown region, and the devices are aged under 
these power conditions. 

Figure 19 shows the behavior of the breakdown voltage of one hundred 
1.2320 diodes across ten years of power aging. No discernible shift in this 
parameter is evident. 

In Fig. 20, a definite upward shift in the reverse-leakage currents can 
be seen. This is the one case where a significant amount of scatter is 
encountered in the magnitudes of the parameter shifts among the 
hundred devices. The data suggests that percent change rather than 
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Fig. 21—Forward voltage in volts at 200 mA for L2320 diodes—before and after 10 years 
of aging at 2 mA in breakdown region. 


incremental (delta) current change is probably the appropriate way to 
regard the behavior of individual devices. A leakage-current magnitude 
of at least 100 nA would be required in any one device to produce even 
the start of any impairment in circuit operation. 

For this particular situation, even the aging of 100 devices for ten years 
still leaves some degree of uncertainty in how best to interpret the be- 
havior in the tails of the distribution where the device would be found 
with potentially the greatest change in leakage current. In consideration 
of this situation and in consideration of the fact that a long cable might 
typically have 400 of these devices compared to the 100 shown here, it 
has been conservatively estimated? that with high (90 percent) confi- 
dence, a 100-uA leakage current would not be reached in any one device 
in less than 60 years.° 

According to Fig. 21, the forward voltages of these diodes have shown 
a small but very consistent downward shift of about 1 mV in ten years 
of power aging. Such a shift is of no significance to system reliability. 
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V. GENERAL COMMENTS 


The regular procedure for producing submarine cable semiconductor 
devices can be divided into three general stages: manufacturing, 
screening, and aging. A selection process is employed in each of these 
stages. The devices in this extended aging program progressed through 
the first two stages in the usual way, but did not have the benefit of 
having completed the standard aging stage of about six months with 
ensuing selection before beginning the extended aging period. Never- 
theless, any of these devices, even after ten years of power aging, could 
still be considered candidates for cable use. 

For the next generation of submarine cable (SG), the first of which was 
laid across the Atlantic Ocean in the early part of 1976, silicon transistors 
have been used throughout. Probably no new systems will ever be de- 
signed with germanium transistors. Nevertheless, the reliability record 
thus far of the germanium transistors used in the SF system is rather 
impressive in its own right. 
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A “DOGLEG” Channel Router. D. N. Deutsch, Des. Automation Conf. Proceed., No. 
13 (June 1976), pp. 425-433. A new algorithm for interconnecting two sets of terminals 
across an intervening channel is presented. The routing is done on two distinct layers, one 
for horizontal paths, the other vertical. Typical channels (300 terminals, 100 nets) are 
routed within one track of the mathematical bound (density) in a few seconds. 


Impact of On-Line Systems on a Literature Searching Service. D. T. Hawkins, Spec. 
Libr., 67 (December 1976), pp. 559-567. Some experiences with on-line literature 
searching systems in the Bell Laboratories Libraries and Information Systems Center 
during 1975 are described. A total of 604 sessions, averaging 15 minutes each, occurred. 
Average cost per session was approximately $22, including connect time, communications, 
and off-line printing. The actual cost of many sessions, however, was $10 or less. On the 
average, on-line searching costs $1/minute and 10¢/off-line print. A search on a specific 
subject took an average of two sessions. Details on types of searches performed, the method 
of charging costs to the user, and some effects on library services are given, as well as some 
of the advantages of using an on-line system for information retrieval. 


LTX—A System for the Directed Automatic Design of LSI Circuits. G. Persky, D. 
N. Deutsch, and D. G. Schweikert, Des. Automation Conf. Proceed., No. 13 (June 1976), 
pp. 399-407. LTX is a minicomputer-based design system for large-scale integrated 
circuit chip layout which offers a flexible set of interactive and automatic procedures for 
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encompasses algorithms for two-dimensional placement, string placement, and channel 
routing. 


Pro—An Automatic String Placement Program for Polycell Layout. G. Persky, Des. 
Automation Conf. Proceed., No. 13 (June 1976), pp. 417-424. PRO is a minicompu- 
ter-based string placement program that reorders cells in polycell rows such that channel 
routing in the adjacent channels needs fewer routing tracks, and is not blocked by cyclic 
constraints. PRO utilizes cell reflections and pairwise neighbor exchanges, simultaneously 
monitoring both channels adjacent to the row being reordered. 


ELECTRICAL AND ELECTRONIC ENGINEERING 


Birefringent Coupler for Integrated Optics: Comment 2. A. Albanese and J. A. Arnaud, 
Appl. Opt., 15 (September 1976), p. 2029. A proof is given of the impossibility to realize 
lossless and one-directional couplers. 


Control of Limit Cycles in Recursive Digital Filters by Randomized Quantization. 
R. B. Kieburtz, V. B. Lawrence, and K. V. Mina, Proc. 1976 IEEE Int. Symp. Circuits 
Systems, April 1976, pp. 624-627. The different types of limit cycles that can occur 
in single second-order sections respond differently to efforts to reduce them. Length 1 
(dc) limit cycles play an important complicating role, especially in cascaded filters. We 
describe a simple method, involving random requantization of multiplier outputs, which 
can reduce or eliminate limit cycles in digital filters. We compare limit cycle and roundoff 
noise in three typical low-pass filters with and without the method. 


Impurity and Substrate Diffusion in a Thin Contact Layer. R. P. Goel and F. E. Bader, 
J. Electrochem. Soc., 123, No. 8 (August 1976), pp. 1242-1245. Based on simplifying 
assumptions, diffusion phenomena of substrate and impurity through a thin contact layer 
have been compared. Results are presented in terms of nondimensionalized parameters 
k (a measure of diffusion rate, time, and contact thickness) and 1/] (film thickness/contact 
thickness). A single master curve is shown to represent each diffusion process. 
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Low-Noise GaAs MESFETS. B. S. Hewitt, et al., Electron. Lett., 172, No. 12 (June 10, 
1976), pp. 309-310. GaAs MESFETs with optimum noise figures of 1.6 dB at 6 GHz 
have been fabricated by projection photolithography. An equation has been developed 
for the calculation of optimum noise figure which gives good agreement between calculated 
and measured values. 


Measurement of Land to Plated-Through-Hole Interface Resistance in Multilayer 
Boards, J. N. Hines, 14th Annual Proceed., 1976 Reliability Phys. Symp., April 20, 1976, 
pp. 132-134. Thermally induced cracks in the bond between an innerlayer land and 
barrel will occur during the lifetime of an MLB if processing is faulty. A technique to detect 
incipient failures that features four-terminal resistance measurements with a minimum 
of ey resistance to overcome the insensitivity of series circuit measurements is pre- 
sented. 


Mechanical-Chemical Technique for Removal of Epitaxial Spikes. L. E. Katz and 
W.C. Erdman, J. Electrochem. Soc., 123, No. 8 (August 1976), pp. 1249-1251. Anew 
mechanical-chemical spike removal technique has been devised which effectively removes 
epitaxial spikes on silicon substrates. A mask damage technique was used to evaluate masks 
subjected to spike removal. Results show effective spike removal as well at no subsequent 
mask damage. 


Millimeter Wave Integrated Circuit Technologies. R. H. Knerr and L. F. Moose, 
Microwave J., 19, No. 11 (November 1976), pp. 23-24. Introduction to a two part series 
of short articles authored by the participants in a panel discussion on “Millimeter Wave 
Integrated Circuit Technologies” at the 1976 IEEE-MTT-S International Microwave 
Symposium. 


A New Digital Echo Canceler for Two-Wire Full-Duplex Data Transmission. K. 
H. Mueller, IREE Trans. Commun., COM-24, No. 9 (September 1976), pp. 956-962. A 
new approach to echo canceling for two-wire full-duplex data transmission is proposed. 
The canceling signal is directly synthesized from the binary data, using a transversal filter 
approach, and the usual multiplications are replaced by additions and subtractions, thus 
allowing efficient operation of a large number of taps as required for the canceling of distant 
echoes. As a specific application, a system processing one sample per baud is discussed 
where timing signals at both communicating stations are assumed to be synchronized. 


Optimization of Digital Postdetection Filters for PSK Differential Detectors. R. 
D. Gitlin and K. H. Mueller, IEEE Trans. Commun., COM-24, No. 9 (September 1976), 
pp. 963-970. A digital transversal postdetection filter for use in differential PSK de- 
tectors is proposed, and a method for optimizing the weighting coefficients is presented. 
It is shown that even though the detector output signal is quadratic, the optimization 
pene can be formulated in a way similar to the well-known MMSE equalizer for linear 
signals. 


MATERIALS SCIENCE 


Application of Weibull-type Analysis to the Strength of Optical Fibers. B. K. Tariyal 
and David Kalish, Mater. Sci. Eng., 27 (January 1977), pp. 69-71. The strength dis- 
tribution for 400 measurements on 0.61-m-length fused silica optical fiber specimens ex- 
hibits a bimodal behavior. The cumulative failure probability is presented for Weibull-type 
distributions and is analyzed by assuming that the strength distribution is comprised of 
two distinct populations having finite lower and upper strength limits. This procedure 
provides the best prediction of the strength in long lengths. 


Many-Body Effects in Transition Metals: Role of the Density of States. G. K. Wer- 
theim and L. R. Walker, J. Phys. F Metal Phys., 6, No. 12 (December 1976), pp. 2297- 
2306. The effects of structure in the density of electron-hole pair excitations is ex- 
amined by numerical integration of the many-body screening formalism. The results are 
applied to Ir and Pt. 


Oxidation State of Tungsten in the Na, WO3 Bronzes. G. K. Wertheim, M. Campagna, 
J.-N. Chazalviel,* H. R. Shankst, Chem. Phys. Lett., 44, No. 1 (November 15, 1976) pp. 
50-52. X-ray photoemission spectra of W 4f electrons in vacuum-cleaved cubic 
Na, WOs3 covering the entire composition range give no evidence for the existence of 
multiple W valence states in the bulk. Strong oxidation effects are observed in air exposed 
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surfaces. * Laboratoire PMC, Ecole Polytechnique, France. f Ames Laboratory, Iowa 
State Univ. 


Thin Film Interaction in Aluminum and Platinum. S. P. Murarka, I. A. Blech, and 
H. J. Levinstein, J. Appl. Phys., 47 (December 1976), pp. 5175-5181. Aluminum and 
Platinum interacted very rapidly with the formation of several intermetallics (namely 
Pt3Ale, PtsAls, PtAle, PtAls and PtAl4) which resulted in large volume changes and loss 
of adhesion. The Al-Pt interaction rate was dependent on the annealing ambient being 
higher in forming gas, argon or helium than in vacuum or air. 


MECHANICAL AND CIVIL ENGINEERING 


Free Vibrations of a Beam-Mass System With Elastically Restrained Ends. R. P. 
Goel, J. Sound Vib., 47, No. 1 (July 1976), pp. 9-14. The vibration problem of a beam 
with an arbitrarily placed concentrated mass and elastically restrained against rotation 
at either end is solved. The effects on eigenfrequencies of the system produced by varying 
the mass, stiffness, and position ratios are presented. 


Tranverse Vibrations of Tapered Beams. R. P. Goel, J. Sound Vib., 47, No. 1 (July 
1976), pp. 1-7. Transverse vibrations of linearly tapered beams, elastically restrained 
against rotation at either end, have been investigated. Results for the first three eigen- 
frequencies, with different values of stiffness ratios and taper ratios are presented. Cases 
of a tapered cantilever beam with a concentrated mass at the free end and spring hinged 
at the other end have also been presented. 


PHYSICS 

Absence of the Low-Temperature Specific-Heat Anomaly in bee *He. D. S. Greywall, 
37, No. 2 (July 12, 1976), pp. 105-107. The specific heat of bec “He has been measured 
along three isochores for temperatures T between 50 mK and the melting curve. The data 
for T < 0.3 K can be described well by the function C, = «aT? + yT* and are inconsistent 
os the existence of the low-temperature specific-heat anomaly reported previously by 
others. 


Digital Normalization of Iodine Filter Structure in Quasielastic Light Scattering. 
K. B. Lyons and P. A. Fleury, Phys. Rev. Lett., 47 (November 1976), pp. 4898-4900. A 
molecular iodine filter may be used to reduce the elastic component in a light scattering 
experiment but also introduces extraneous structure into the spectrum. A technique is 
eepored here to compensate quantitatively for this structure by a digital normalization 
scheme. 


Enhanced Optical Fluorescence by Resonant Méssbauer Excitation. C. P. 
Lichtenwalner, H. J. Guggenheim, and L. Pfeiffer, Phys. Lett., 564A, No. 2 (March 8, 1976), 
pp. 117-118. Enhanced optical Eu?* fluorescence in Eu-doped CaF» crystals is ob- 
served following resonant absorption of Méssbauer y-rays by Eu** nuclei. The effect is 
due to scintillation fluorescence from Eu?t ions excited by conversion electrons from 
decaying 21.6 keV 15!Eu. 


Particle Brownian Motion Near a Hydrodynamic Instability. J. B. Lastovka and J.-P. 
Boon,* Phys. Rev., 14, No. 4 (October 1976), pp. 1583-1586. We consider the alteration 
of the normal Brownian motion of seed particle brought about by the enhanced velocity 
fluctuations that exist near the threshold of a hydrodynamic instability. We utilize our 
results to quantitatively predict the fluctuation induced effects for the threshold region 
of the Bénard convective instability. * Université Libre de Bruxelles. 


Pressure-Induced Valence Change in Cerium Phosphide. A. Jayaraman, W. Lowe, 
L. D. Longinotti, and E. Bucher. The pressure-volume relationship for cerium mo- 
nophosphide (CeP) has been determined to 200 kbars pressure. From the data, we conclude 
that there is an electronic transition involving a change in the valence state of Ce from 3+ 
towards the 4+ state near 100 kbar pressure. Qualitative observations of the reflectivity 
of CeP under pressure support the valence transition. 


ABSTRACTS 1011 


The Use of Our Specially Designed Nonachromatic Microscope Objective to Ex- 
amine Tracks in Nuclear Emulsions. J.S. Courtney-Pratt, Opt. Eng., 15, No. 4 (July- 
August 1976), pp. 379-383. We have shown that our specially designed nonachromatic 
objective allows us to produce sharp images of particle tracks in nuclear emulsions, for 
visual study and/or photographic recording, with a lateral resolution of about one micron, 
and a depth of field many times greater than can be achieved with conventional objec- 
tives. 
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B.S.T.J. Briefs 


Adhesive Sandwich Optical Fiber Ribbons 
By M. J. SAUNDERS and W. L. PARHAM 
(Manuscript received April 15, 1977) 


In 1974 a proposal was made to put optical fibers together into easily 
handled units for optical communication purposes.! This proposal 
suggested “the use of fiber ribbons consisting of linear arrays of fibers 
embedded in a thin, flexible supporting medium as components of a cable 
for fiber transmission systems.” This note is a brief description of fiber 
ribbons made by sandwiching fibers between two layers of polyester- 
backed adhesive (adhesive sandwich ribbons). 

The machine for making these ribbons has evolved around two 
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Fig. 1—Schematic diagram of adhesive sandwich ribbon machine. 
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motor-driven pressure rollers. Figure 1 is a schematic diagram of the 
machine. Twelve fibers from payout reels are directed through the fiber 
guide where they are made contiguous. The fibers then pass between the 
two pressure rollers and are sandwiched between two layers of polyester 
backed adhesive tape supplied from large reels. After emerging from the 
pressure rollers, the sandwich is cut by a blade cutter to a width some- 
what larger than the width of the 12 fibers. The excess tape that has been 
cut away from the fiber-carrying ribbon is attached to a reel that supplies 
the power to pull the tape through the system. The ribbon with the fibers 
is either wound on a reel or is permitted to fall, under zero tension, into 
a container so that loss measurements can be made. Currently, we are 
making kilometer lengths of ribbon at speeds of about 0.2 m/s. However, 
shorter lengths of ribbon have been made at a speed of 0.8 m/s. We have 
made about 100 adhesive sandwich ribbons, including ribbons for the 
Atlanta Fiberguide System Experiment,’ varying in length from 0.2 to 
2 km. 

The added loss of fibers in these ribbons, caused by microbending,® 
is quite small. The average added loss of the fibers in the ribbons used 
to make the two cables, each 1 km in length, for the Atlanta Fiberguide 
System Experiment was 0.95 + 0.07 dB/km for 137 fibers and 0.55 + 0.05 
dB/km for 132 fibers.4 Some fiber breakage occurred, caused both by 
defects in the machine and by the fibers being wound improperly on the 
reels. 

We thank M. J. Buckler for making the loss measurements, E. D. Knab 
for much of the design work, and R. Sabia for originally suggesting the 
concept of the adhesive sandwich ribbon. 
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